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SECTION  1 


INTRODUCTION 

The  deposition  and  characterization  of  thin  films  has  received  enormous  attention  owing  to 
their  importance  in  the  fabrication  of  microelectronic  devices  and  solar  cells.  Considerable 
attenti'.n  has  been  given  to  silicon  films  deposited  on  quartz  or  other  substrates  and,  more  recently, 
to  in  situ  silicon  doping  of  compound  semiconductors  such  as  GaAs.  One  process  for  film 
deposition  is  to  induce  a  gas  phase  reaction  by  either  pyrolytic  or  photolytic  means,  the  film  growth 
being  determined  mainly  by  gas  phase  decomposition  kinetics.  Such  chemical  vapor  deposition 
(CVD)  processes  are  often  slow  and  inefficient,  particularly  if  laser  irradiation  is  required  to  induce 
the  decomposition  reactions. 

Another  route  toward  thin  film  production  involves  the  plasma  enhanced  deposition  of  silicon 
(Reference  1)  from  a  gas  mi.xture  containing  a  silicon- bearing  molecule  such  as  silane  (SiH4), 
disilane  (312^6)  ^  halogen  substituted  silicon  compound  such  as  silicon  tetrafluoride  (SiF4).  In 

particular,  glow  discharges  of  silane  or  silane-noble  gas  mixtures  have  attracted  considerable 
interest  as  a  means  of  producing  device -quality  hydrogenated  amorphous  silicon  (a-Si.H)  films 
(Reference  2).  The  growth  of  such  films  depends  on  the  complicated  chemistry'  that  is  occurring  in 
the  cathode  sheath  region  of  the  typical  capacitively  coupled  rf  discharge.  A  good  deal  of  progress 
has  been  made  over  the  past  few  years  in  understanding  the  electron  energy  distribution  and  the 
primary  collision  processes  that  occur  in  the  sheath  region  (References  3  •  5).  However,  a 
quantitative  description  of  all  the  major  reactions  that  occur  in  the  discharge  region,  including 
lectron  impact  dissociation  and  ionization,  and  the  chemical  role  of  both  cations  and  anions 
(formed  by  dissociative  attachment)  is  required  to  permit  reliable  modeling  of  rf  discharge  sheaths. 
A  goal  of  this  research  program  was  to  identify  the  ion  chemistry  occurring  in  silane  discharges, 
and  to  predict  both  the  kinetic  reaction  rates  and  product  distributions.  Emphasis  was  placed  on  the 
development  of  sound  theoretical  methods  to  predict  the  collisional  energy  dependence  of  the 
reaction  rates  and  to  examine  the  role  of  dissociative  charge  transfer  in  silane  ion-molecule 
reactions. 

Several  qualitative  pictures  of  both  ion  and  neutral  fragment  production  in  silane  discharges 
have  been  developed.  Haller  (Reference  6)  has  analyzed  the  current  status  of  ion  species  of  the 
type,  Si^Hn,"*",  using  mass  spcctrometric  techniques.  Secondary  ion-silane  reactions  to  form 
Jisilane  ions  were  proposed  to  explain  the  observed  ion  currents.  A  separate  model  was  proposed 
by  Turban,  et  al  (Reference  7)  in  which  primary  silane  ions  arc  formed  by  low-energy  electron 
impact  dissociation,  followed  by  secondary  ton  molecule  reactions.  Chapman  and  Gallagher 
(Reference  5),  however,  have  criticized  the  applicability  of  this  model  to  cathode  sheath  conditions 
where  it  is  believed  that  high  energy  electron  processes  (10- 100  eV)  govern  the  overall  reaction 
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kineticb.  Their  interpretation  of  the  discharge  ion  chemistry  is  that  electron  impact  dissociation  of 
silane  occurs  as  a  primary  process  forming  SiHn+  (n=0,3),  with  n=2  or  3  heavily  favored  over 
more  highly  stripped  silane.  Neutral  SiH^  radical  production  as  a  primary  process  must  also  be 
considered.  Based  on  the  known  rate  of  H  atom  abstraction  from  silane  and  the  direct  electron 
impact  dissociation  reaction,  SiH3  radicals  are  believed  to  predominate  in  the  discharge  region. 
SiH2''^,  and  more  highly  stripped  silane  cations,  are  probably  driven  to  SiH3  by  fast  reactions  of 
the  type  SiHn"^  +  SiH4  SiH^+i  +  SiH3‘^. 

The  largest  uncertainties,  at  present,  in  the  chemistry  of  the  discharge  sheath  region,  are  with 
the  secondary  ion-molecule  reactions  occurring  in  pure  silane-noble  gas  mixtures.  The  possible 
silane  ion- molecule  reactions  are  well  known  (Reference  8)  but  in  most  cases,  the  rate  coefficients 
.ne  known  experimentally  only  in  the  low-energy  region  (<  1  eV)  or  have  been  estimated  using  the 
Langevin-Gioumousis-Stevenson  formula  for  the  collision  cross  sections.  This  simple  model 
predicts  that  the  ion-molecule  reaction  rate  coefficients  are  independent  of  collision  energy. 
Dissociative  charge  traiisfer  rates  for  collisions  between  SiHn'*"  ions  and  silane  are  mainly 
unknown  at  present  since  knowledge  of  the  detailed  mechanism  of  the  ion-molecule  collision  is 
required  t  o  determine  branching  into  the  dissociation  channels.  The  somewhat  simpler  reactions 
between  silane  and  H+,  H3'^  or  the  noble  gas  ions,  which  are  more  amenable  to  theoretical 
analysis,  are  also  mainly  unknown  at  present  in  a  mechanistic  sense.  In  addition,  a  preliminary 
analysis  of  the  role  of  silicon  hydride  anions,  SiHn’,  has  recently  been  reported  (Reference  4). 

This  study  suggested  that  electron  dissociative  attachment  to  silane  may  exhibit  a  large  cross 
section,  particularly  if  the  silane  molecule  is  vibratior.ally  or  rotationally  excited.  The  role  of  anion 
chemistry  in  describing  the  kinetic  processes  in  the  sheath  region  is  uncertain  at  present.  As  an 
example,  Garscadden  (Reference  9)  and  Haaland  (Reference  10)  have  dismissed  the  possible  role 
of  silane  anions  in  plasma  assisted  CVD  devices.  A  significant  production  rate  for  SiHn"  would 
greatly  influence  our  understanding  of  the  overall  chemistry  since  fast  mutual  neutralization 
reactions  between  cations  and  anions  could  occur. 

The  intent  of  this  technical  program  was  to  address  these  uncertainties  I.i  the  chemistry  of 
silane  ion-molecule  reactions  through  a  theoretical  analysis  of  the  governing  reaction  rates.  Since 
dissociative  reactions,  with  or  without  accompanying  charge  transfer  may  occur,  our  approach  was 
based  on  detailed  quantum  mechanical  studies  of  the  reaction  surfaces  to  yield  all  possible  reaction 
product  channels.  Subsequent  to  this  analysis  of  the  reaction  energetics,  cross  sections  for  the 
various  ion-molecule  processes  were  calculated  within  a  quantum  framework.  During  the  course 
of  this  research  investigation,  several  new  areas  were  examined  at  an  exploratory  level.  These 
areas  mainly  involved  the  formation  and  reactions  of  negatively  charged  species  whose  chemistry 
and  relative  importance  in  silane  plasmas  is  uncertain.  In  contrast  to  positive  ions,  which  are  swept 
to  the  walls  by  the  ambipolar  fields  in  the  discharge,  negative  ions  tend  to  become  trapped  with 
resultant  long  residence  times.  Garscadden  (Reference  1)  has  suggested  that  negative  ions  may 
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form  clusters  through  a  series  of  ion-molecule  reactions  similar  to  that  discussed  by  Mandich, 
Reents  and  Jarrold  (Reference  1 1)  and  Raghavachari  (Reference  12),  for  positive  ion  cluster 
growth.  It  is  not  clear  whether  sequential  clustering  reactions  of  negative  silane  ions  will  encounter 
a  bottleneck  at  a  iinall  cluster  size  or  whether  the  anion  chemistry'  is  sufficiently  different  that 
growth  continues  to  the  critical  sizes  needed  for  spontaneous  formation  of  large  silicon  particulates. 
Additional  research  effort  in  silane  negative  ion  chemistry  is  recommended. 

The  general  composition  of  this  report  is  as  follows.  In  Section  2,  we  present  a  description 
of  the  mathematical  methods  which  were  employed  in  this  research.  Included  in  Section  2  are  sub¬ 
sections  which  deal  with  the  construction  of  electronic  wavefunctions,  the  calculations  of 
expectation  properties,  the  calculation  of  electronic  wavefunctions  using  ab  initio  methods,  a 
description  of  the  development  of  R-matrix  codes  and  a  description  of  several  electronic  structure 
codes  that  were  utilized  during  the  course  of  this  research  program.  The  calculated  results  and 
pertinent  discussion  are  presented  in  Section  3.  Based  on  these  present  studies,  we  have  identified 
several  areas  for  further  technical  work.  Recommendations  are  also  presented  in  Section  3  and  a 
concluding  summary  of  this  program  is  given  in  Section  4.  A  list  of  all  technical  presentations  and 
publications  that  have  resulted  from  this  research  program  are  given  in  Appendix  A.  Abstracts  of 
the  pertinent  published  papers,  reports  and  technical  talks  are  collected  in  Appendix  B.  Appendix 
C  contains  a  preprint  of  a  paper  to  be  submitted  for  publication. 
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SECTION  2 


MATHEMATICAL  BACKGROUND 

Central  to  these  theoretical  studies  are  the  actual  quantum-mechanical  calculations  which 
must  be  carried  out  for  the  atomic  and  molecular  species.  In  particular,  the  methodology  for 
computing  the  minimum  energy  reaction  pathways  and  testing  for  vibrational  stability  must  be 
carefully  analyzed.  For  added  clarity,  various  aspects  of  these  calculations  are  discussed  below. 

Much  evidence  on  diatomic  and  polyatomic  systems  indicates  the  inadequacy  of  minimum 
basis  sets  for  constructing  quantitatively  correct  molecular  wavefunctions  (References  13  and  14). 
This  means  inner-shell  and  valence-shell  orbitals  of  quantum  numbers  appropriate  to  the  atoms 
(Is,  2s,  2p,  for  C,  N,  O;  etc.).  The  main  deficiency  of  the  minimum  basis  set  is  its  inability  to 
properly  describe  polarization  and  the  change  of  orbital  shape  for  systems  which  exhibit  large 
charge  transfer  effects.  V'alues  of  the  screening  parameters  for  each  orbital  can  either  be  set  from 
atomic  studies  or  optimized  in  the  molecule,  the  latter  approach  is  indicated  for  studies  of  higher 
precision.  When  high  chemical  accuracy  is  required,  as  for  detailed  studies  of  the  ground  or  a 
particular  excited  state  of  a  system,  a  more  extended  basis  must  be  used.  Double-zeta  plus 
polarization  functions  or  optimized  MO’s  are  usually  required  as  a  minimum  representation  for 
reliable  calculated  results  of  chemical  accuracy. 

The  chosen  basis  sets  give  good  results  only  when  used  in  a  maximally  flexible  manner.  This 
implies  the  construction  of  perturbation  expansions  or  the  use  of  Cl  wavefunctions  with  all  kinds 
of  possible  orbital  occupancies,  so  that  the  correlation  of  electrons  into  overall  states  can  adjust  to 
an  optimum  form  at  each  geometrical  conformation  and  for  each  state.  Except  when  well-defined 
pairings  exist,  as  for  closed  shell  and  exchange  dominated  systems,  a  single-configuration  study 
(even  of  Hartree-Fock  quality)  will  be  inadequate. 

Proper  electronic  states  for  systems  composed  of  light  atoms  should  possess  definite 
eigem  alues  of  the  spin  operator  as  well  as  an  appropriate  geometrical  symmetry.  The  geometri¬ 
cal  symmetry  can  be  controlled  by  the  assignment  of  orbitals  to  each  configuration,  but  the  spin 
state  must  be  obtained  by  a  constructive  or  projective  technique.  Formulas  have  been  developed 
(Reference  15)  for  projected  construction  of  spin  states  from  orthogonal  orbitals,  and  programs 
implementing  these  formulas  have  been  in  routine  use  at  UTRC  for  several  years.  One  of  the  least 
widely  appreciated  aspects  of  the  spin-projection  problem  is  that  the  same  set  of  occupied  spatial 
orbitals  can  sometimes  be  coupled  to  give  mure  than  one  overall  state  of  a  given  spin  quantum 
number.  It  is  necessary  to  include  in  calculations  all  such  spin  couplings,  as  the  optimum  coupling 
will  continuously  change  with  changes  in  the  molecular  conformation.  This  is  especially  important 
in  describing  degenerate  or  near-degenerate  excited  electronic  states. 


4 


In  the  sections  below,  we  describe  the  several  mathematical  approaches  that  are  applicable 
to  calculation  of  the  electronic  structure  of  molecules,  to  the  calculation  of  potential  energy 
surfaces  for  chemically  reacting  systems,  and  to  the  subsequent  calculation  of  radiation  and 
collisional  processes.  Since  several  different  approaches  are  indicated,  we  describe  their  expected 
areas  of  applicability. 

2.1  Method  of  Ab  Initio  Calculation 

2.1.1  Born-Opponheimer  Separation 

For  a  system  of  n  electrons  and  N  nuclei,  and  considering  only  electrostatic  interactions 
between  the  particles,  we  have  for  the  total  Hamiltonian 


36  =  36el-Xr— Vi  + 


2Mt 


N  N 


N  n 


I  Iv<,  v^  +  222;v,-v,+  XX  ’i'^i 
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and  where  mg,  me ,  Mj,  are  the  masses  of  the  electron,  atom  and  combined  system  mass,  respec¬ 
tively.  Now  since  the  ratios  mg/ma ,  and  me/Mj  are  both  small,  (2- 10"^  -  5  •  lO"'*)  we  can  effect 

a  separation  of  the  electronic  and  nuclear  coordinates  treating  the  total  wavefunction  as  a  product 
of  a  nuclear  and  an  electronic  part.  We  have 


V'Cfn.RN)  =  Xil^k(RNVk(fn,RN)  (3) 

k 

where  Rn)  is  an  electronic  wavefunction  parametric  in  the  nuclear  coordinates  as  given  in 

Equation  (3)  and  ji'kCRN)  are  nuclear  motion  wavefunctions  which  satisfy  (neglecting  terms  of  the 
order  of  mg/  ma ) 
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The  cross  term  in  Vq  can  be  eliminated  by  a  proper  change  of  variables  and  Equation 
(4)  then  reduces  to  a  3N-3  dimensional  Schrddinger  equation. 

For  most  systems,  where  the  velocity  of  motion  of  the  nuclei  is  slow  relative  to  the  electron 
velocity,  this  decoupling  of  electron  and  nuclear  motion  is  valid  and  is  referred  to  as  the  adiabatic 

approximation.  Equation  (3)  thus  defines  an  electronic  eigenstate  Rn).  parametric  in  the 

nuclear  coordinates,  and  a  corresponding  eigenvalue  Ek(RN)  which  is  taken  to  represent  the 
potential  energy  curve  or  surface  corresponding  to  state  k. 

2.1.2  Variational  Methods 

By  an  ab  initio  method  is  meant  one  that  starts  from  a  zero-order  Hamiltonian  which  is 
exact  except  for  relativistic  and  magnetic  effects,  and  which  involves  the  evaluation  of  electronic 
energies  and  other  relevant  quantities  for  wavefunctions  which  are  properly  antisymmetrized  in 
the  coordinates  of  all  the  electrons.  For  a  system  containing  n  electrons  and  M  nuclei,  the 
zero-order  Hamiltonian  depends  parametrically  on  the  nuclear  positions  and  is  of  the  form 


1  N  n  M  M 

a6  =  -iyv?-y  y-^+  y  y  ^ 

2  ifl  ifl  jfl  Fi  -  Rj  1  ^r<j  Ri  -  Rj  1  ii<  j  h  -  Tj 


(5) 


where  Z,  and  R,  are  the  charge  and  position  of  nucleus  i,  Fj  is  the  position  of  electron  j,  and  Vj^ 

is  the  Laplacian  operator  for  electron  j.  All  quantities  are  in  atomic  units,  i.e.  lengths  in  bohrs, 
energies  in  hartrees  (1  hartree  =  2  Rydbergs).  The  many-electron  wavefunction  consists  of  one,  or 
a  linear  combination,  tp  =  I.  of  terms  of  the  form 


i=I 

where  each  is  a  spatial  orbital,  A  is  the  antisymmetrizing  operator,  Oj  is  the  spin-projection 
operator  for  spin  quantum  number  S,  and  Of,tn  is  a  product  of  a  and  one-electron  spin 
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functions  of  magnetic  quantum  number  Mj.  No  requirement  is  imposed  as  to  the  double  occu¬ 
pancy  of  the  spatial  orbital,  so  linear  combinations  of  the  form  given  by  Equation  (6)  can  describe 
a  completely  general  wavefunction.  The  spatial  orbitals  0^,  may  be  whatever  basis  orbitals  have 

been  introduced,  arbitrary  linear  combinations  thereof,  or  specific  linear  combinations  deter¬ 
mined  pursuant  to  the  particular  calculational  method  in  use. 

The  spatial  orbitals  0^, ,  the  spin  functions  and  the  coefficients  of  different 
linear  combination  of  ^1)^  is  used,  may  be  explicitly  determined  by  invoking  the  variational 

principle.  Various  specific  methods  are  described  below  for  determining  wavefunctions.  However, 
we  should  first  observe  that  the  adequacy  of  an  ab  initio  calculation,  or  for  that  matter  any  energy 
calculation,  will  depend  crucially  upon  the  extent  to  which  the  wavefunction  can  be  qualitatively 
appropriate.  Some  of  the  considerations  surrounding  the  choice  of  a  wavefunction  are  the 
following: 

(i)  Necessity  that  the  wavefunction  possess  sufficient  flexibility  to  be  able  to  describe 
dissociation  to  the  correct  atomic  and  molecular  fragments  as  various  internuclear  separations  are 
increased; 

(ii)  Maintenance  of  equivalent  quality  of  calculation  for  nuclear  geometries  differing  in  the 
nature  or  number  of  chemical  bonds; 

(iii)  Ability  to  describe  degenerate  or  near-degenerate  electronic  states  when  they  are 
pertinent; 

(iv)  Ability  to  describe  different  electronic  states  to  equivalent  accuracy  v/hen  their  Interrela¬ 
tion  (e.g.,  crossing)  is  relevant,  in  particular,  ability  to  de.scribe  ionic-valence  state  mixing; 

(v)  Ability  to  represent  changes  in  the  coupling  of  electron  spins  as  bonds  are  broken  or 
reformed. 

The  foregoing  considerations  indicate  that  it  will  often  be  necessary  to  consider  wavefunc¬ 
tions  with  more  than  a  minimum  number  of  singly-occupied  spatial  orbitals,  and  that  there  will 
be  many  potential  curves  or  surfaces  for  which  a  wavefunction  consi.sting  of  a  single  %  cannot 

suffice.  It  will  then  be  necessary  to  allow  mixing  of  ipfi  with  different  degrees  of  orbital  spatial 

occupancy  so  as  to  obtain  smooth  transitions  ’"'oni  the  occupancies  characteristic  of  separated 
atoms  or  molecules  (or  ions)  to  those  characteristic  of  a  compound  system  or  a  different  fragmen¬ 
tation. 

Another  implication  of  the  considerations  surrounding  the  choice  of  a  wavefunction  is 
related  to  the  treatment  of  electron  spin.  Not  only  is  it  necessary  to  require  that  the  wavefunction 
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be  an  ei,'»enfunction  of  and  Sz  but  it  is  also  necessary  to  take  account  of  the  fact  that  under 
many  conditions,  there  will  be  more  than  one  spin  eigenfunction  of  given  S  and  ms.  The  different 
spin  eigenfunctions  correspond  to  different  couplings  among  the  individual  spins.  Since  reactive 
processes  involve  the  breaking  and  forming  of  electron-pair  bonds,  they  must  necessarily  be 
accompanied  by  reorganizations  of  the  spin  coupling.  A  failure  to  take  account  of  this  will  lead  to 
qualitatively  inappropriate  wavefunctions. 

In  Hartree-Fock  calculations  ip  is  restricted  to  a  single  which  is  assumed  to  consist 
as  nearly  as  possible  of  doubly-occupied  orbitals.  The  orbitals  are  then  selected  to  be  the 
linear  combinations  of  basis  orbitals  best  satisfying  the  variational  principle.  Writing 
(pftt  =  2  a^ccv .  the  av,  are  determined  by  solving  the  matrix  Hartree-Fock  equation 


F;ivai,i  €iY  SAvSyi  (each  A)  (7) 

V  V 


where  is  the  orbital  energy  of  4>fti  • 


The  Fock  operator  F,iv  has  been  thoroughly  discussed  in  the  literature  (Reference  16)  and  de¬ 
pends  upon  one-  and  two-electron  molecular  integrals  and  upon  the  .  This  makes  Equa¬ 
tion  (7)  nonlinear  and  it  is  therefore  solved  iteratively.  UTRC  has  developed  programs  for  solving 
Equation  (7)  for  both  closed  and  open-shell  systems,  using  basis  sets  consisting  of  either  Slater- 
type  or  Gaussian-type  atomic  orbitals.  Examples  of  their  use  are  in  the  literature  (Reference  17). 

In  configuration  interaction  calculations,  the  overall  wavefunction  has  more  than  one  term, 
ipf, ,  and  the  c^  are  determined  by  invoking  the  variational  principle  to  obtain  the  secular 

equation 


X(^/<*’“Ws,v)cv  =  0  (each  fi) 
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Equation  (8)  is  solved  by  matrix  diagonalization  using  either  a  modified  Givens  method  (Refer¬ 
ence  18)  or  a  method  due  to  S'havitt  (Reference  19)  or  Raffenetti  (Reference  20). 

The  matrix  elements  and  S;,v  may  be  leduced  by  appropriate  operator  algebra  to  the 
forms 

n'P4ri.RN)|3fi(RN)Pi  n'-yvi(r;.RN))  (10) 

P  i=l  i=l 

S;,v  =  I  /p  {Om  I  %):  \  ^ )  ( 1  ]  %i(ri,  Rn)  I  p  I  fl  I'viCri,  Rn)  )  (11) 

P  i=l  i=l 

where  P  is  a  permutation  and  €p  its  parity.  Tire  sum  is  over  all  permutations.  (9ni|OsPl0m)  is  a 

“Sanibei  coefficient”  and  the  remaining  factors  are  spatial  integrals  which  can  be  x^actured  into 
one-  and  two-electron  integrals.  If  the  0^,  ^re  orthonormal,  Equations  (10)  and  (11)  become  more 

tractable  and  the  and  may  be  evaluated  by  explicit  methods  given  in  the  literature  (Ref¬ 
erence  15).  Computer  programs  have  been  developed  for  carrying  out  this  procedure,  and  they 
have  been  used  for  problems  containing  up  to  106  total  electrons,  10  unpaired  electrons,  and  sev¬ 
eral  thousand  configurations. 

The  Cl  studies  described  above  can  be  carried  out  for  any  orthonormal  set  of  for  which 

the  molecular  integrals  can  be  calculated.  Programs  developed  by  UTRC  make  specific  provision 
for  the  choice  of  the  as  Slater-type  atomic  orbitals,  as  Gau.>sian-type  orbitals,  as  symmetry 

molecular  orbitals,  as  Hartree-Fock  orbitals,  or  as  more  arbitrary  combinations  of  atomic 
orbitals. 

The  one-  and  two-electron  integrals  needed  for  the  above  described  method  of  calculation 
are  evaluated  for  STO’s  by  methods  developed  by  this  Center  (Reference  21).  For  Gaussian 
orbitals,  either  the  Carnegie-Mellon  integral  package  (Reference  22)  or  the  integral  routines 
incorporated  in  the  HONDO  or  GAMESS  programs  (Reference  23)  can  be  employed. 

2.1.3  Configuration  Selection 

Using  a  double-zeta  plus  polarization  basis  set  of  on*. -electron  functions,  a  typical  system 
can  easily  have  of  the  order  of  10^  configurations  in  full  Cl  (that  resulting  from  all  possible  orbital 
occupancies).  It  is  therefore  essential  ti>  identify  and  use  the  configurations  describing  the 
significant  part  of  'he  wavefunction.  There  are  several  ways  to  accomplish  this  objective.  First, 


one  may  screen  atomic -orbital  occupancies  to  eliminate,  configurations  with  excessive  formal 
charge.  Alter"  rtively,  in  a  molecular-orbital  framework,  nne  may  eliminate  configurations  with 
excessive  numbers  of  anti-bonding  orbitals.  A  third  possibility  is  to  cany  out  an  initial  screening 
of  ccnfigurations,  rejecting  those  whose  diagonal  energies  and  mtei action  matrix  elements  do  not 
satisfy  energy  significance  criteria. 

Another  common  method  of  classifying  configurations  is  to  exami  «c  the  total  number  of 
orbitals  in  the  wavefunction  that  differ  from  the  SCF  reference  wavefunclion.  We  write  the  Cl 
.vavefunction  as 


'*'■■■  (12) 

«  j 

where  a  single  excitation  function,  0  f  differs  from  the  SCF  reference  %  by  one  orbital  and  the 
double  excitation  function  by  two,  and  so  on.  The  Cl  coefficients,  C,,  are  then  determined 

variationally  to  yield  the  lowest  possible  total  energy.  In  the  limit  of  a  complete  basis  set  (N-»oo), 
and  where  all  possible  substitutions  are  included  in  0ci  >  the  variational  energy  approaches  the 
correct  nonrelativistic  Bor"  Oppenheimer  result.  Errors  arise  from  a  truncation  of  the  functions 
used  to  determine  the  SCF  refeience  wavefunction  and  from  the  i.’.ucation  of  the  excitation  func 
tions  series  in  0ci  •  The  reference  wavefunction  %  will  typically  be  the  .„me  for  a  Cl  or  Many 
Body  Perturbation  Theory  (MBPT)  calculation,  however,  for  a  molecule  even  as  small  as  water, 
0C1  becomes  a  function  of  a  very  large  number  of  basis  functions.  Additionally,  the  refers.';. e 

state  itself  may  be  multi-dimensional  in  order  to  describe  systems  such  as  diradicals  or  ..vstci.u 
with  orbital  degeneracy.  Because  of  thia,  one  must  truncate  the  expansion  and  eliminate  uninipor 
tant  configurations.  In  general  one  can  show  that 


=  0 


(13) 


In  words,  the  matrix  elements  between  the  reference  wavefunction  and  triple  and  quadruple  exci¬ 
tations  is  zero.  Thus  to  first  order,  only  single  and  double  excitations  contribute.  Because  of  this, 
many  Cl  calculations  attempt  to  include  all  single  and  double  excitations  in  the  expressions  for 
fu  •  To  go  beyond  this,  generally  more  than  one  reference  wavefunction  is  usi,d.  Programs  to  han¬ 
dle  configurations  on  all  the  above  criteria  are  available  at  UTRC. 
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Other,  potentially  more  elegant  methods  of  configuration  choice  involve  formal  approaches 
based  on  natural  orbital  (Reference  24)  or  multiconfiguration  SCF  (Reference  25)  concepts.  To 
implement  the  natural-orbital  approach,  an  initial  limited  Cl  wavefunction  is  transformed  to 
natuial  orbital  form,  and  the  resulting  natural  orbitals  are  used  to  form  a  new  CL  The  desired 
result  is  a  concentration  of  the  bulk  of  the  Cl  wavefunction  into  a  smaller  number  of  significant 
terms.  The  multiconfiguration  SCF  approach  is  more  cumbersome,  but  in  principle  more  . 
effective.  It  yields  the  optimum  c  rbital  choice  for  a  preselected  set  of  configurations.  This 
approach  works  well  when  a  small  i.  jmber  of  dominant  configurations  can  be  readily  identified. 
The  method  is  described  briefly  below. 

2.1 .4  Multiconfiguration  -  Seif  Consistent  Field  Method  (MC-SCF) 

The  Hartree-Fock  self  consistent  field  method  has  been  proven  to  be  a  powerful  tool  for  the 
calculation  and  understanding  of  many  ground  state  properties  of  molecules  in  .he  vicinity  of 
their  equilibrium  structure.  Howe'.'c  •,  in  most  cases  the  one  configuration  Hartree-Fock  approxi¬ 
mation  is  not  adequate  to  pr’peily  describe  the  dissociation  of  molecular  bonds.  Also,  many 
excited  states  cannot  be  repr.  sented  by  a  single  configuration  wave  function.  In  order  to  calculate 
properties  for  such  states,  or  to  investigate  the  formation  of  molecular  bonds,  one  often  needs 
multiconfiguration  wave  functions  for  which  both  the  linear  coefficients,  C,  of  the  configuration 
expansion 


f  =  jQV'i 

I 


(14) 


as  well  as  the  set  of  orthonormal  molecular  orbitals  {^j}  ,  from  which  the  configurations  ip,  are 
constructed,  are  optimized  according  to  the  variational  principle.  As  is  well  known,  this  “MC- 
SCF”  problem  presents  many  more  difficulties  than  the  simple  one  configuration  Hartree-Fock 
case  and  much  work  has  been  devoted  to  obtaining  convergent  solutions  during  the  last  decade. 

The  difficulties  mainly  arise  from  the  fact  that  for  general  MC-SCF  wave  functions  the 
energy  is  not  invariant  with  respect  to  rotations  between  occupied  orbitals.  Hence,  instead  of  a 
relatively  simple  pseudo- eigenvalue  equation  :h».  one  determinant  case,  the  set  of  coupled  Fock 
equations 


(15) 

j  j 
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with  the  hermiticity  conditions 


^ij 


=  eji 


(16) 


has  to  be  solved.  e,j  are  Lagrange  multipliers  which  account  for  the  orthonormality  constraints 
imposed  on  the  orbitals.  The  Fock  operators  F,j  depend  on  the  orbitals  (^,J  and  the  set  of  Cl  co¬ 
efficients  [C,] .  In  analogy  to  the  one  determinant  case  many  attempts  have  been  made  to  solve 

these  equations  iteratively  by  keeping  the  Fock  operators  fixed  in  each  iteration  step.  Then  the 
Lagrange  multipliers  can  be  expressed  by  coupling  operators  constructed  such  that  the  Fock 
equations  are  transformed  into  pseudo-eigenvalue  equations  yielding  the  improved  orbitals.  These 
are  used  in  a  second  step  to  determine  new  Cl  coefficients  by  diagonalizing  the  Cl  matrbc 

.  The  convergence  of  these  algorithms,  however,  has  often  been  found  to  be  poor. 

A  second  group  of  MC-SCF  methods  is  based  on  the  generalized  Brillouin  theorem.  In  these 
methods,  the  orbital  changes  are  derived  from  the  coefficients  of  a  Cl  expansion  consisting  of  the 
MC-SCF  wavefunction  and  alt  one-electron  singly  excited  configurations.  A  computer  program 
(ALIS)  implementing  this  method  has  been  developed  by  Ruedenberg,  et  al.  (Reference  26).  A 
somewhat  more  elegant  program  (GAMESS)  which  also  incorporates  analytical  gradients  is  also 
available  (Reference  23).  More  recently,  various  methods  have  been  proposed  which  are  based  on 
direct  minimization  cf  the  energy,  avoiding  the  Fock  operators  altogether.  Such  methods  are  now 
being  studied  in  our  laboratoiy  and  will  be  incorporated  into  our  existing  computer  programs  if 
they  prove  to  be  highly  efficient. 

2.1 .5  Many-Body  Perturbation  Theory  (MBPT) 

In  MBPT  we  again  begin  with  the  SCF  wavefunction  {%)  as  our  reference  and  attempt  to 
account  for  Eq,  the  correlation  energy.  The  concept  of  excitation  functions  described  in  the  above 
section  on  Cl  calculations  carries  over  to  MBPT  calculations.  In  MBPT  one  can  write  the 
wavefunction,  V'p .  as 


%  =  e>o) 


(17) 
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where  T  is  an  excitation  operator  defined  as 


T  =  Ts  +  Td  +  Tt  +  Tq...  (18) 

where  S,  D,  X  ^nd  Q  refer  to  single,  double,  triple,  and  quadruple  substitutions  respectively.  One 
can  write  Tn,  in  general,  where  n  refers  to  the  number  of  excitations  as 


T 


T„  = 


(19) 


where  a,b,c,...  are  excited  orbitals  and  i,j,k,...  are  orbitals  occupied  in  (po  -  The  total  energy  is  now 
given  by 


Embpt  =  (^o|Ke'^|0o)  (20) 

To  evaluate  Embpt»  the  tf  ■”  from  above  must  be  determined.  An  equivalent  expression  that 
makes  the  perturbation  expansion  clearer  is 


Embpt  =  ^  (^o|K[(Eo-3£o)‘*36j''|0o)  (21) 

k=0  '  ^ 

where  the  sum  is  over  only  so-called  linked  diagrams,  36o  has  eigenfunctions  (f>o .  and  the  expan¬ 
sion  is  of  orders  in  the  perturbation,  V  =  (a6-36o)  •  The  k  =  0  term  gives  the  reference  energy 
and  for  k  >  0  correlation  corrections  are  included.  In  practice  the  MBPT  total  energy  is  calcu¬ 
lated  by  truncating  the  T  operator  expansion  and  projecting  3£e^=[^o}  onto  the  appropriate  n- 

space.  This  leads  to  a  set  of  nonlinear  coupled  equations  for  the  tf  coefficients  which 

correspond  to  the  Cl  expansion  coefficients.  The  equations  are  solved  iteratively  and  Embpt 
evaluated.  In  practice  T4  is  an  upper  limit  that  corresponds  to  quadruple  substitutions.  The  scries 
is  an  oscillatory  convergent  sum,  which  in  practice  has  proven  to  be  at  least  as  accurate  as  Eci 
with  single  and  double  excitations  included. 

The  best  possible  method  to  use  would  be  a  full  Cl  (all  possible  configurations)  with  a 
complete  basis  set.  However,  since  the  number  of  configurations  is  proportional  to  (nf ,  where  n 
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is  the  number  of  basis  functions  and  i  is  the  level  of  excitation,  it  would  be  prohibitive  to  even 
include  all  single  and  double  excitations.  This  truncation  causes  the  loss  of  size  consistency,  which 
implies  that  the  energy  calculated  for  A  and  B  as  a  molecular  system,  but  dimensionally  far  apart, 
is  the  sum  of  the  energy  calculated  for  A  and  B  separately.  In  a  size  extensive  calculation  the 
energy  i^  proportional  to  the  size  of  the  system.  These  properties  are  very  important  if  one  wishes 
to  compute  correct  relative  energies  on  a  potential  energy  surface,  a  necessary  criteria  for  deHning 
the  reactive  pathways  of  interest  in  this  research  program.  MBPT,  on  the  other  hand,  is  guaran¬ 
teed  to  have  the  correct  size -dependence  because  the  expansions  contain  only  the  so-called  linked 
diagrams.  In  addition,  because  of  the  computational  efficiency  of  MBPT,  calculations  can  be 
performed  up  to  fourth  order  in  the  perturbation  expansion  and  include  single,  double,  triple,  and 
quadruple  excitations  in  the  calculation  of  the  correlation  energy.  Such  calculations  are  usually 
performed  with  at  least  a  split  valence  plus  polarization  basis  set.  For  the  light  element  com¬ 
pounds  which  exhibit  ionic  bonding,  the  inclusion  of  diffuse  basis  functions  and  possibly  higher 
polarization  d  or  f-functions  may  be  required  for  a  more  quantitative  treatment. 

2.2  Ab  Initio  Gaussian  Wavefunction  Electronic  Structure  Codes 

Owing  to  the  complexity  of  evaluating  multicenter  electron  repulsion  integrals  over  Slater- 
type  (exponential)  orbitals,  various  groups  have  adopted  a  computational  approach  to  electronic 
structure  calculations  based  on  gaussian  orbitals.  A  highly  developed  series  of  programs,  named 
GAUSSIAN  8. 1  (X  =  2, 6,  8),  is  available  from  Camegie-Mellon  University  (Reference  22).  A 
second  series  of  programs  has  evolved  from  the  original  version  of  Dupuis  and  King’s  (Refer¬ 
ence  27)  HONDO  code.  This  code  has  been  further  developed  by  separate  groups  as  GAMESS 
(Reference  23),  CADPAC  (References  28)  and  HONDO  7  (Reference  29).  Finally  a  new  code 
named  COLUMBUS  (Reference  30  and  31)  has  been  developed  by  the  Ohio  State/'Argonne/Bat- 
telle  group.  UTRC  has  been  designated  as  a  beta  test  site  for  this  new  development.  A  brief 
description  of  the  features  of  these  codes  follows. 

2.2.1  GAUSSIAN  88 

GAUSSIAN  88  is  a  connected  system  of  programs  for  performing  ab  initio  molecular  orbital 
(MO)  calculations.  It  represents  further  development  of  the  GAUSSIAN  101761^1^’%  systems 
already  published.  The  contributors  to  this  program  include;  M.  J.  Frisch,  J.  S.  Bin’Kley,  H.  B. 
Schlegel,  K.  Raghavachari,  C.  F.  Melius,  R.  L.  Martin,  J.  J.  P.  Stewart,  F.  W.  Bobrowic^  C.  M. 
Rohifing,  L.  R.  Kahn,  D.  J.  Defrecs,  R.  Seeger,  R.  A.  Whiteside,  D.  J.  Fox,  E.  M.  Fleuder,  and  J. 
A.  Pople.  GAUSSIAN  8X  was  originally  implemented  on  the  chemistry  department  DEC  VAX 
11/780  computer  at  Carnegie- Mellon  University.  Since  then  this  program  has  been  installed  on  a 
number  of  different  computers. 

GAUSSIAN  88  was  designed  with  a  transparent  input  data  stream,  making  this  program 
very  user  friendly.  All  of  the  standard  input  is  frce-format  and  mnemonic.  Reasonable  defaults 
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for  input  data  have  been  provided,  and  the  output  is  intended  to  be  self-explanatory.  Mechanisms 
are  available  for  the  sophistivated  user  to  override  defaults  or  interface  their  own  code  to  the 
GAUSSIAN  system.  In  this  respect,  we  intend  to  utilize  GAUSSIAN  88  as  a  fundamental 
framework  for  several  applications.  Options  are  being  incorporated  into  this  code  to  provide 
capabilities  beyond  Hartree-Fock  and  various  perturbation  theory  options. 

The  capabilities  of  the  GAUSSIAN  88  system  include: 

a)  Calculation  of  one-  and  two-electron  integrals  over  s,  p,  d,  and  f  contracted  gaussian 
functions.  The  basis  functions  can  either  be  cartesian  gaussians  or  pure  angular  mo¬ 
mentum  functions  and  a  variety  of  basis  sets  are  stored  in  the  program  and  can  be 
requested  by  name. 

b)  Self-consistent  field  calculations  for  restricted  closed-shell  (RHF),  unrestricted  open- 
shell  (UHF),  and  open-shell  restricted  (ROHF)  Hartree-Fock  wavefunctions  as  well  as 
those  types  of  multiconfigurational  wavefunctions  that  fail  within  the  Generalized  Va¬ 
lence  Bond-Perfect  Pairing  (GVB-PP)  formalism. 

c)  Evaluation  of  various  one-electron  properties  of  the  Hartree-Fock  wavefunction,  in¬ 
cluding  Mulliken  population  analysis,  multipole  moments,  and  electrostatic  fields. 

d)  Automated  geometry  optimization  to  either  minima  or  saddle  points,  and  analytical  or 
numerical  differentiation  to  produce  force  constants,  polarizabilities,  and  dipo’*  de¬ 
rivatives.  This  feature  can  be  used  to  develop  minimum  energy  reaction  paths  along  a 
complicated  many  dimensional  potential  energy  surface. 

e)  Correlation  energy  calculations  using  M(J)ller-Plesset  perturbation  theory  carried  to 
second,  third,  or  fourth  order. 

f)  Correlation  energy  calculations  using  configuration  interaction  (Cl)  with  either  all  dou¬ 
ble  excitations  (CID)  or  all  single  and  double  excitations  (CISD). 

g)  Correlation  energy  calculations  using  coupled  cluster  theory  with  double  substitutions 
(CCD). 

h)  Correlation  energy  calculau'ons  using  quadratic  convergence  SCF  (QCSCF).  This  is  a 
new  highly  efficient  size-consistent  method  recently  developed  by  Pople. 

i)  Analytic  computation  of  the  nuclear  coordinate  gradient  of  the  RHF,  UHF,  ROHF, 
GVB-PP,  MP2,  CID  and  RCISD  energies. 

j)  Computation  of  force  constants  (nuclear  coordinate  second  derivaUves), 
polarizabilities,  hyperpolarizabilities,  dipole  derivatives,  and  polarizability  derivatives 
either  analytically  or  numerically. 
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k)  Harmonic  vibrational  analysis. 

l)  Determination  of  intensities  for  vibrational  transitions  at  the  HF,  MP2,  and  Cl  levels. 

m)  Testing  the  SCF  wavefunctions  for  stability  under  release  of  constraints. 

n)  Correlated  electron  densities  and  properties. 

o)  Minimum-energy  pathway  following  from  products  to  reactants  through  a  transition 
state-intrinsic  reaction  coordinate  finder. 

2.2.2  GAMESS 

A  wide  range  of  quantum  chemical  computations  are  possible  using  GAMESS  (Refer¬ 
ence  23),  a  refinement  by  M.  Schmidt  and  S.  Elbert  of  the  original  HONDO  code  (Reference  27). 

The  capabilities  of  this  code  include: 

a)  Calculations  of  RHF/UHF/ROHF/GVB-SCF  molecular  wavefunctions. 

b)  Calculations  of  multiconfiguration  SCF  (MCSCF)  wavefunctions. 

c)  Calculations  of  Cl  wavefunctions  using  the  unitary  group  method. 

d)  Optimization  of  molecular  geometries  using  an  energy  gradient  in  terms  of  Cartesian 
or  internal  coordinates. 

e)  Searches  for  potential  energy  surface  saddle  points. 

f)  Tracing  the  intrinsic  reaction  path  from  a  saddle  point  to  reactants  or  products. 

g)  Computation  of  normal  modes  and  vibrational  frequencies. 

h)  Calculation  of  the  following  properties: 

1.  dipole,  quadrupole,  and  octupole  moments 

2.  electrostatic  potentials 

3.  electric  field  and  electric  field  gradients 

4.  electron  density  and  spin  density 

5.  Muiliken  and  Ldwdin  population  analysis 

6.  localized  orbitals  by  the  Boys  method 

7.  virial  theorem  and  energy  components. 
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GAMESS  is  a  synthesis,  with  many  major  modifications,  of  several  programs.  A  large  part  of  the 
program  is  from  HONDO.  For  pure  sp  basis  sets,  the  HONDO  symmetry  and  supermatrix  proce¬ 
dure  has  been  adapted  from  GAUSSIAN  76  and  GRADSCF  integrals,  both  for  the  SCF  and  gra¬ 
dient  parts.  The  GVB  section  is  a  heavily  modified  version  of  GVBONE.  A  Boys  [.realization 
algorithm  is  implemented  from  a  heavily  modified  version  of  Streitweiser’s  QCPE  program. 

The  Cl  module  is  based  on  Brooks  and  Schaefer’s  unitary  group  program  which  was 
modified  to  run  within  GAMESS,  using  a  Davidson  eigenvector  method  written  by  S.  T.  Elbert. 
The  MCSCF  module  is  a  NevH'ton-Raphson  procedure,  developed  at  NRCC,  based  on  the  unitary 
group  Cl  package.  The  intrinsic  reaction  coordinate  pathfinder  was  written  at  North  Dakota  State 
University. 

2.2.3  CADPAC 

The  Cambridge  Analytic  Derivatives  Package  (CADPAC)  (Reference  28)  is  a  group  of 
programs  which  has  been  under  development  at  Cambridge  University,  UK.  It  originated  as  a 
version  of  Dupuis  and  King’s  HONDO  program.  From  its  initial  state  as  an  SCF  gradient 
package  for  closed-shell  and  UHF  wavefunctions,  this  program  has  been  extensively  modified 
with  many  of  the  old  features  being  enhanced  and  many  new  features  being  added.  The  input  data 
is  now  in  free  format  with  a  ‘keyword’  system,  to  make  the  program  easier  to  use.  The  integral 
routines  use  essentially  the  same  methods  as  those  in  HONDO,  being  based  on  the  Rys  polyno¬ 
mial  method,  but  have  been  extended  to  cover  f-functions.  These  routines  have  also  been 
vectorized  to  take  advantage  of  the  much  greater  availability  of  supercomputers  such  as  the  Cray 
in  recent  years.  The  SCF  programs  are  by  now  a  blend  of  techniques,  but  still  contain  a  few 
features  from  the  initial  HONDO  program,  particularly  in  the  way  symmetry  is  handled.  There 
have  been  many  modifications  to  improve  efficiency  and  the  .''.ddition  of  level-shifting  and 
damping  techniques  and  the  implementation  of  the  DIIS  method  lo  aid  convergence.  A  high-spin 
open-shell  SCF  program,  and  a  completely  general  open-shell  SCF  program  have  also  been 
Included. 

The  gradient  routines  have  been  considerably  altered  from  those  in  the  original  HONDO, 
the  method  having  been  changed  to  one  which  is  more  efficient  and  easier  to  extend  to  higher 
order  derivatives.  These  routines  also  work  for  f-functions.  The  original  program’s  capabilities  for 
the  calculation  of  force-constants  by  numerical  differentiation  of  gradients  and  for  the  optimiza¬ 
tion  of  geometries  remain  essentially  intact.  However  there  is  now  a  choice  of  two  optimization 
algorithms  with  the  inclusion  of  Schlegel's  method. 

In  addition  to  the  above  modifications  a  range  of  extra  facilities  were  developed.  These 
included  a  mo’"‘  nowerful  method  of  calculating  one-electron  properties  and  analyzing  molecular 
charge  distributions,  and  various  ‘post-SCF’  stages,  beginning  with  a  4-index  transformation. 
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These  new  sections  include  M(j)ller-PIesset  perturbation  theory  for  total  energies  and  molecular 
properties,  and  coupled  Hartree-Fock  calculations  of  molecular  polarizabilities.  The  polarizability 
routines  can  calculate  dynamic  properties  at  real  or  imaginary  frequencies,  and  obtain  dispersion 
coefficients.  It  is  also  possible  to  use  CHF  theory  to  obtain  the  perturbations  due  to  nuclear 
displacements.  These  can  be  used  to  obtain,  analytically,  all  the  dipole  and  quadrupole  moment 
derivatives  of  a  molecule.  They  also  form  part  of  the  most  important  addition  to  the  package,  the 
section  which  calculates  analytic  second  derivatives  of  the  energy.  This  is  a  powerful  technique 
whose  speed  and  accuracy  represents  a  considerable  improvement  over  numerical  differentiation. 

The  capabilities  of  the  latest  CADPAC  code  (Version  4)  include 

a)  the  evaluation  of  one-  and  two-electron  integrals  over  contracted  cartesian  gaussian 
basis  functions  of  type  s,  p,  d,  or  f. 

b)  SCF  calculations  for  closed  shell,  open-shell,  UHF  and  generalized  open-shell  tech¬ 
niques. 

c)  calculation  of  one-electron  properties  for  these  types  of  wavefunction,  including  a  dis¬ 
tributed  multipole  analysis. 

d)  calculation  of  the  gradients  of  the  SCF  energy. 

e)  use  of  the  gradients  for  automatic  geometry  optimization,  and  for  the  calculation  of 
force  constants  by  numerical  differentiation.  There  is  a  choice  of  two  optimization  al¬ 
gorithms. 

f)  transformation  of  the  integrals  from  the  atomic  orbital  to  the  molecular  orbital  basis. 

g)  M<})ller-Plesset  perturbation  theory  calculations  to  third  order  in  the  energy  and  sec¬ 
ond  order  in  the  one-electron  properties. 

h)  coupled  Hartree-Fock  calculations  of  polarizabilities,  including  frequency  dependence, 
and  magnetizabilities. 

i)  coupled  Hartree-Fock  calculations  of  the  perturbation  due  to  nuclear  displacements. 

j)  calculation  of  the  dipole  and  quadrupole  moment  derivatives. 

k)  calculation  of  the  second  derivatives  (force  constants)  of  the  energy  by  analytic  meth¬ 
ods. 

l)  analytic  calculations  of  polarizability  derivatives. 

m)  calculation  of  infrared  and  Raman  intensities,  and  the  study  of  vibrational  circular 
dichroism. 
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n)  calculation  of  MP2  gradients,  dipole  moment  derivatives,  polarizabilitites  and  force 
constants  using  analytic  algorithms. 

o)  spin-projected  UHF  MP2  energies 

2.2.4  HONDO 

The  HONDO  program,  (Reference  27)  originally  developed  by  Dupuis  and  King  at  NRCC 
has  recently  been  refined  and  updated  by  M.  Dupuis  at  IBM-Kingston  (Reference  29). 

The  following  features  are  available  in  the  present  version  (HONDO  7)  of  the  program: 

a)  Single  configuration  self-consistent-field  wavefunctions  (closed  shell  RHF,  spin  unre¬ 
stricted  UHF,  restricted  open  shell  ROHF),  generalized  valence  bond  GVB  and  gen¬ 
eral  multiconfiguration  self-consistent-field  MCSCF  wavefunctions,  and  configuration 
interaction  Cl  wavefunctions  can  be  calculated. 

b)  The  electron  correlation  correction  to  the  energy  of  closed  shell  RHF  wavefunctions 
can  be  calculated  by  means  of  M<|)ller-Plessel  (MP)  perturbation  theory  applied  to 
second-,  third-,  and  fourth-order  (with  or  without  the  effects  of  triple  excitations). 

c)  The  effective  core  potential  approximation  can  be  used. 

d)  Optimization  of  molecular  geometries  using  the  gradient  of  the  energy  with  respect  to 
nuclear  coordinates  is  possible  with  all  but  the  Cl  and  MP  wavefunctions.  Optimiza¬ 
tion  can  be  carried  out  in  the  cartesian  space  or  in  the  internal  coordinate  space  with 
the  possibility  of  freezing  some  cartesian  or  internal  coordinates. 

e)  The  force  constant  matrix  in  the  cartesian  space,  and  the  vibrational  spectrum,  includ¬ 
ing  infrared  and  Raman  intensities,  can  be  calculated  with  all  but  the  Cl  and  MP 
wavefunctions. 

f)  Calculation  of  the  dipole  moment  and  polarizability  derivatives  with  respect  to  the  nu¬ 
clear  coordinates  is  possible,  for  use  with  a  previously  calculated  force  constant  ma¬ 
trix.  The  force  constant  matrbc  can  be  transformed  to  the  internal  coordinate  basis. 

g)  Transition  state  structures  can  be  determined  with  all  but  the  Cl  and  MP  wavefunc¬ 
tions  by  taking  advantage  of  the  energy  gradients. 

h)  The  ‘Intrinsic  Reaction  Coordinate’  (IRC)  pathway  can  be  determined  with  all  but  the 
Cl  and  MP  wavefunctions  by  taking  advantage  of  the  energy  gradients. 

i)  Molecular  energies  for  several  points  on  a  potential  energy  surface  can  be  calculated  in 
a  single  run. 


19 


j)  Non-gradient  optimization  of  basis  function  exponents  is  possible.  The  source  code 
can  be  modified  to  carry  out  optimization  of  other  non-linear  parameters,  for  example, 
contraction  coefficients  and  geometrical  parameters. 

k)  The  following  electronic  properties  can  be  extracted  from  the  wavefunction: 

1.  Dipole  moment 

2.  Quadrupole  moment 

3.  Mulliken  population,  bond  order  and  valency  analyses 

4.  Spin  density  maps 

5.  Electron  density  maps 

6.  Electrostatic  potential  maps 

7.  Localized  orbitals  via  Boys’  method 

8.  Static  dipole  polarizability 

9.  Static  first  and  second  hyperpolarizabilities. 

l)  The  potential  due  to  finite  point  charges  for  a  classical  representation  of  an  environ¬ 
ment,  or  a  uniform  electric  field  can  be  incorporated  into  the  one-electron  Hamil¬ 
tonian. 

2.2.5  COLUMBUS 

The  COLUMBUS  code  (Reference  30  and  31)  is  a  continuing  development  based  on  a  joint 
project  at  Ohio  State  University  (I.  Shavitt),  University  of  Karlsruhe  (R.  Ahlrichs),  and  Argonne 
National  Laboratory  (R.  Shepard).  The  unique  features  of  this  code  are  the  incorporation  of 
relativistic  or  non-re!ativistic  core  potentials  to  permit  analysis  of  heavy  atom  molecular  systems 
and  the  inclusion  of  a  sophisticated  Cl  package  based  on  the  unitary  group  approach  of  Shavitt 
(Reference  32).  The  basic  programs  included  in  the  COLUMBUS  code  are  as  follows. 

1.  AO  Integrals 

This  is  R.  Pitzer’s  version  of  the  integral  package  from  HONDO  with  partial  vectorization 
of  auxiliary  integral  routines.  This  package  can  handle  up  to  g-functions. 

2.  Integral  Hansformation 

The  transformation  algorithm  is  written  over  all  the  orbitals  which  constitute  a  shell,  using 
contraction  coefficient  Jefined  over  the  primitive  basis  functions.  Limited  vectorization  is 
possible  over  the  innermost  loops  for  dyadic  operations. 
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3.  SCF  Gradients 

These  routines  are  based  on  the  HONDO  versitii  but  include  checks  on  integral  symmetry 
to  avoid  operations  over  zero  or  near  zero  elements.  At  this  time,  only  a  first  derivative 
analysis  has  been  included. 

4.  SCF  Energy 

The  SCF  routines  are  typical  of  those  found  in  the  molecular  structure  codes.  They  incor¬ 
porate  level  shifting,  damping  and  the  incorporation  of  Pulay’s  DIIS  convergence  accelera¬ 
tion  procedure. 

5.  MCSCF  Analysis 

The  MCSCF  package  is  due  to  R.  Shepard  and  incorporates  extensive  vectorization.  The 
output  vectors  can  be  taken  directly  from  this  package  and  transferred  to  the  Cl  program. 

6.  Multireference  Direct  Cl 

This  set  of  routines  is  based  on  the  graphical  representation  of  the  unitary  group  ap¬ 
proach  (GUGA)  for  constructing  a  Cl  wavefunction.  This  code  differs  from  conventional 
Cl  in  that  a  matrix  representation  of  the  hamiltonian  is  not  explicitly  computed.  The  uni¬ 
tary  group  approach  can  be  much  more  efficient  in  most  applications  since  ~80%  of  the 
hamiltonian  matrbc  elements  are  usually  zero  but  must  still  be  included.  This  leads  to  the 
storage  and  diagonalization  problems  associated  with  large,  sparse  matrices.  However,  lit¬ 
tle  work  has  been  done  on  analysis  of  excited  states  with  the  same  molecular  symmetry  as 
the  ground  state.  An  approach  that  looks  promising  is  to  define  an  approximate  vector 
corresponding  to  the  excited  state  wavefunction  and  to  iteratively  improve  upon  this  solu¬ 
tion  using  the  direct  Cl  contraction.  The  uncertainties  in  the  procedure  are  applications  to 
situations  where  degeneracies  or  near-degeneracies  arise.  Further  studies  of  this  case  are 
in  progress. 

2.3  Spin-Projected  Unrestricted  Hartree-Fock  Method 

The  unrestricted  Hartree-Fock  (UHF)  method  developed  by  Pople  and  Nesbet  (Refer¬ 
ences  33  and  34)  yields  the  best  single-determinant  approximation  to  the  exact  wavefunction  for 
an  atomic  or  molecular  system.  Such  a  wavefunction  incorporates  correlation  by  allowing  orbitals 
of  different  spin  to  adjust  to  spatially  different  forms,  thus  breaking  the  symmetry  restrictions  of 
the  conventional  (RHF)  method  (References  35  and  36).  It  is  necessary,  however,  to  project  from 
such  a  wavefunction,  a  properly  antisymmetrized  spin  and  angular  momentum  state  in  order  to 
define  eigenstates  and  eigenenergies  corresponding  to  observable  spectroscopic  states. 

Let  \a),  [(3)  refer  to  the  doubly-occupied  molecular  orbitals  (MO’s),  (pa,  (pp  •  Let  |yM<5},|r) 
refer  to  the  singly  occupied  MO’s.  We  assume  that  all  MO’s  have  been  subjected  to  a 
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transformation  to  orthogonal  form.  Let  Cj^  be  the  coefficient  associated  with  the  permutation 
y  * — *  d  (this  permutation  is  (-1)  times  the  overlap  of  the  permuted  spin  eigenfunction  with  the 
original  spin  eigenfunction),  and  adopt  the  convention  Cyy  =  -1.  Then  the  expectation  value  of 
the  nonrelativistic  Hamiltonian 


M  =  2U(i)  +  Z  V(i.j) 

i  i<j 


is  given  by 


(36)  =  22HU|a)  +  +  2  IMvM 

a  Y  afi 

afi  ay  ay 


(22) 


(23) 


where  i,  j  refer  to  electron  numbers,  a,  y,  d  etc.,  refer  to  MO’s  and  the  sums  run  over  all 
a,  /S.... 

The  equations  for  determining  the  optimum  MO’s  can  be  derived  from  the  following 
variational  form 


(k)  -  2  5;  {a  1)8)  -  Y^{€ya{a  I  y)  +  ^ay  (v  I  a})  -  Z  (y  1 

afi  ap  yd  . 


=  0 


(24) 
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These  equations  are 


(U  +  2Jc-K,  +  Jo-iKo)|«)  = 


(25) 


where 


(U  +  ZI<,-K<,  +  +  2C^l(«|V|y)a)  =  2e/i,l^)  +  lejvli  )  (26) 

6  fid 


l((5lV|y)(5)=|  dP0i(r)V(r;r)^y(r)^<5(r) 


(W|V|W')  =  I  dFW'-(i)V(r,?)W'(P) 


Jc  =  2mm 

p 


Jo  =  2{^\m 


(27) 


Kc|W  =  Xl(i5|V|W)/S)  ) 

P  6 


Left-multiplying  Equation  (Z'?)  by  (  y|  and  Equation  (26)  by  ( a|  ,  these  equations  can  be  com¬ 
bined  to  yield 


cya 


=  XI^Cie+lKv'SIVISa) 


(28) 


Noting  that 


ZCj^<5r|Vliy)  =  2;c*((5r|V|zy). 


(29) 
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Equations  (25)  and  (26)  can  be  recast  as  follows.  From  Equation  (25), 

6 

with  (30) 


M  =  X(2C4,+  l)l'5)(f('5|V|r)| 

dr 


From  Equation  (26), 


+  i)\D){fir\'/\6y)  =  OcM’iv) 

p  px 


(?1) 


with 


Oc=2l/’)(^li  Oo=El«)(<S|  (32) 

p  d 


Using  Equations  (31)  and  (32)  in  (25)  and  (26),  we  find 


(U  +  2Jc-Kc  +  Jo-^Ko-|M)|a)  =  Y.^paW) 

Z  Z.  p 


(33) 


(U  +  2Jc-Kc  +  Jo-|Ko  +  |Mt-OcM0ly>=  > 


(34) 
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Finally,  Equations  (33)  and  (34)  can  be  rearranged  to  yield  the  same  Hermitian  operator  on  the 
left-hand  side  with  the  resulting  definition  of  the  one-electron  eigenenergies  for  the  closed  and 
open-shell  eigenstates: 


ea'  =(«|U  +  2Jc-Kc  + Jo-^Kola  } 


(35) 


6/  =(y|U  +  2Je-Kc  +  Jo-|Ko  +  ^M^ly) 


(36) 


Equations  (35)  and  (36)  represent  the  correct  one-electron  energy  expressions  for  the  spin- 
projected  eigenstate. 

This  formalism  has  been  developed  to  extend  the  conventional  HF  method  to  include 
split-shell  correlation  and  proper  spin  and  symmetry  projections.  It  is  being  incorporated  into  a 
computer  program  using  gaussian-type  orbitals  (GTO’s)  as  the  elementary  basis  functions.  A 
crucial  feature  of  this  method  is  that  dissociation  always  follow.^  the  lowest  energy  pathwaj 
thereby  permitting  a  proper  description  of  bond  foimation  and  breakage.  In  contrast,  RHF  or 
MC-SCF  methods  often  exhibit  improper  dissociation  character  or  exhibit  size  inconsistencies 
owing  to  correlation  energy  changes  in  going  from  molecular  geometrie.s  to  separated  atom- 
molecule  or  atom-atom  dimensions. 

2.4  Dissociative  Recombination  Calculations 

Tne  .neory  of  the  capture  of  an  electron  by  a  positive  molecular  ion  has  been  developed  for 
both  direct  attachment  processes  (Reference  37),  and  for  several  possible  indirect  processes  which 
involve  the  formation  of  an  intermediate  excited  Rydberg  state  belonging  to  the  molecular  ion 
core  configuration  (References  38-40).  Let  V'f  (£,»  represent  the  continuum  wavefunction 
describing  the  free  electron  with  energy  e  plus  the  n-1  electron  bound  system  of  the  molecular 
ion.  Let  (t,>  represent  the  wavefunction  for  an  eigenstate  of  the  neutral  molecule  which 
can  be  written,  in  the  Born-Oppenheimer  approximation,  as  the  product  of  an  electronic  and 
nuclear  wavefunction  in  the  form 


V;d(r.R)  =  t/»§'(r.R)^d(R) 


(37) 


The  cross  section  for  electron  capture  in  dissociative-recombination  is  determined  by  the 
asymptotic  form  of  the  nuclear  wavefunction  ^d(R)-  Let  the  incident  flux  of  electrons  per  unit 
area  be 


2n^h 


de 


(38) 


where  (56  is  a  measure  of  the  energy  spread  of  the  incoming  beam  and  k  is  the  wave  number 
associated  with  the  incoming  electron.  The  outgoing  flux  of  atoms  is  given  by 


M 


l&(R)l^ 


(39) 


where  —  is  the  relative  velocity  of  the  separating  neutral  atoms.  Equations  (37)  and  (39)  yield 
for  the  cross  section 


o(e,Q)=  lim^^^I^d(R)l2 

R-*  00  flk 


(40) 


Equation  (40)  must  finally  be  averaged  over  all  rotational  orientations  as 


<7(e) 


Q)dQ 


(41) 


Assuming  that  the  total  cross  section  can  be  written  as  the  product  of  a  resonant  capture  cross 
section  and  a  survival  factor,  in  the  form 


o(€)  =  <7cap(e)S(e) 


(42) 
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we  have 


(r,R,e)lV(R,e)|v;d(r,R))2  (43) 

where  the  factor  tWo/w+  is  the  ratio  of  the  electronic  degeneracies  of  the  neutral  and  ionic  states, 
respectively.  The  capture  width,  Tc  (R,  e) ,  is  determined  by  integrating  over  the  coordinates  of 
the  electronic  wavefunctions 


rc(R.e)  =  2ji  {rpi(s,R)ii€XR)\%,i-€\ri^Ul'^)Y 


(44) 


Thus,  the  capture  cross  section  can  be  written  as 


/  Wo  \  1 


(45) 


or  equivalently 


(5+(R)V'i(R)«'i(r.R,e)lV(R,e)|v-d(f.R)?„(R)): 


2m€  CO+ 


(46) 


The  nuclear  wavefunction,  ^o>  is  energy  normalized  so  that  asymptotically 


(47) 
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Equation  (46)  can  be  cast  into  the  computational  form 


(7(e,V)cmVsec 


1.38188  X 
e(a.u.) 


{ri(R)i 


ry\R.e) 


l&(R) )' 


(48) 


where  the  electron  energy  e  and  the  capture  width  Fc  are  in  atomic  units. 

Assuming  a  Maxwellian  temperature  distribution  for  the  electrons,  the  rate  coefficient  can 
be  written  as 


«v(Te)  =  j  *  ''el  *  (49) 

or  equivalently 

9/9  f  “ 

The  capture  width  Fc  can  be  calculated  by  examining  the  high  members  of  the  Rydberg  series  of 
neutral  states  which  have  the  structure  of  a  ground  state  molecular  ion  coupled  to  an  electron  in 
diffuse  hydrogenic  orbitals  of  large  effective  principal  quantum  numbers.  A  direct  calculation  of 
the  capture  width  involves  the  knowledge  of  the  continuum  wavefunction  for  the  electron  as  a 
function  of  the  interparticle  coordinates.  This  approach  is  computationally  very  difficult  and  at 
the  present  stage  of  development  would  probably  lead  to  errors  of  at  least  the  same  magnitude  as 
an  extrapolation  procedure  of  the  corresponding  neutral  Rydberg  states. 

2.4.1  Indirect  Models 

An  indirect  model  for  calculating  electron-molecular  ion  recombination  coefficients  can  also 
be  formulated.  Here  we  postulate  the  formation  of  a  collision  complex  (AB*  for  purposes  of 
illustrating  this  for  a  diatomic  system)  during  the  collision  and  its  subsequent  decay  to  form  A*, 
B**,  or  AB^  +  e  (References  38-40).  The  approaches  differ  by  the  nature  of  the  (AB*)  complex. 
The  high-energy  approximation  assumes  that  all  of  the  energy  of  the  incident  electron  goes  into 
exciting  one  of  the  electrons  in  the  AB^  core,  resulting  in  the  capture  of  the  incident  electron  into 
a  doubly  excited  state  of  AB.  This  mechanism  is  unlikely  to  be  applicable  to  AB  recombination 
at  thermal  energies  because  of  the  energy  required  to  promote  an  electron  out  of  the  AB  *  core.  A 
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second  approach  is  more  applicable  to  thermal  energy  collisions,  since  it  requires  that  the 
electron-molecular  ion  collision  results  in  a  vibrationally  excited  mole'-ulai  ion  which  captures  an 
electron  in  the  Coulomb  field  of  AB'*'.  This  collision  complex  is  equivalent  to  a  vibrationally  or 
rotationally  excited  Rydberg  state  of  AB.  Experimental  evidence  for  the  existence  of  such  states 
and  their  autoionization  has  been  independently  obtained  from  photoionization  spectra  of  AB, 
within  a  few  electron  volts  of  the  threshold  (Reference  41).  Wth  this  model,  the  prediction  of  the 
atomic  products  of  the  recombination  reaction  reduces  to  calculating  the  Rydberg  states  of  AB 
that  are  most  likely  to  be  populated  during  the  initial  recombination,  and  also  those  atomic  states 
to  which  they  subsequently  decay. 

As  contrasted  with  the  direct  recombination  process  which  is  governed  by  the  configuration 
interaction  strength,  the  nonadiabatic  coupling  of  the  electron  and  nuclear  motion  ‘  the  dominant 
mechanism  in  both  stages  of  the  indirect  recombination  process.  In  the  first  stage,  it  is  the 
vibrational  (or  rotational)  excitation  of  the  AB'*'  core  which  results  in  the  initial  capture  of  the 
incident  electron.  The  subsequent  decay  of  the  collision  complex  is  due  to  a  transition  to  a 
repulsive  molecular  state  which  crosses  the  Rydberg  state.  It  is  then  the  nonadiabatic  coupling  of 
the  electronic  and  nuclear  motion  which  results  in  the  transition  between  these  two  molecular 
states.  Actually,  it  is  the  same  term  in  the  Hamiltonian  for  the  entire  molecular  system  that  causes 
the  key  transition  in  both  stages  of  the  recombination  process.  This  becomes  obvious  in  the 
quantum  mechanical  formulation  of  the  problem.  A  possible  competing  process  for  the  decay  of 
the  AB*  complex  is  the  radiative  transition  to  lower-lying  states  of  the  molecule.  The  importance 
of  this  effect  can  be  estimated  from  calculations  of  the  band  transition  probabilities. 

Implicit  in  this  model  of  dissociative-recombination  processes  is  the  existence  of  three 
excited  electronic  states  of  AB.  Since  this  recombination  process  is  envisioned  as  occurring  in  two 
stages,  we  are  concerned  with  only  two  of  these  states  at  a  time.  The  state  common  to  both  stages 
of  the  recombination  process  is  the  highly  excited  Rydberg  state  (AB*),  which  is  bound  with 
respect  to  dissociation  to  AB'^  and  a  free  electron.  For  the  calculation  of  the  initial  rate  of 
recombination,  a  wavefunction  for  a  continuum  Rydberg  state  is  needed.  This  wavefunction  has 
the  asymptotic  behavior  of  a  free  electron  moving  in  the  field  of  AB”*". 

To  determine  the  transition  rates  in  the  two  stages  of  the  recombination  process,  it  is 
necessary  to  identify  the  terms  in  the  Hamiltonian  that  are  responsible  for  the  transition.  The 
Hamiltonian  for  the  complete  molecular  system  with  respect  to  the  center  of  mass  is 


36  =  36o- 


y 
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where  R  is  the  vector  joining  nucleus  B,  is  the  Elamiltonian  for  the  system  when  the  nuclei  are 
held  fixed,  ft  is  the  reduced  mass  of  the  two  nuclei,  M  is  the  sum  of  the  nuclear  masses,  and  the 
summations  extend  over  all  of  the  electrons.  Since  all  of  the  electronic  eigenfunctions  under 
discussion  are  implicitly  dependent  on  the  internuclear  separation  |  R  |  ,  it  is  the  term 

3£  -  -  which  plays  a  dominant  role  in  coupling  the  electronic  and  nuclear  motion 

giving  rise  to  the  appropriate  transition  (References  42-45).  This  becomes  more  obvious  when  we 
write  the  total  wavefunction  for  the  collision  complex,  AB*,  as 


V'c*=;fc(re,|Rl)Fc(R)  (52) 

where  Xc  (Fe,  |R|)  is  the  electronic  wavefunction  for  the  appropriate  Rydberg  state  of  AB  at 

internuclear  separation  |R|  and  Fc(R)  is  the  wavefunction  associated  with  the  vibration  and 

rotation  of  the  nuclei.  In  a  similar  manner,  the  total  wavefunction  associated  with  the  repulsive 
state  is 


Vr  =  xxn,  |R|)F#) 


(53) 


where  Ff(R)  is  the  wavefunction  for  the  neutral  atomic  products  as  they  separate.  With  the  use 

of  Fermi’s  golden  rule,  the  transition  rate  for  the  decay  of  the  collision  can  be  obtained  by 
evaluating  the  matrix  elements; 


=yl{V'ciK'IV'r)l^ 


(54) 


The  initial  rate  of  recombination  into  the  highly  excited  Rydberg  states  can  be  calculated  in  a 
similar  manner. 

2.5  Dissociative  Attachment  and  Vibrational  Excitation  Calculations 

A  model  that  has  been  quite  successful  in  explaining  the  dissociative  attachment  (DA)  and 
vibrational  excitation  (VE)  of  diatomic  molecules  is  the  resonance  model  in  which  the  projectile 
electron  is  temporarily  trapped  by  the  target  molecule.  The  molecular  anion  (or  the  resonant 
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state)  thus  formed  has  a  finite  lifetime  and  it  can  either  autodetach  leading  to  VE  of  the  molecule 
or,  if  the  lifetime  is  sufficiently  large,  it  can  lead  to  DA  forming  a  neutral  atom  and  an  atomic 
anion. 


A  few  points  of  this  model  are  the  following;  first,  the  dissociative  electron  attachment  and 
the  vibrational  excitation  of  the  molecule  are  two  possible  decay  channels,  apart  from  electronic 
excitation  etc.,  resulting  from  a  particular  resonance  state.  Thus  a  calculation  of  the  cross  section 
for  the  dissociative  electron  attachment  to  a  molecule  will  provide  resonant  contributions,  of  that 
particular  resonant  state,  to  the  cross  sections  for  vibrational  excitation  of  the  molecule  and  vice 
versa.  Second,  in  explaining  vibrational  excitation  by  this  resonance  model,  it  is  implicitly 
assumed  that  the  transition  from  the  resonant  state  V~(R)  to  the  neutral  state  Vo(R)  after 
autodetachment  is  a  Franck-Condon  transition,  that  is,  an  instantaneous  transition  with  no 
change  in  nuclear  velocities  or  positions.  This  is  a  so-called  local  complex  potential  model.  This  is 
true  only  if  the  energy  of  the  projectile  electron  is  much  greater  than  the  vibrational  spacing.  At 
very  low  impact  energies  or  if  the  vibrational  spacing  of  the  molecule  is  relatively  large,  a 
description  of  the  DA  and  VE  proce.'^ses  using  a  non-local  complex  potential  for  the  resonant 
state  is  essential  since  the  neutral  molecule,  for  its  vibrational  excitation,  can  accept  only  quanta 
of  vibrational  energy.  Third,  in  the  cases  of  some  molecules  it  might  be  possible,  for  certain 
electron  energies,  to  form  more  than  one  intermediate  resonant  state.  Alternatively,  the  resonant 
anion  state  may  decay  into  more  than  one  electronic  state  of  the  neutral  molecule.  In  such  cases, 
the  total  width  T  (R)  is  the  sum  of  various  partial  widths  -  each  partial  width  corresponding  to  a 
certain  transition  between  the  resonant  anion  state  (or  states)  and  the  neutral  molecular  state  (or 
states). 


2.5.1  Dissociative  Attachment 

The  theory  of  the  capture  of  an  electron  through  the  p.ocess  of  dissociative  attachment  (DA) 
to  a  molecule  parallels  that  developed  above  for  the  dissociative  recombination  process.  The 
process  of  dissociative  attachment  to  a  neutral  molecule,  XY,  can  be  written  as 


e  +  XY 


X  +Y- 
X-  +  Y 


(55) 


where  X,  Y  are  the  dissociation  fragments.  The  cross  section  for  dissociative  attachment  is 
obtained  by  comparing  the  flux  of  the  incoming  electrons  per  unit  area  to  the  outgoing 


31 


ion-neutral  pair  flux  defined  by  the  asymptotic  continuum  wavefunction.  We  find,  for  energy 
normalized  asymptotic  functions, 


<7da 


/rk? 


lim  |^j(R)| 

R-*oo 


(56) 


where  k,  is  the  wave  number  associated  with  the  incoming  electron  and  ^  is  the  relative  velocity 
of  the  asymptotic  ion-neutral  pair.  The  partial  wave,  satisfies  the  radial  equation: 


'■  ^  dr^ 


V-(R)-E)fj(R) 


=  -  X  ^vje(R)Ivj(R) 

V 


6vVi“i^ 


I 


dR'^vj(R')V(R')^vj(R') 


(57) 


Here: 


Vvje(R)  =  (  VVje(R)  1 3fiei  -  E  I  0d(R)  )  (58) 


V'vje(R)  is  the  wavefunction  of  the  neutral  molecule  in  vibrational  state  V  in  the  presence  of  the 
scattering  electron.  0d(R)  is  the  wavefunction  of  the  anion  at  short  internuclear  separations  as 

described,  for  example,  by  a  stabilization  calculation.  If  the  electron  energy  is  large  or  the 
vibrational  spacing  is  small,  one  can  replace  E  -  Ey  by  the  local  classical  electron  energy  (in  a.u.): 


V-(R)-Vo(R) 


kI(R) 

2 


(59) 


32 


Then  using  the  closure  property: 


EfvXRKvjCR')  =  a(R-R')  («,) 

V,J 

Equation  (57)  reduces  to  a  local  width  model: 

=  Cv>j(R)[ro(R)/2;rko(R)]'/' 

where  ^v,  j  (R)  is  the  wavefunction  of  ti.u  initial  rovibrational  state  of  the  target  molecule.  The 
resonant  wavefunction  ^j(R)  can  be  obtained  by  numerical  integration  of  Equation  (61)  using,  for 
example,  the  Numerov  technique. 

2,5.2  Vibrational  Excitation 

Vibrational  excitation  of  a  molecule  can  be  written  schematically  as 


e(ki)  +  XY(Vi)-e(kf)  +  XY(Vf) 


(62) 


When  the  electron  is  far  from  the  molecular  target,  we  have  V[e  +  XY(V,)]  =  0.  The  initial 
wavefunction  can  then  be  written  as 


=  A(ki)e-“'>’'^V'^  Xv, 


(63) 


where  A(k,)  is  the  initial  plane  wave  amplitude,  is  just  the  electronic  wavefunction  of  the 

unperturbed  target.  In  a  similar  fashion,  we  can  write  for  the  asymptotic  wavefunction  after 
vibrational  excitation 
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(64) 


V'f  =  A(kf)e-*'''>^ 

Conservation  of  energy  requires 

E(Vi)  +  ^  =  E{Vf)  +  ^  (65) 

Finally  the  cross  section  for  vibrational  excitation  can  be  written  as 

<7v,-*v,  =  ^  I  dk|  1  Tif  1 2  (66) 

where  the  transition  matrix  element  from  the  initial  to  the  final  excitation  state  is  obtained  from 

Equation  (67)  can  be  further  decomposed  into  resonant  and  nonresonant  contributions  and  a 
partial  wave  expansion  corresponding  to  the  outgoing  plane  wave. 

A  computer  program  for  calculating  dissociative  attachment  cross  sections  has  been 
developed  based  on  the  local  width  model.  This  model  is  generally  valid  when  the  electron 
collision  energy  is  large  relative  to  the  vibrational  and  rotational  spacings  of  the  neutral  target. 

For  the  problem  of  interest  in  this  research,  the  dissociative  attachment  of  electrons  to  silane,  the 
local  width  approximation  should  be  valid  for  collision  energies  greater  than  ~  1.0  eV.  This  is 
approximatelj’  five  times  the  vibrational  quantum  of  the  fundamental  mode  in  silane.  An  upper 
limit  is  determined  by  the  convergence  properties  of  the  outgoing  plane  wave  representation  of  the 
dissociation  species,  X*  +  Y  or  X  +  Y*.  In  practice,  this  upper  limit  is  10  to  20  eV,  depending  on 
the  extent  of  the  long-range  potential.  For  dissociative  attachment  at  lower  collision  energies,  the 
non-local  width  model  described  by  Bardsley  (Reference  46)  must  be  employed.  In  either  method, 
the  calculation  of  the  complex  part  of  the  interaction  potential  (capture  width)  is  the  most  difficult 
aspect  of  the  theoretical  analysis. 
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2.6  Discussion  of  Charge  Transfer  Calculations 

Low-to-intermediate  energy  ( <  lOOeV)  ion-atom  and  ion-molecule  charge  transfer 
reactions  play  an  important  role  in  many  physical  processes.  Until  recently,  comparatively  little 
effort  has  been  devoted  to  theoretical  studies  of  the  appropriate  cross  sections.  In  the  past,  both 
theoreticians  and  experimental  physicists  have  found  it  easier  to  study  high-energy  (>  1  keV) 
collisions.  At  these  energies,  the  two  colliding  particles  preserve  their  identities,  ard  it  is  possible 
to  treat  the  interaction  between  them  as  a  perturbation.  There  is  no  guarantee  that  this  procedure, 
known  as  the  Born  approximation,  will  always  converge  to  the  correct  result  (Reference  47).  As 
the  energy  of  the  colliding  particles  decreases,  it  is  necessary  to  take  account  of  the  distortion 
these  particles  undergo  during  the  collision.  The  method  of  perturbed  stationary  states  was 
developed  for  calculating  charge  transfer  and  electronic  excitation  cross  sections  in  relatively  slow 
collisions  between  heavy  particles  (Reference  48).  The  method  has  been  presented  in  both  wave 
and  impact  parameter  formalisms.  In  the  first  of  these,  the  entire  system  is  treated  quantum-me- 
chanically,  while  in  the  latter  the  nuclei  are  assumed  to  behave  as  classical  particles,  traveling 
along  straight  line  trajectories,  and  the  time-dependent  Schrodinger  equation  is  solved  to 
calculate  the  probability  of  various  types  of  electronic  transitions  (References  49-51).  Forcing  the 
particles  to  travel  along  straight  lines  limits  the  validity  of  the  impact  parameter  method  to 
collisions  of  several  hundred  electron-volts  oi  greater  (Reference  52).  The  wave  formulation  of  the 
method  of  perturbed  stationary  states  appears  to  be  one  practical  method  of  calculating  thermal 
energy  charge-exchange  cross  sections.  A  semi-classical  close-coupling  method  (Reference  53), 
based  on  an  averaging  scaling  procedure,  also  offers  utility  for  low  to  intermediate  collision 
energies. 

Approximate  quantal  treatments  of  reactive  charge  transfer  at  intermediate  energies  fall 
generally  into  three  basic  categories: 

(a)  Bom  approximation  ignoring  distortion 

(b)  Distorted  wave  Bom  approximations  (adiabatic  perturba¬ 
tions) 

(c)  Quantal  impulse  approximation 

(a)  In  the  Born  approximation,  distortion  of  the  target  by  the  projectile  ion  is  ig¬ 
nored  and  the  reactive  and  product  wavefunctions  are  replaced  by  their 
asymptotic  forms.  Since  no  explicit  mechanism  for  either  charge  transfer  or 
dissociation  is  possible  within  this  end-state  approximation,  the  method  is  not 
particularly  useful  for  ion-molecule  collisions  at  low  collision  energies.  The 
qualitative  behavior  of  the  Born  approximation  does,  however,  indicate  that 
the  cross  sections  should  decrease  very  rapidly  with  increasing  collision 
energy. 
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(b)  The  distorted  wave  Born  approximation  is  based  on  the  exact  solution  for  the 
reactive  cross  section  in  terms  of  the  T-matrix  elements: 


(68) 


where  Xfi  solution  of  the  Schrodinger  equation  using  incoming  boundary  conditions  and  Xa  is 
the  solution  of  the  complete  Hamiltonian  for  the  outgoing  reaction  channel,  including  all  possible 
internal  product  states.  V)?'  is  the  reaction  perturbation  potential  which  includes  distortion, 

polarization,  etc.  In  the  distorted  wave  Born  approximation,  Xfi  and  Xa  are  replaced  by  simpler 

approximations,  usually  corresponding  to  a  single  internal  state  only.  Choi  and  Tang  (References 
54  and  55)  have  discussed  the  numerical  problems  associated  with  evaluating  matrix  elements 
arising  from  this  approximation  and  several  applications  to  reactive  collisions  have  been  reported 
(References  56  and  57). 

A  still  simpler  approximation  is  available  by  reducing  the  matrix  elements  given  by  Equation 
(68)  to  a  product  of  an  electronic  matrix  element  coupling  the  reactant  and  product  states  and  a 
Franck-Condon  overlap  of  the  corresponding  nuclear  wavefunctions: 


T^  —  Mel  {xp\Xa') 


(69) 


This  so-called  overlap  or  Franck-Condon  model  was  introduced  by  Halavee  and  Shapiro 
(Reference  58)  and  has  been  described  in  detail  by  Schatz  and  Ross  (Reference  59). 

(c)  In  the  impulse  approximation,  the  operator  of  the  T-matrix  element  given  by 
Equation  (68)  is  replaced  by  an  approximation  which  describes  the  elastic 
scattering  corresponding  to  a  collision  between  the  projectile  ion  and  the  final 
free  particle  of  the  target  molecule.  This  approximation  is  equivalent  to  the 
assumption  that  the  collisions  between  the  particles  are  sudden  so  that  the  col¬ 
lision  time  is  much  shorter  than  the  vibrational  period  of  the  target  molecule. 
Henglein,  et  al.  (Reference  60)  have  argued  that  for  ion-molecule  collisions  of 
several  eV,  this  model,  also  called  the  spectator  stripping  mechanism,  should 
be  valid.  Bates,  et  al.  (Reference  61)  have  described  a  classical  treatment  of 
this  impulse  approximation. 

In  the  following  section,  a  very  powerful  method  for  solving  the  reactive  coupled-channel 
collision  problem  is  described.  This  approach  is  often  referred  to  as  the  R-matrbc  or  reactance 
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matrix  propagation  method.  It  generally  yields  a  fast,  efficient  .solution  and  automatically 
produces  unitarity  in  the  S-matrices.  By  freezing  open  channels  or  including  the  effects  of  closed 
channels,  several  levels  of  approximation  are  possible  within  the  same  computational  framework. 

2.7  R-Matrix  Propagator  Solution  for  Charge  Exchange  Reactions 

The  theory  of  heavy  particle  charge  exchange  is  complicated  by  the  fact  that,  in  general,  no 
simple  natural  coordinate  system  exists  for  which  a  uniform  description  of  reactant  and  product 
states  can  be  written.  Considering  the  reaction: 


H+  +  SiH4 

't 

D 

->  SiH3+ 

+  Hi 

g 

1 

L  J 

r 

r 

SiH,+ 


+  2H 


(70) 


the  internal  coordinates  describing  SiH4  cannot  be  uniformly  mass-scaled  to  represent  SiH3‘'' 
formed  as  a  product.  A  natural  coordinate  system  for  such  a  collision  has  been  described  by 
Marcus  (Reference  62),  but  this  yields  very  complicated  forms  for  the  kinetic  energy  operator  and 
presents  many  computational  difficulties. 

Recently,  Stechel,  et  al.  (Reference  63)  have  described  an  approach  in  which  two  separate 
nonorthogonal  frameworks  are  employed  to  describe  reactant  and  product  states.  The  total 
scattering  wavefunction  is  written  as: 


Ip  =  V'“(R>ea)  + 


(71) 


where  a  and  y  represent  reactant  and  product  states,  respectively.  Expanding  each  as  a  linear 
combination  of  wavefunctions  of  the  type: 

^  (72) 

j 

they  show  that  a  unique  connection  between  reactant  and  product  states  can  be  achieved, 
provided  the  basis  functions,  {0'(Pi)}  are  taken  to  be  nonorthogonal,  i.e.,  5“"''  =  <  ^  0. 
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We  note  that  Equation  (72),  in  molecular  language,  is  a  product  wavefunction  of  the 
valance-bond  type.  We  also  note  that  an  SCF  or  MCSCF  framework  for  Equation  (72)  may  not  be 
appropriate  for  treating  the  collisional  aspects  of  the  reactant  channel  due  to  the  orthogonality 
constraints  that  must  be  imposed  on  such  wavefunctions.  The  use  of  non-orthogonal  basis 
functions  is  simplified  through  the  application  of  projection  algebra  techniques  for  the  construc¬ 
tion  of  wavefunctions  of  overall  symmetry.  Efficient  computer  codes  for  implementing  such 
techniques  has  been  developed  in  this  Research  Center  over  the  past  years  and  can  be  utilized  in 
constructing  the  proper  asymptotic  reactant  and  product  state  channels. 

The  solution  of  the  reactive  charge  exchange  problem,  including  the  product  dissociation 
channels,  can  be  accomplished  using  a  direct  Numerov  integration  of  the  Schrodinger  equations 
arising  from  the  basis  representation  of  Equation  (72).  Alternatively,  the  R-matrix  propagator 
method,  originally  formulated  by  Diestler  (Reference  64),  can  be  employed.  This  approach  relates 
the  known  asymptotic  solution  in  the  reactant  channel  to  that  in  the  product  channel  through  a 
series  of  recursive  steps.  Defining  the  normal  derivative  to  the  surface  R  =  constant  by 
VnVKR.P).  we  seek  a  propagator  operator  of  the  form: 


V<R,p)  =  RV„v<R,p) 


(73) 


The  explicit  matrix  representation  of  this  operator  can  be  written  as; 


=  -  r  I  VnV' 


+  r  ^  VnV' 


(74) 
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Stechel,  et  al.  (Reference  63),  have  given  explicit  expressions  for  constructing  the  components  of 
this  operator,  both  in  the  asymptotic  and  overlap  regions  of  the  interaction.  The  initial  boundary 

condition  in  the  reactant  channel  is  R'  =  1/2, ,  where  A,  is  the  eigenvalue  of  the  reactant  target 

(in  this  case,  SiH4).  The  accuracy  in  developing  the  final  R-matrix  depends  on  the  chosen  step 
size,  h„  which  is  taken  as  an  adjustable  parameter  that  can  be  optimized  during  the  course  of  the 
propagator  solution. 
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We  have  recently  implemented  this  R-propagator  analysis  at  UTRL,  using  a  modified 
version  of  the  code  developed  by  Stechel,  et  al.  (Reference  63).  One  application  has  been  the 
low-energy  collisional  de-excitation  reaction: 


N 


N 


(76) 


which  occurs  along  a  curve-crossing  potential  surface.  The  results  of  this  study  have 

been  reported  (Reference  65).  A  similar  study  of  the  kinetics  of  the  charge  transfer  reaction, 

+  N2  -+  NO"’"  +  N,  has  been  carried  out  (Reference  66). 

In  summary,  we  have  found  that  a  description  of  the  energetics  of  charge  transfer  reactions 
in  a  valence-bond  framework,  coupled  with  an  analysis  of  the  kinetics  of  these  collisions  using  a 
non -orthogonal  basis  for  reactant  and  product  channels,  within  the  R-propagator  method  of 
solution,  represents  an  optimum  analysis  for  the  low  to  intermediate-energy  kinetics. 

2.8  Coilisionaily  Induced  Transitions 

Consider  the  charge  transfer  process: 


A+  +  B  -  C+  +  D 


(77) 


where  A'^  is  a  charged  atomic  or  molecular  ion  reactively  scattered  by  a  target  gas  B  to  yield 
products,  C'*'  +  D.  If  the  reactant  and  product  channels  can  be  separately  decoupled,  it  is  often 
possible  to  treat  Equation  (77)  as  a  pseudo-curve  crossing  problem.  For  near-adiabatic  collisions, 
this  type  of  system  can  be  handled  using  a  semi-classical  approach,  directly  solving  the  second 
order  differential  equations  which  couple  the  adiabatic  levels,  possibly  of  different  angular 
momentum. 

The  adiabatic  states  and  V'2(r,R)  are  eigenstates  of  the  non-relativistic 

Hamiltonian,  36 ,  for  a  fixed  internuclear  separation,  R.  We  have 


39 


at^i(r,R)  =  ei0i(r,R):  i  =  1,2 


(78) 


In  the  two-state  approximation  we  have,  using  impact  parameter  formulation. 


t/j  =  ci(t)^i(r,R)  e 


+  C2(t)^2(f.R)  e  Jo 


^1. 


cjdr 


(79) 


where  R  =  R(t)  and  is  defined  by  the  collision  trajectory. 

We  now  require  the  inner  product  of  ip  and  0,  to  vanish  over  the  electronic  coordinates. 


( (p2\%-ih^\ip)  =  0 


(80a) 


(80b) 
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Combining  Equations  (78),  (79),  and  (80),  we  have 


3ci  ,  ,  3^1 
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Converting  from  differentiation  with  respect  to  time,  to  velocity  and  angular  momentum,  we  have 


30i  _  3R  30i 
3t  "  "iT  3R 


h  3t 
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dVe 

dt 


Lr0i 


(82) 


where  is  the  angular  velocity  of  the  internuclear  axis  and  Lj  is  the  angular  momentum 
coupling  operator.  The  first  term  on  the  right  hand  side  of  Equation  (82)  leads  to  the  well  known 
Landau-Zener  (Reference  67)  solution  for  states  of  identical  molecular  symmetry.  For  such  cases 
the  second  term  in  Equation  (82)  vanishes. 

For  the  systems  under  study  here,  where  the  second  term  on  the  right  hand  side  of  Equation 
(82)  contributes  in  Equation  (81),  we  are  led  to  the  following  coupled  equations: 
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aci  .  i  ^ 


=  0 


(83a) 


^  <  <P2lLrl<p^  >  e 

dt  A  dt 


-ij  (ei-e2)dr 


=  0 


(83b) 


Assuming  a  linear  dependence  of  €  on  R  near  the  crossing  point  of  the  collision,  we  have 


£2-^1  =  b^(R-Rx)  =  at 


(84) 


d(€2-€l)  , 
dR 


(85) 


Assuming  also  that  <  ^i|Lt|<^2  >  =  <  Lt  >  is  essentially  constant  over  the  dominant  region  of 
the  collision,  we  have 
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(86a) 


(86b) 


Equations  (86a)  and  (86b)  have  been  solved  numerically  by  Russek  (Reference  68)  for  the  case  of 
curve-crossing  states  with  AA  =  ±  1.  For  small  values  of  angular  velocity,  Russek  shows  that 
Equation  (86)  reduces  to  a  standard  Landau-Zener  form.  For  large  velocities,  the  general  solution 
of  Equation  (86)  must  be  employed. 

When  more  than  one  product  channel  is  possible,  a  matrix  of  coupled  equations  of  the  type 
given  by  Equation  (83)  arises  and  a  multi-surface  mechanism  is  required  for  a  complete  descrip¬ 
tion  of  such  systems.  Often  a  single  reaction  channel  dominates  and  electron  or  H  transfer 


42 


occurs  with  nearly  unit  probability  along  a  single  reaction  surface.  For  such  systems  where  a  curve 
crossing  is  evident,  the  harpoon  (Reference  69)  or  spectator  stripping  model  (Reference  70)  is 

often  useful.  This  leads  to  the  simplest  reaction  cross  section,  o  =  jtR^,  where  R  is  the  effective 
separation  for  a  crossing  on  the  reaction  path.  In  a  next  higher  order  approximation,  the  reaction 
cross  section  is  given  by 


a  =  2jtj  Pij(b)bdb  (87) 

where  P,j  is  the  single  crossing  transition  probability  ( s  1.0)  which  can  often  be  approximated  by 

the  Landau  -  21ener  formula.  An  extension  of  this  simplified  treatment  to  include  internal 
vibrational  states  has  been  given  by  Bauer,  et  al.  (Reference  71). 

2.9  Classical  Ion-Molecule  Collisions:  The  Langevin-Gioumousis-Stevenson 
Reaction  Rate 

In  the  simplest  development  of  ion-molecule  kinetics,  it  is  often  assumed  that  the  reaction 
will  occur  at  the  first  collision  between  the  reactant  pair,  provided  that  the.e  is  at  least  one 
exothermic  reaction  cl  'el  available.  Simple  ion-molecule  reaction  rates  have  therefore  been 
equated  with  ion-molecule  collision  rates.  Experimentally,  many  ion-molecule  reactions  do  occur 
at  every  collision  between  the  ion-molecule  pair,  but  there  are  many  exceptions  to  this  simple 
case.  It  is  now  established  that  in  many  cases,  ion-molecule  collisions  may  be  unreactive  even  if 
an  exothermic  channel  is  available  (Reference  72).  Competing  reaction  channels  or  the  presence  of 
barriers  in  the  reaction  path  (ex;  O'*"  +  N2)  may  also  cause  a  drastic  reduction  in  the  collisional 
reaction  rate  (References  73  and  74). 

An  elementary  treatment  of  ion-molecule  reaction  kinetics  can  be  developed  from  a  classical 
description  of  the  scattering  cross  section.  For  scattering  by  a  spherical  potential,  the  reaction 
cross  section  can  be  written  as 


0(E)  =  271  f  P(b,E)bdb  (88) 

Jo 

where  b  is  the  classical  impact  parameter  and  P(b,  E)  is  the  probability  that  the  reaction  will 
occur  for  a  given  energy,  E,  and  impact  parameter,  b.  For  ion-molf’cule  reactions,  a  model  can  be 
constructed  based  on  the  assumption  that  for  values  of  the  impact  parameter  less  than  a  critical 
orbiting  parameter,  bg,  the  reaction  probability  can  be  taken  as  unity.  This  critical  orbiting  impact 
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parameter,  bo,  is  defined  so  that  ion  orbits  for  which  b  <  bo,  spiral  in  close  to  the  target  and  orbits 
for  which  b>bocome  no  closer  than  h^/  'll.  For  the  simplest  possible  case  of  a  spherical 

molecule  with  isotopic  polarizability  and  with  no  permanent  moment,  interacting  with  a  point 
charge,  the  interaction  potential  takes  the  form 


V(R)  = 


In  this  case  the  orbiting  parameter,  bo,  is  found  to  be 


bo 


4e^a 
.  - 


(89) 


(90) 


where  is  the  reduced  mass  of  the  colliding  pair  and  v  is  their  initial  relative  velocity.  Using  bo  as 
an  upper  integration  limit  in  Equation  (88),  we  find  for  the  reaction  cross  section: 


aiy)  -  nhl-  — 

V 


e^g 
L  n 


1/2 


(91) 


Equation  (91)  has  been  restated  by  Gioumousis  and  Stevenson  (Reference  75)  in  terms  of  the 
energy,  E,  of  the  primary  ion  of  mass,  m,  as 


o(E)  =  7t 


2e^am 

.”^e“J 


1/2 


(92) 


which  shows  an  inverse  square  root  behavior  of  the  reaction  cross  section  with  collision  energy. 
Assuming  a  Maxwellian  distribution  function,  f(v),  for  both  the  ions  and  neutral  molecules,  the 
rate  coefficient  corresponding  to  the  cross  section  given  by  Equation  (91)  becomes 


k  =  f(v)ovdv  =  2;r  ^  sec  (93) 

Equation  (93)  is  the  well-known  Langevin  formula  (Reference  76)  for  the  reaction  rate,  valid  for 
close  encounters  and  typically  of  the  order  of  10"^  cm^/sec.  The  validity  of  Equation  (93)  is 
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strengthened  for  systems  where  the  polarizability  of  the  molecule  is  large  and  the  reduced  mass  is 
small,  so  that  few  collisions  occur  outside  of  the  orbiting  radius.  The  Langevin  rate,  which  is 
predicted  to  be  independent  of  the  collision  energy,  is  thus  most  valid  for  collisions  with  small 
target  gases  such  as  He  +  H2.  For  large  target  molecules  and  ion-molecule  systems  with  many 
internal  degrees  of  freedom,  including  target  dissociation,  the  Langevin  rate  will  predictably  be 
too  large.  This  simple  model  also  is  incapable  of  predicting  the  product  distribution  of  an 
ion- molecule  reaction.  It  should  be  noted  that,  in  this  regard,  the  Langevin  formula  is  formally 
invalid  since  the  requirement  of  unitary  in  the  scattering  matrix  (principle  of  detailed  balancing)  is 
violated  for  this  classical  treatment. 


SECTION  3 


DISCUSSION  OF  RESULTS 


The  main  focus  of  our  current  studies  was  on  understanding  the  role  that  positive  ion- 
molecule  dissociative  charge  transfer  reactions  play  in  the  discharge  chemistry  pertinent  to 
deposition  of  amorphous  silicon  thin  films.  A  preliminary  analysis  of  the  production  of  negative 
silane  ions  by  dissociative  attachment  and  negative  ion-molecule  dissociative  charge  transfer 
reactions  was  also  undertaken.  A  summary  of  our  results  from  these  investigations  follows. 

3.1  Dissociative  Charge  Transfer 

3.1.1  Thermochemistry  of  SiHn  Species 

The  first  part  of  our  study  of  positive  ion  silane  reactions  was  concerned  with  an  examination 
of  the  electronic  structure,  predicted  equilibrium  geometries  and  thermochemistry  of  both  neutral 
and  ionic  silane  species.  After  some  computer  experimentation  with  basis  sets  and  the  appropriate 
level  of  theory,  it  was  decided  that  a  globally  accurate  analysis  of  the  reaction  surfaces  could  be 
accomplished  using  the  MP2/6-31  +  i-G**  level  of  theory  (References  22  and  77).  At  this  level,  a 
split  valence  descripi  jn  of  the  valence  electrons  is  employed ,  augmented  by  d-orbital  polarization 
functions  on  Si  and  p-orbital  polarization  functions  on  H.  In  addition,  diffuse  sp  shells  are  added 
to  give  a  more  balanced  treatment  to  negative  ions.  This  choice  of  basis  was  also  employed  by 
Curtiss  and  Pople  (Reference  78)  in  earlier  work  on  SiHjj  species.  The  Mpller-Plesset  (MP2)  level 
of  theory  is  not  really  adequate  for  an  accurate  thermodynamic  analysis  but  was  a  compromise 
choice  since  large  regions  of  the  reaction  surfaces  (multiple  geometries)  needed  to  be  examined.  In 
Table  1,  we  give  tlie  calculated  geometry  and  electronic  energy  for  species  of  the  silane/hydrogen 
system  at  the  MP2/6-31-t-Kj**  level  of  theory.  These  calculated  data  were  combined  in  isodesmic 
reaction  groups  to  estimate  heats  of  formation.  Our  resultant  enthalpies  at  0  K  are  given  in  Table  2 
where  we  have  also  included  the  known  thermochemistry  of  non-silane  species.  This  table  is  the 
basis  of  our  energetic  analysis  for  the  ion-molecule  reactions  in  the  silane/hydrogen  system. 

3.1.2  H+  +  SiH4 

Our  next  task  involved  the  generation  of  the  potential  energy  surface  for  the  ion-molecule 
reaction: 


H+  +  SiH4  ->  H2-i-SiH3+ 


(94) 
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Table  1.  Calculated  Energies  of  the  Silane/Hydrogen  System 


MP2/6-31++G**  Level  of  Theory 


Point 


Species 

State 

Group 

R(Si-H) 

Angle 

E(MP2) 

SiH4 

'Ai 

Td 

1.4719 

-291.353112 

SiHs 

1.4723 

107.67 

-290.711666 

SiH2 

lAi 

Civ 

1.5079 

92.55 

-  290.097098 

SiH2 

3Bi 

C2v 

1.4699 

118.27 

-  290.074758 

SiH 

2n 

Oov 

1.5144 

-  289.483907 

SiH4+ 

(not  bound  as  T^  or  C3V 

structure,  weakly  bound  in  C5) 

SiH3+ 

lAi* 

Dsh 

1.4550 

-  290.426400 

SiH3+ 

3E" 

D3h 

1.5349 

-290.228128 

SiH2+ 

2Ai 

C2v 

1.4684 

119.97 

-289.206614 

SiH+ 

1S+ 

C«V 

1.4945 

-  289.206614 

SiH3- 

lAi 

C3v 

1.5281 

120.72 

-  290.740824 

SiH2- 

2Bi 

C2v 

1.5353 

93.09 

-290.124127 

SiH- 

32:- 

Cc«v 

1.5441 

-289.517930 

Energies  are  in  hartrees  ( 1  h  =  27.21 161  eV),  geometries  are  in  Angstroms  and  degrees. 
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Table  2.  Enthalpies  for  Silane/  Hydrogen/Noble  Gas  Species 


Species 

AH°fo  (kcal/mol) 

AH°fo  (eV) 

Ip(eV) 

SiH4  [lAj] 

10.5 

0.455 

11.00 

SiHs  [2Ai] 

47.6 

2.064 

8.01 

SiH2[lAi] 

65.5 

2.845 

9.15 

SiH2  [3Bi] 

86.6 

3.755 

8.24 

SiH  [2n] 

89.6 

3.885 

7.91 

Si  [3P] 

106.6 

4.623 

8.151 

SiH4+  [2A'] 

264.2 

11.457 

SiH3+  [lA-] 

232.4 

10.078 

SiH2+  [2Ai] 

276.6 

11.995 

SiH+  [11+] 

272.0 

11.795 

Si+  [2P] 

294.5 

12.771 

H2  [‘Xg+J 

0.0 

0.0 

15.426 

H[2S] 

51.63 

2.239 

13.598 

H3+  [lA'J 

265.82 

11.527 

H2+[2Xg+] 

355.73 

15.426 

H+  - 

365.72 

15.838 

Ar[lS] 

0.0 

0.0 

15.759 

Ar+  [2p] 

363.41 

15.759 

He  [IS] 

0.0 

0.0 

24.587 

He+  [2S] 

566.99 

24.587 

ArH+  [11+] 

270.73 

11.740 

HeH+  [1E+] 

322.68 

13.993 
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over  a  wide  range  of  internuclear  geometries.  This  surface  was  first  examined  at  the  HF/6-31G* 
level  of  theory  and  subsequently  at  the  correlated  MP2/6-31G*  level.  The  most  interesting  feature 
of  the  global  surface  for  this  reaction  is  the  abstraction  of  an  H"  along  the  H+-H-SiH3  bond  axis. 

In  Figure  1  we  show  this  region  of  the  reaction  surface,  plotting  the  energy  of  the  reactants  as  a 
function  of  the  expanding  H-SiH3  bond,  with  the  H+-SiH3  distance  as  a  parameter.  The  Si-H 
bond  angles  and  the  three  Si-H  bonds  not  involved  in  the  H"*"  attack  were  optimized  for  each  fixed 
IH-H  SiH3  geometry.  Basically,  this  region  of  the  surface  illustrates  the  energy  change  as  the  H' 
is  stripped  from  SiH4  and  attaches  to  the  colliding  H+.  At  the  Hartree  Fock  (HF)  level.  Figure  1 
indicates  that  the  H"  transfer  is  endothermic  (transfer  does  not  occur)  until  the  H''’-SiH3  distance  of 
closest  approach,  b  ,  is  less  than  approximately  4.5  A.  For  collisions  with  b  <  4.5  A,  the 
probability  of  the  H'  abstraction  reaction  (94)  is  unity  since  the  surface  is  uniformly  attractive. 

The  H+  +  SiH4  reaction  surface  was  also  calculated  globally  at  the  MP2/6-31G*  level  of 
theory.  This  surface,  which  is  shown  in  Figure  2,  was  chosen  as  our  prototype  surface  for 
analyzing  the  validity  of  R-matrix  calculations.  The  surface  shown  in  Figure  2  is  an  extension  and 
refinement  of  the  HF  surface  shown  in  Figure  1.  The  kinetics  of  the  ion-molecule  reaction  given 
by  Equation  (94)  were  computed  using  the  minimum  energy  reaction  path  along  the  global  surface 
given  in  Figure  2.  The  intrinsic  reaction  path  hamiltonian  method,  as  described  by  Diestler 
(Reference  64)  was  utilized  for  most  of  our  kinetic  studies  of  positive  ion-molecule  reactions. 

In  order  to  study  this  same  reaction  without  imposing  high  collisional  symmetry,  we  have 
calculated  the  H+  +  SiH4  reaction  in  Cj  symmetry  where  only  a  reflection  plane  for  the  SiH5'‘' 
complex  is  imposed.  The  calculated  results  are  shown  in  Figure  3.  Here  the  strength  of  the  H" 
abstraction  reaction  is  clearly  illustrated.  In  reaction  sequence  c),  we  already  see  the  proton 
attraction  causing  a  rotation  of  the  silane  molecule  to  yield  a  shorter  H-H+  bond  distance.  By 
sequence  O7  the  H"  abstraction  has  occurred  and  the  resulting  SiH3‘'’  ion  is  collapsing  to  a  trigonal 
form.  No  constraints  at  all  were  placed  on  bond  lengths  or  bond  angles  in  this  reaction  sequence. 
.At  the  reaction  sequence  illustrated  in  f),  a  large  electron  density  has  suddenly  built  up  around  the 
reactant  ion.  This  resembles  the  harpooning  type  of  reaction  that  occurs  in  ionic-covalent  curve¬ 
crossing  situations  but  the  analogy  should  not  be  taken  too  literally  since  only  a  single  reaction 
surface  is  involved  in  this  H'*'  +  SiH4  collision.  This  reaction  is  dramatically  portrayed  in  Figure  4 
which  illustrates  the  electric  potential  contour  plots  for  this  same  sequence  of  geometries.  Both  the 
rotation  of  the  SiH4  molecule  in  sequence  c)  and  d),  and  the  H"  abstraction  by  sequence  f)  are 
identifiable  in  these  contour  plots.  The  rapid  concentration  of  contours  about  the  reactant  H'*'  ion 
illustrates  the  electron  charge  buildup.  We  find  no  evidence  of  a  direct  reaction  pathway  leading  to 
SiH2+  formation  in  Cj  symmetry. 

The  potential  energy  reaction  surfaces  for  the  various  product  channels  of  the  reaction,  H+  + 
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Figure  3.  Reaction  Sequence  for  H+  +  SiH4  ->  H-j  +  SiH 
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Figure  4.  Electric  Potential  (H+  +  SiH4  — >  H2  +  81113+). 


) 


Silij,  \vere  next  examined  using  correlated  wavefuncdons.  The  calculations  were  carried  out  at  the 
MP2  level  of  theory  using  a  Gaussian  6-31++G**  basis,  which  includes  both  polarization 
functions  and  diffuse  functions  on  silane  and  hydrogen  atoms.  The  product  channels  SiH+  +  H2 
and  SiH3+  are  shown  in  Figure  5  at  the  MP2  level.  We  find  a  potential  energy  ridge  of  ~  6.8  eV 
separating  these  product  channels,  similar  to  a  6.0  eV  barrier  found  at  the  HF  level  of  theory. 

SiH+  is  thus  not  an  important  product  for  collision  energies  less  than  ~  2.3  eV  [6.8  eV  (barrier)  - 
4.5  eV  (exothermicity)].  The  dissociation  channel,  SiH2'’'  Si+  +  H2,  is  shown  in  the  co.ntour 
plot  given  in  Figure  6.  An  activation  barrier  of  2.3  eV  is  found  at  the  MP2  level  of  theory,  similar 
to  a  3.0  eV  barrier  found  at  the  SCF  level.  In  agreement  with  the  SCF  results,  the  MP2  minimum 
energy  reaction  path  is  for  a  perpendicular  approach  of  H2  to  Si"*".  The  contours  are  somewhat 
smoother  at  the  MP2  level  with  the  plateau  of  the  reaction  surface  at  =  1.6  A,  Rsi-H2  " 

A.  Minimum  energy  reaction  paths  for  the  addition  of  H  to  Si"*",  SiH'*’  and  SiH2'’’  are  shown  in 
Figure  7.  There  are  no  apparent  activation  barriers  for  the  addition  (or  elimination)  of  H  atoms  to 
silane  ions. 

The  diabatic  potential  energy  surfaces  describing  both  the  reactant  and  product  channels  for 
reaction  (1)  were  next  calculated.  The  long  range  attractive  potential  for  the  reactant  channel  was  fit 
to  an  ion-induced  dipole  form,  V(R)  =  -  ot/2R^.  The  polarizability  of  silane  was  taken  to  be  4.62 
A^  ,  in  agreement  with  Haaland's  (Reference  10)  recommendation.  The  long  range  form  for  the 
product  channel  was  scaled  to  the  asymptotic  value  calculated  for  a  H2-SiH3'‘’  separation  of  5.5 
bohrs.  The  correct  long  range  form  should  be  : 


V(R)  =  - 


e^  K-«x) 
R'*  3 


P2  (cos0)- 


e  K+gx) 

2R^  3 


(95) 


where  ttn  and  aj_  are,  respectively,  the  parallel  and  perpendicular  components  of  the  polarizability 
tensor  for  the  H2  molecule.  Tire  averaging  of  the  anisotropy  of  the  polarizability  mainly  affects  the 
elastic  cross  section  for  H2  +  SiH3'*'  collisions  and  has  little  effect  on  the  reactant  channel.  The 
resultant  diabatic  curves  are  shown  in  Figure  8  where  Hy  is  the  coupling  element  connecting 
reactant  to  product  channels.  For  comparison,  the  calculated  global  adiabatic  surface  for  this 
reaction  is  shown  in  Figure  2.  Diagonalization  of  the  diabatic  potentials  shown  in  Figure  8  leads  to 
the  minimum  energy  pathway  representing  the  valley  on  the  global  adiabatic  surface. 

Reactive  scattering  cross  sections  were  calculated  for  the  reaction,  H+  -f  Silij  H2  + 
SiH3*,  using  the  R  matrix  method.  The  input  diabatic  interaction  potentials  illustrated  in  Figure  8 
were  employed.  Our  calculated  cross  sections,  as  a  function  of  center-of-mass  collision  energy, 
are  given  in  Table  3  and  illustrated  in  Figure  9.  We  find  that  the  calculated  cross  sections  generally 
follow  the  Langevin  rate  energy  dependence,  but  arc  10-40%  higher  in  ab.solute  value,  for  collision 
energies  0.1  eV.  Tliis  suggests  that  the  effective  reactive  collisional  diameter  is  underestimated 
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MP2/6-31++G**  level  of  theory 
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Figure  5.  Potential  Energy  Surface  for  SiH3+  ->  SiH+  +  H2. 


Figure  6.  Contour  Plot  for  Si+  +  H2 
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INTERNUCLEAR  SEPARATION.  R  (a.u.) 


Figure  8.  Diabatic  Potential  Energy  Curves 
for  the  Reaction:  Sillj  +  II’*'  ->  Sin3+ +  H2 


Table  3.  Calculated  Reaction  Cross  Section 


SiH4  +  H+  ->  H2  +  SiH3+ 

Energy  (eV)  Energy  (hartrees)  Jmax  Sigma  (A^) 


0.01 

0.00037 

70 

384.7 

0.03 

0.0011 

90 

274.5 

0.10 

0.0037 

150 

182.2 

0.30 

0.011 

210 

108.5 

1.00 

0.037 

330 

59.8 

3.00 

0.11 

530 

29.0 

10.00 

0.37 

890 

13.4 

30.00 

1.1 

1230 

8.2 

100.00 

3.7 

2270 

5.4 

i 
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Calculated  Cross  Section 


by  the  Langevin  model.  Further  studies  of  this  point,  including  sensitivity  of  the  R-matrix  results 
to  the  calculated  input  potentials,  are  suggested  as  possible  extensions  of  the  present  investigation. 

3.1.3  SiH4+  Structure 

As  an  addendum  to  our  studies  of  the  H+  +  SiH4  reaction  surfaces,  a  study  of  the  stability  of 
SiH4+  was  carried  out.  Haaland  and  Rahbee  (Reference  79)  report  negative  evidence  for  the  silane 
cation  using  Fourier  transform  mass  spectrometry.  Recently,  Pople  and  Curtiss  (Reference  80) 
have  reported  that  SiH4+  is  stable  as  a  donor-acceptor-complex  of  H2  with  SiH2‘‘‘  (Cs  structure)  at 
the  HF/6-31G*  level  of  theoiy.  Previous  theoretical  studies  have  indicated  no  stability  for  SiH4+ 
as  a  T^i  or  C3V  structure.  The  Cg  structure  reported  by  Pople  and  Curtiss  should  appear  at  1 1.22  eV 
which  is  in  excellent  agreement  with  the  result  of  1 1.0  eV  reported  by  Berkowitz,  Greene  and  Cho 
(Reference  81).  In  order  to  gain  more  insight  into  this  Cg  structure  for  SiH4+,  a  structure 
optimization  calculation  was  carried  out  at  the  MP2/6-31++G**  level  of  theory.  We  find  an 
optimized  structure  very  similar  to  that  reported  by  Pople  and  Curtiss  as  indicated  below. 


R(Si-Hi,2) 

=  1.459  A 

R(Si-H3) 

=  1.859  A 

R(Si-H4) 

=  I.8I9A 

a 

=  120.9  ° 

P 

=  89.3° 

Y 

=  24.1° 

£  =  -290.9592  8775  hartrees 
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We  find  that  this  Cg  structure  has  slightly  more  stability  (shorter  Si-H2  bonds)  at  the  MP2 
correlated  level  of  theory  relative  to  the  Hartree  Fock  study.  Our  calculated  themiodynamics  for 
this  ion  are: 


SiH4  — >  SiH4'*’  e 

Ip  =  10.7  eV 

(96) 

SiH4+  SiH2+  +  H2 

AH  =  +  10.8  kcal/mol 

(97) 

These  calculated  data  are  in  excellent  agreement  with  the  experiments  of  Berkowitz,  et  al. 
(Reference  10)  and  the  recent  study  by  Borlin,  Heins  and  Jungen  (Reference  82). 

3.1.4  Thermochemistry  of  SiF^  Species 

The  chemical  reactions  involving  SiPj,  and  SiFn+  species  are  computationally  more  difficult 
to  analyze  than  the  corresponding  hydride  systems  owing  to  their  more  complex  electronic 
structure.  After  some  experimentation  with  basis  sets,  it  was  decided  to  restrict  our  smdies  to  the 
use  of  the  6-3 IG*  split- valence  shell  +  polarization  basis  set  at  the  MP2  level  of  theory.  Although 
absolute  thermochemistry  is  not  well  described  at  this  level,  recent  studies  reported  by  Weber  and 
Armentrout  (Reference  83)  and  older  data  on  SiFn  species  reported  in  JANAF  (Reference  84) 
could  be  used  to  calibrate  our  data,  which  should  be  accurate  for  isodesmic  processes.  In  Table  4, 
we  give  the  calculated  geometry  and  electronic  energy  for  species  of  the  silane/fluorine  system  at 
the  MP2/6-31G*  level  of  theory.  These  calculated  data  were  combined  in  isodesmic  reaction 
groups  and  combined  with  best  estimate  experimental  values,  where  available,  to  predict  a 
consistent  set  of  heats  of  formation.  Our  resultant  enthalpies  at  0  K  are  given  in  Table  5.  This 
table  is  the  basis  for  our  energetic  analysis  for  ion-molecule  reactions  in  the  silane/fluorine  system. 

3.1.5  Positive  Ion-Molecule  Reactions 

The  thermochemical  values  for  the  SiH^  and  SiF„  systems  given  in  Tables  2  and  5, 
respectively,  can  be  employed  to  construct  the  energetics  of  the  product  channels  for  silane  ion- 
molecule  reactions.  These  data  are  collected  in  Table  6  which  illustrates  the  exothermic  character 
for  most  of  the  reaction  channels.  The  following  reactions  were  examined: 


H+  +  SiH4  (98) 

H3+  +  SiH4  (99) 

Ar^  +  SiH4  (100) 

He+  +  CH4  (101) 

H3+  +  SiF4  (102) 
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Table  4.  Calculated  Energies  of  the  Silane/Fluorine  System 


MP2/6-31G*  Level  of  Theory 


Point 


Species 

State 

Group 

R(Si-H) 

Angle 

E(MP2) 

SiF4 

lAi 

Td 

1.5810 

109.472 

-  687.740660 

SiFs 

^3v 

1.5993 

110.777 

-  587.983952 

SiF2 

lAi 

C2v 

1.6160 

100.894 

-488.321054 

SiF2 

3Bi 

C2v 

1.6174 

115.923 

-  488.206088 

SiF 

2n 

Coov 

1.6272 

-388.585225 

SiF4+ 

(not  bound  as  T^j  or  C3V 

structure,  weakly  bound  in 

Cs) 

SiFs-*- 

lAi- 

C»3h 

1.5410 

-587.664495 

SiF2+ 

2Ai 

C2v 

1.5557 

119.850 

-  487.946747 

SiF+ 

lz+ 

Qov 

1.5592 

-  388.353729 

SiF+ 

a3n 

Qov 

1.5775 

-388.166157 

Energies  are  in  hartrees  (Ih  =  27.21161  eV),  geometries  are  in  Angstroms  and  degrees. 
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Tables.  Enthalpies  For  Silane/Fluorine/Noble  Gas  Species 


Species 

AHfO  (kcal/mol) 

AHfO  (eV) 

yeV) 

SiF4[lAi] 

-  384.96 

-  16.68 

[15.7] 

SiF3  pAi] 

-239. 

-  10.4 

[9.3] 

SiF2[lAi] 

-  139. 

-6.0 

[10.7] 

SiF2[3Bi] 

SiF  [2n] 

-  13. 

-0.6 

7.54 

Si  [3P] 

106.6 

4.62 

8.15 

Si2  [%-] 

140.320 

6.03 

SiF4+  [2A'] 

[-24] 

[-1.0] 

SiF3+  [U'] 

[-  26.2] 

[- 1.1] 

SiF2+  [2Ai] 

109.0 

4.73 

SiF+  [1Z+] 

170.4 

7.39 

Si+  [2p] 

294.5 

12.77 

F2  V] 

0.0 

0.0 

15.43 

HF  [lZ+] 

-65.13 

-2.82 

16.06 

F[2P] 

18.47 

0.80 

17.43 

F2+ 

361.7 

15.68 

HF+  [2l+] 

305.2 

13.23 

F+  [3p] 

420.25 

18.22 

Ar[lS] 

0.0 

0.0 

15.76 

Ai+  [2p] 

363.41 

15.76 

He  [IS] 

0.0 

0.0 

24.59 

He+  [2S] 

566.99 

24.59 

Uncertain 
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Table  6.  Product  Channels  for  Silane  Ion-Molecule  Reactions 


t 


AH  (kcal/moll 


1) 


2) 


3) 


4) 

5) 

6) 


7) 


H+ 

+ 

SiH4 

0.0 

H2 

+ 

SiH3+ 

-143.8 

2H2 

+ 

SiH+ 

-103.7 

H 

+ 

SiH4+ 

-60.4 

Ho  +  H 

+ 

SiH2+ 

-48.0 

2H 

+ 

SiH3+ 

-40.5 

2H2  +  H 

+ 

Si+ 

-30.1 

H2+2H 

+ 

SiH+ 

-0.44 

H3+ 

+ 

SiH4 

-> 

0.0 

2H2 

+ 

SiH3+ 

-43.9 

3H2 

+ 

SiH+ 

-3.8 

He+ 

+ 

CH4 

0.0 

He 

+ 

CH4+ 

-276.0 

He  +  H 

+ 

CH3+ 

-236.9 

He  +  H2 

+ 

CH2+ 

-219.7 

He  +  2H2 

+ 

C+ 

-121.4 

He  +  2H 

+ 

CH2+ 

-116.4 

He  +  H2  +  H 

+ 

CH+ 

-112.4 

He  +  H2  +  2H 

+ 

C+ 

-18.1 

He+3H 

+ 

CH+ 

-9.1 

H3+ 

+ 

SiF4 

0.0 

HF  +  H2 

+ 

SiF3+ 

+27.5 

H+ 

+ 

SiF4 

0.0 

HF 

+ 

SiF3+ 

-72.4 

He+ 

+ 

SiH4 

-> 

0.0 

He 

+ 

SiH4+ 

-313.0 

He  +  H2 

+ 

SiH2+ 

-300.9 

He  +  H 

+ 

SiH3+ 

-293.5 

He  +  2H2 

+ 

Si+ 

-283.0 

He  +  H2  +  H 

+ 

SiH+ 

-253.4 

He  +  2H 

+ 

SiH2+ 

-197.6 

He  +  H2  +  2H 

+ 

Si+ 

-179.7 

He  +  3H 

+ 

SiH+ 

-150.1 

He  +  4H 

+ 

Si+ 

-76.5 

Ar+ 

+ 

SiH4 

-> 

0.0 

Ar 

+ 

SiH4+ 

-141.5 

Ar + H2 

+ 

SiH2+ 

-97.3 

Ar  +  H 

+ 

SiH3+ 

-89.9 

Ar  +  2H2 

+ 

Si+ 

-79.4 

Ar  +  H2  +  H 

+ 

SiH+ 

-49.8 
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In  Figure  9  we  illustrated  the  dissociative  charge  transfer  cross  section  for  reaction  (98)  calculated 
as  a  single  product  channel  problem.  Other  possible  product  channels  are  listed  in  Table  6.  Our 
previous  studies  of  the  potential  energy  surfaces  for  silane  ions  indicated  an  activation  barrier  of 
6.8  eV  for  the  decomposition  of  SiH3‘'’  to  SiH'*'  +  H2.  In  addition,  we  find  a  barrier  of  2.3  eV  for 
the  decomposition  of  SiH2'^  to  Si"*"  +  H2.  There  appeared  to  be  no  barrier  for  hydrogen  atom 
extraction  from  either  SiH3+  or  SiH2'''.  These  data  were  incorporated  into  our  R-matrix  code  for 
calculation  of  the  reaction  cross  sections.  The  multi-channel  data  for  reaction  (98)  are  shown  in 
Figure  10,  along  with  the  Lange vin  rate  for  this  reactant  pair.  We  find  that  product  channels  (2), 
(5)  and  (6)  all  exhibit  greatly  reduced  cross  sections  for  low  energy  collisions.  This  effect  is 
primarily  due  to  the  barriers  in  the  product  channel  for  these  three  branching  routes.  Note  that  at 
higher  energies  (>  10  eV),  all  of  the  cross  sections  appear  to  be  large  and  approach  the  Langevin 
limit.  There  is  some  scatter  at  these  higher  energies  which  we  believe  is  due  to  numerical 
instability  in  describing  the  high  frequency  oscillatory  behavior  of  the  outgoing  wave  packet. 

For  the  H3+  +  SiHL;  system,  an  unexpected  sensitivity  of  the  calculated  charge  transfer  cross 
sections  to  the  collision  energy  was  found.  The  potential  energy  along  the  minimum  reaction  path 
for  H3+  +  SiH4  2H2  +  SiH3+  is  shown  in  Figure  11.  Our  calculated  cross  sections  for  this 
reaction  are  displayed  in  Figure  12  where  we  see  both  quantum  structure  and  a  sharp  falloff  in  the 
cross  section  for  £<-011  <  2  eV.  The  quantum  structure  is  not  unexpected  and  was  also  found  to 
some  degree  in  our  earlier  studies  of  the  prototype  H+  +  SiH4  H2  +  SiH3+  reaction.  Several 
tests  were  carried  out,  however,  to  better  understand  the  falloff  in  the  calculated  cross  section  at 
low  collision  energies.  Since  the  energy  defect  for  H3+  +  Silij  ->  2H2  +  SiH3+  is  small  (AE  = 
-  43.9  kcal/mol),  we  suspected  that  there  was  a  phase  cancellation  in  the  reaction  product  channel. 
Close  inspection  of  the  potential  curves  shown  in  Figure  1 1  indicated  a  srong  dependence  of  the 
product  channel  to  the  repulsive  nature  of  Vf(R)  near  the  asymptotic  region  of  V,(R).  This 
dependence  is  clearly  illustrated  in  Figure  13  which  shows  that  this  falloff  in  the  calculated  reaction 
cross  sections  is  not  found  if  Vf(R)  is  shifted  downward  by  -  5  kcal/mol.  The  H+  +  SiH4  reaction 
is  also  shown  in  Figure  13  for  comparison.  We  believe  that  this  sensitivity  of  our  calculated  cross 
sections  to  the  repulsive  interaction  potential  is  peculiar  to  the  H3+  +  SiH4  system  which  exhibits  a 
small  energy  defect.  More  studies  are  required  to  further  test  this  hypothesis. 

The  Ar*’  SiH4  ion-molecule  reaction  (100)  is  especially  interesting  since  several 
experimental  groups  (Refs.  85-87)  have  reported  that  this  reaction  is  inefficient  at  low  collisional 
energies.  The  long-range  interaction  potentials  correlating  with  the  several  possible  product 
channels  are  shown  in  Figure  14.  W'e  have  examined  the  interaction  potential  as  a  function  of  SiH4 
orientation.  For  collision  along  the  SiH  bond,  reaction  occurs  as  a  charge  transfer  event  at  large  (~ 
10  A)  Ar*"  SiH4  separations.  For  collisions  along  the  H-Ar-H  bisector  the  SiH4  target  is  distorted 
to  Cg  symmetry  with  an  energy  barrier  of  ~  53  kcal/mol,  relative  to  its  equilibrium  geometry.  This 
orientation  results  in  a  very  small  reaction  rate  for  low  energy  (<  1  eV)  collisions,  in  agreement 
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Figure  10.  Calculated  Cross  Sections  for:  H+  +  SiH^.  — >  products 


+  SiH4 


(■n'B)  A6j0U3  leiiuajod 
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Figure  11.  Potential  Energy  Along  the  Intrinsic  Reaction  Path 
for  H^+  +  SiHd  2H2  +  SiH,+ 


Reaction  Cross  Section,  a(A) 
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with  experimental  observations.  However,  this  suggests  that  reaction  (99)  should  proceed 
efficiently  at  higher  collision  energies.  Clearly  more  work  is  required  to  understand  these  reaction 
systems.  A  model  for  the  Ar^  +  SiH4  reaction  proceeding  along  the  Ar-H-Si  bond  can  be 
developed  based  on  the  energetics  of  the  charge  transfer  processes. 


Ar+  +  SiH4 

[Ar+-H-SiH3]* 

(103) 

[Ar+-H-SiH3]* 

[Ar-H’^-SiH3]*-* 

(104) 

[Ar-H-^-SiH3]** 

->  [Ar-H-SiH3+]** 

(105) 

In  (103),  the  P  ■•♦■-SiH4  complex  is  formed  at  very  long  range:  R(Ar^-SiH4)  ~  lOA.  Step  (104) 
illustrates  charge  transfer  from  Ar*'  to  H.  Finally,  step  (105)  illustrates  charge  transfer  from  to 
SiH3.  A  key  feature  of  this  model  is  that  the  bond  ener^  _  of  ArH+  (144.3  kcal/mol)  is  equal  to  the 
proton  affinity  of  SiH3.  This  resonant  situation  exists  only  for  the  Ar'*’  +  SiH4  system.  The  model 
can  be  looked  at  as  a  sort  of  super-exchange  of  the  electron  fi'om  Ar+-SiH3+  to  Ar  +  SiH3+.  This 
picture  suggests  that  SiH3'''  may  be  readily  formed  in  collisions  at  this  orientation  but  with  such 
low  kinetic  energy  that  these  ions  may  not  be  efficiently  detected. 

A  different  explanation  is  offered  by  Fisher  and  Armentrout  (Reference  85)  and  by  Haaland 
(Reference  10)  who  suggest  that  there  is  a  poor  match  in  the  ionization  potential  of  Ar  with  either 
ionization  of  the  2t2  or  3ai  orbital  of  neutral  SiH4.  In  contrast,  they  sugge  •.  that  Xe'*  +  SiH4  and 
Kr^  -K  SiH4  reactions  should  proceed  with  high  efficiency  through  ionization  of  the  2t2  level  of 
SiH4,  leading  to  primarily  SiH3'*'  and  SiH2'’'  products.  Further,  the  reactions  of  Ne'*'  +  SiH4,  and 
He+  +  SiH^,  which  have  a  better  energy  match  through  ionization  of  the  3ai  orbital  of  SiH4,  are 
predicted  to  yield  primarily  SiH'^  and  Si'*’  as  products.  Although  these  arguments  correlate 
somewhat  with  the  observed  product  distributions,  there  is  no  real  justification  in  using  an  orbital 
energy  criteria  to  correlate  energy  resonance  in  a  system.  One  conclusion  of  this  orbital  energy 
resonance  model  would  be  that  SiH"^  and  Si"*"  ion  should  be  formed  with  lower  efficiency  in  all  of 
the  rare  gases,  as  compared  to  Ar,  since  the  energy  defect  between  the  3sli  orbital  of  SiH4  and  the 
rare  gas  ionization  potential  is  smallest  with  Ar.  Based  on  the  data  of  Fisher  and  Armentrout 
(Reference  85),  this  is  clearly  not  the  case,  suggesting  that  a  more  complete  picture  of  the  actual 
potential  energy  surfaces  is  required  to  explain  the  observed  anomalous  behavior  of  the  Ar’’  +  SiH4 
system.  In  addition,  the  positive  energy  dependence  for  the  observation  of  all  SiHn’’’  product  ions 
for  <;  )llisional  energies  >  ~  2.0  eV  h  s  no  explanation  within  the  Fisher  and  Armentrout  model. 

Foi  collisions  along  the  H-Ar-H  bisection,  we  find  that  the  SiH4  target  is  distorted  to  Cg 
sym.metry  with  a  large  energy  barrier,  relative  to  its  equilibrium  geometry.  The  resonant  model  of 
charge  transfer  described  above  does  not  appear  to  be  correct  for  this  orientation  and  long  range 
charge  transfer  does  not  occur.  Clearly,  a  more  detailed  study  of  the  several  potential  energy 
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reaction  surfaces  is  required  to  understand  thw  physics  of  the  Ar*"  +  SiH4  system.  At  this  point, 
there  is  no  satisfactory  explanation  for  the  observed  efficiency  of  this  reaction.  However,  further 
experimental  studies  are  suggested  to  test  whether  SiH3'*'  ions  may  be  formed,  but  with  very  low 
kinetic  energies. 

The  He'*'  +  CH4  ion-molecule  reaction  (reaction  6)  was  examined  in  detail.  Our  analysis  is 
based  on  studies  of  the  He"*"  -1-  CH4  reaction  surface.  Our  initial  studies  of  this  surface,  using  a  6- 
31 IG**  basis,  indicated  that  the  interaction  potential  is  very  flat  up  to  a  He-C  bond  distance  of  ~ 

5.4  A.  At  this  point  there  is  a  sudden  electron  jump  (harpooning)  from  CH4  to  the  He'*'  center 
leaving  a  CH4'*'  radical  to  relax  from  an  initial  CH4  geometry.  This  first  step  in  the  reaction  is 
illustrated  in  Figure  15  which  shows  the  sudden  change  in  the  electric  potential  as  the  electron 
transfers  from  CH4  to  He"*".  Tlie  subsequent  reaction  is  the  decomposition  of  CH4'''  to  form  CH2'*' 
H2.  This  occurs  smoothly  on  the  surface  of  CH4'’',  as  illustrated  in  Figure  16.  A  hydrogen 
molecule  is  released  directly  from  the  initial  geometry,  leaving  a  methylene  cation  as  the  basic 
charged  fragment.  The  electric  potential  at  various  geometries  along  this  decomposition  path  are 
illustrated  in  Figure  16. 

A  second  decomposition  pathway  for  CH4'‘'  exists  via  the  ^Aj  surface.  Initial  studies  indicate 
that  an  H2'*'  ion  is  first  formed  as  the  C-H  bond  length  increases,  followed  by  an  electron  jump  to 
again  form  CH2'*'  +  H2.  The  electric  potential  at  various  geometries  along  this  decomposition  path 
for  CH4'*'  is  illustrated  in  Figure  17.  We  find  that  a  methylene  fragment  appears  to  fomi  with  the 
simultaneous  expulsion  of  a  H2’*'  cation.  At  R  ~  2Re  for  CH2,  a  back  charge  transfer  occurs,  again 
yielding  the  thermodynamically  more  stable  CH2‘’'  +  H2  products.  The  calculated  reaction  rate  for 
He'*'  +  CH4  is  close  to  Langevin  with  predicted  products  of  ^2*”  (most  abundant)  and  with  CH'*' 
and  H'*'  resulting  from  back  charge  transfer  to  the  decomposing  methylene  radical.  This  analysis  is 
in  complete  agreement  with  the  experimental  data  of  Chatham,  Hils,  Robertson  and  Gallagher 
(Reference  87). 

Finally,  an  analysis  of  the  H'*'(H3'*')  +  SiF4  reaction  (reaction  102)  was  carried  out. 
Calculations  of  the  lowest  potential  energy  surface  for  H"*"  +  SiF4  clearly  illustrate  that  reaction 
(102)  proceeds  as  a  fluoride  anion  abstraction  reaction  with  either  FH  or  H3'*'  as  the  reactant  ion. 
This  observation  is  similar  to  that  for  the  H'*'  (H3'*')  -f  SiH4  system.  Only  two  channels  are  open 
for  the  H3'''  +  SiF4  reaction  below  50  kcal/mol  collision  energy.  There  is  no  apparent  barrier  for 
HF  fomation  and  the  reaction  rate  is  predicted  to  be  close  to  the  Langevin  rate.  The  dominant 
product  channel  is  predicted  to  be  SiF3'*'  +  HF. 
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3.2  SiHn+  +  SiH4  Clustering  Reactions 


In  addition  to  our  studies  of  silane  dissociative  charge  transfer  reactions,  an  analysis  of  the 
Si'*'  +  SiH4  clustering  reaction  was  initiated.  Mandich,  et  al.  (Reference  11)  have  experimentally 
studied  the  silicon  cluster  growth  reactions  starting  with  Si'*'  +  SiBLj.  They  find  that  the  addition 
sequence  stops  at  the  Si5Hio'’'  structure.  In  a  parallel  theoretical  study,  Raghavachari  (Reference 
12)  has  reported  ab  initio  calculations  on  these  silicon  growth  reactions.  The  structures  were 
optimized  at  the  HF/6-31G*  level  of  theory.  The  initial  ion-molecule  complex  that  is  formed  in  the 
Si"^  -t-  SiH4  reaction  is  reported  to  be  a  near  linear  attack  of  Si"*"  along  an  Si-H  bond  as  shown  in 
Figure  18.  The  bridge  structure,  also  shown  in  Figure  18,  is  reported  to  be  ~  1.5  kcal/mol  higher 
in  energy.  Tlie  barrier  between  these  two  structures  has  not  been  calculated  but  is  estimated  to  be 
low  owing  to  the  extended  Si''‘-SiH4  bond.  Either  of  these  structures  is  stabilized  by  the 
polarization  forces  of  Si'*'  on  SiHL^.  These  result  in  a  long-range  [-ct/R4],  attractive  force  balanced 
by  exchange  repulsion.  The  relative  stability  of  the  two  structures  shown  in  Figure  18  is  therefore 
critically  dependent  on  the  value  of  the  polarizability  of  SiH4.  Structure  2  should  be  favored  for  a 
larger  value  of  the  polarizability.  It  therefore  becomes  important  to  carefully  examine  the 
polarizability  of  SiH4. 

Haaland  (Reference  10)  has  pointed  out  that  the  calculated  polarizability  of  SiH4  is  too  small 
(3.60  A^)  even  at  the  PriP2/6-31++G**  level  of  theory.  The  best  experimental  estimate  for  a  is 
4.62  A^,  according  to  Haaland.  We  iiave  computed  the  polarizability  using  several  basis  sets  at  the 
MP2  level  of  theory  available  in  the  CADPAC  program.  The  results  are  summarized  in  Table  7 
which  illustrates  that  the  use  of  diffuse  basis  functions  on  both  Si  and  H  is  required  for  a 
quantitative  estimate  of  the  polarizability.  Both  the  MP2/6-31G**(sp,d,Sn)  and  MP2/6- 
31  lG**(sp,d,Sn)  calculations  yield  a  calculated  value  for  the  polarizability  within  90%  of 
Haaland's  estimate.  The  calculated  value  at  the  MP2/6-31G*  level  of  theory  (a  =  3.21  A^)  is 
significantly  smaller  and  therefore  underestimated  the  long  range  attractive  force  in  the  Si'*’  +  SiH4 
reaction.  This  adversely  affects  the  relative  energetics  of  structures  1  and  2  shown  in  Figure  18. 

To  assess  the  role  of  the  polarization  potential  in  determining  the  relative  stability  of  these 
loose  charge  transfer  complexes,  a  series  of  calculations  was  undertaken.  Using  the  6-31G**(sp, 
d,Sn)  basis,  optimized  UHF  calculations  were  carried  out  for  both  the  Cj  Si-H  bond  structure  and 
the  C2v  bridge  structure.  The  results  are  given  in  Table  8  which  again  indicate  that  the  calculated 
energies  are  too  close  at  the  UHF  level  to  choose  between  these  two  clustering  mechanisms. 
However,  all  calculations  at  a  correlated  level  indicate  that  the  C2v  bridge  structure  is  lower  in 
energy.  These  new  data  shed  doubt  on  the  suggested  Si2H4'''  equilibrium  structure  (Reference  12). 
Insertion  of  an  Si  atom  into  the  long  Si-H  bond  in  the  Cs  ion-molecule  structure  yields  the  (Cj) 
equilibrium  structure  shown  at  the  top  of  Figure  19.  However,  since  the  C2v  bridge  structure 
shown  at  the  bottom  of  Figure  19  is  lower  in  energy  at  long  Si-H2  bond  lengths,  this  may  actually 


77 


Table  7.  Calculated  MP2  Polarizabilities  of  SiH4 


Method 

E(hartrees) 

a(A3) 

MP2/6-31G* 

1.4829 

-291.3168497 

3.20 

MP2/6-31G** 

1.4726 

-291.3498595 

3.17 

MP2/6-31++G** 

1.4719 

-291.3531119 

3.59 

MP2/6-31G**(sp,d.Sn) 

1.4730 

-291.3619018 

4.12 

MP2/6-311++G** 

1.4737 

-291.4967731 

3.74 

MP2/6-311G**(sp,d.Sn) 

1.4759 

-291.5030978 

4.21 

EXP. 

4.62 

«sp  =  0.025795;  =  0.15; 

a.  =0.036 
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Table  8.  Energetics  of  the  Si'*'  +  SiH4  Clustering  Reaction 


Method 

Geometry 

Energy 

Rl 

R2 

% 

0 

Cg  Structure 

UHF/6-31G** 

1.460 

1.601 

1.802 

180. 

-  579.80980 

UHF/6-31G**(sp,d,Sn) 

1.461 

1.600 

1.784 

180. 

-579.81459 

MP2/6-31G** 

1.460 

1.588 

1.753 

180. 

-579.98130 

MP4/6-31G**//UHF 

-  580.02340 

C2v  Structure 

UHF/6-31G** 

1.457 

1.530 

1.971 

-  579.80650 

UHF/6-31G**(sp,d,Sn) 

1.457 

1.534 

1.975 

-579.81249 

MP2/6-31G** 

1.458 

1.537 

2.148 

-  579.98620 

MP4/6-31G**//UHF 

-  580.02720 

Bond  lengths  in  A,  bond  angles  in  degrees,  energy  in  hartrees. 


Bond  Structure 


Bridge  Structure 


,  Si+  -  SiH4  Ion-Molecule  Complexes. 


represent  the  more  stable  equilibrium  conformation.  The  bonding  in  higher  order  clusters  should 
now  be  reexamined  in  light  of  this  finding. 


3.3  e  +  SiH4  CoIIisional  Dissociative  Attachment 

Ab  initio  calculations  of  the  dissociative  reaction  surfaces  for  negative  ion  formation  in  silane 
were  undertaken.  In  Table  9,  we  show  the  possible  product  branching  for  e  +  SiH4.  Two 
possible  kinetic  routes  must  be  considered  for  conditions  that  exist  in  the  cathode  sheath  region  of  a 
silane  glow  discharge.  The  first  assumes  that  an  unstable  resonance  state  of  Silit"  is  formed  by 
primary  electron  capture.  Subsequently,  this  unstable  negative  ion  can  decay  into  the 
thermodynamically  accessible  product  channels  which  correspond  to  dissociation  products.  Table 
9  illustrates  that  reactions  (2)  -  (6)  fit  this  criteria  but  that  (7)  and  (8)  do  not.  Dissociation  of  SiH4‘ 
should  exhibit  thresholds  of  1.27, 2.44, 2.78  and  3.09  eV  for  formation  of  SiH2",  SiH3",  Si"  and 
H",  respectively,  provided  that  there  are  no  barriers  for  dissociation  in  the  exit  channels.  Although 
taken  with  limited  energy  resolution,  the  data  of  Potzinger  andLampe  (Reference  88)  and 
Ebinghaus,  Krauss,  Muller-Duys:  and  Neuert  (Reference  89)  appear  to  be  in  agreement  with 
.hese  predictions  for  formation  ot  oiH2"  and  SiH3",  although  a  unique  reaction  pathway  cannot  be 
determined  from  these  experiments.  Srivastava  (Reference  90)  has  also  observed  resonances  at  ~ 
2.5  eV  w  hich  he  attributes  to  Sin2"  and  SiH3"  formation.  The  process  of  negative  ion  formation 
following  decay  of  a  resonance  state  of  SiH4"  was  first  described  by  Haaland  (Reference  91). 

In  order  to  analyze  the  possible  dissociation  products  for  low  energy  collisions,  the  reaction 
surfaces  governing  the  dissociative  attachment  (DA)  of  electrons  to  SiH4  were  analyzed  using  ihe 
intrinsic  reaction  coordinate  procedure  of  Schmidt,  Gordon  and  Dupuis  (Reference  92).  In  Figure 
20,  we  show  the  energetics  of  the  reaction  coordinate  leading  to  dissociation  of  SiH4  to  SiH3  +  H. 
At  the  HF/6-31++G**  level  of  theory,  we  find  a  bond  dissociation  energy  of  3.27  eV.  At  the 
MP2/6-31++G**  level,  a  bond  energy  of  3.85  eV,  very  close  to  the  experimental  value,  is 
predicted.  Shown  on  Figure  20  are  data  for  the  SiH3"  +  H  and  SiH3  +  H"  reaction  surfaces.  We 
find  an  avoided  crossing  of  these  two  surfaces,  very  close  to  the  intemuclear  separation  where  the 
ion  curves  cross  tlie  neutral  Sill3  +  H  +  e  surface,  but  probably  in  the  region  where  the  ion 
potentials  are  complex.  Thus,  the  excited  state  ion  channel,  leading  to  SiH3  +  H"  may  depopulate 
through  this  avoided  crossing  leaving  SiH3"  +  H  as  the  only  significant  products  of  the  DA 
reaction.  Experimental  data  by  Haaland  (Reference  91)  and  by  Potzinger  and  Lampe  (Reference 
88)  indicate  that  H"  is  probably  not  formed  by  direct  DA  of  e  +  SiH4. 

To  facilitate  analysis  of  the  a.symptotic  behavior,  we  have  calculated  the  polarizability  of  SiH3 
at  the  HF/6-31++G**  level  and  find  a  =  3.95  A^.  Tlie  corresponding  value  of  the  polarizability  of 
SiH4  at  this  level  of  theory  is  3.25  A^.  The  leading  term  of  the  long-range  potential  of  SiH3  +  H" 
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Table  9.  Thermodynamics  of  Negative  Ion  Formation  in  Silane 


e  +  SiH4[lAi] 


Reaction 

SiH4-[2T2] 

^  SiH2-[2Bi]  +  H2[l2:g+] 

SiH3-[Ui]  +  HpS] 

Si-[4S]  +2H2[lSg+] 

SiH3[2Ai]  +  H-[lS3 

->  Si[3p]  +  H2[J£g+] 
+  H[2S]  +  H-[^S] 

SiH-[3S-]  +  H2[l2:g+] 

siH[2n]  +  H2[i2:g+] 


AH  (eV) 

(+3.98) 

(1) 

+1.266 

(2) 

+2.442 

(3) 

+2.783 

(4) 

+3.094 

(5) 

+3.269 

(6) 

+  H[2S] 

+4.392 

(7) 

+  H-[lS] 

+4.915 

(8) 
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and  SiH3*  +  H  can  be  calculated  according  to  the  fonn;  V(R-^)  =  -  a/2R4.  These  asymptotic 
curves  are  shown  on  Figure  20  which  clearly  indicate  that  an  avoided  crossing  between  these  two 
channels  must  occur. 

Preliminary  studies  of  dissociative  attachment  O^A)  of  electrons  to  silane  were  undertaken 
through  analysis  of  the  possible  product  channels  and  quantum  coupling  with  the  reactants.  The 
lowest  adiabatic  attachment  channel  correlates  to  SiH2"pB|]  +  H2[X^Zg+]  at  1.266  eV.  The 
SiH4*  resonance  state  corresponding  to  this  channel  appears  to  be  of  the  Feshbach  type  with  a 
small  capture  width  from  reactants.  In  addition,  as  studies  described  below  will  indicate,  there  is  a 
large  barrier  on  the  reaction  surface  leading  to  SiH2"  +  H2  products.  The  second  DA  channel 
corresponds  to  the  formation  of  SiH3‘(’Ai]  +  H[2S]  at  2.442  eV.  This  channel  corresponds  to  a 
shape  resonance  but  also  would  appear  to  have  a  small  DA  cross  section  for  vibrationally  cold 
SiH4.  Srivastava  (Reference  90)  finds  a  resonance  in  the  SiH3'  +  H  channel  at  ~  2.7  eV  but  this 
experimental  energy  appears  to  be  the  same  for  formation  of  Si",  SiH",  SiH2'  and  SiH3",  which 
should  exhibit  quite  different  threshold  behavior.  Recently,  Wan,  Moore  and  Tos.sell  (Reference 
93)  have  reported  electron  scattering  cross  sections  for  silane  in  the  range  0.2-»12.0  eV.  They 
observe  a  broad  maximum  in  the  total  electron  scattering  cross  section  at  ~  2.9  eV  which  they 
attribute  to  a  virtual  transition  to  the  lowest  unoccupied  o*(t2)  level  in  SiHj.  This  feature  has  a 
large  cross  section  (~  50  A‘)  and  may  be  related  to  the  features  seen  by  Srivastava  in  dissociative 
attachment  (DA)  at  ~  2.5  eV.  Based  on  these  known  low-energy  quantum  processes,  our 
conclusion  is  that  low-energy  direct  electron  dissociative  attachment  processes  should  have  low 
probability. 

The  second  kinetic  route  for  DA  of  e  +  Silij  was  examined  by  looking  at  the  excited  state 
surfaces  of  SiH4'.  The  DA  mechanism  for  this  route  may  be  the  initial  capture  of  an  electron  into  a 
high  l>ing  SiH4'  state  with  subsequent  statistical  branching  into  the  allowed  product  channels.  For 
a  mechanism  of  this  type,  the  product  distribution  is  uncorrelated  with  the  initial  dynamics  of  the 
reaction  channel  and  standard  DA  cross  section  codes  are  not  useful.  However,  the  data  of 
Srivastava  suggests  that  a  stochastic  formation  of  product  states,  all  formed  from  the  same  siate(.s) 
of  SiH4"  lying  at  8-9  eV  is  occurring. 

In  order  to  identify  the  excitation  channels  for  e  +  SiH4,  extensive  configuration  interaction 
(Cl)  studies  of  Sil^'  were  next  carried  out  to  identify  the  excited  state  reaction  surfaces  that  play  a 
role  in  dissociative  attachment  of  electrons.  In  Table  10,  we  show  the  e.\citation  spectrum  .‘or 
51114"  calculated  at  the  equilibrium  geometry  of  SiH4  (R  =  1.477  A).  The  basis  employed  in  these 
calculations  is  that  due  to  Handy  (basically  6-3 IG**  augmented  by  a  diffuse  sp  shell  (a  = 
0.025795)  and  a  diffuse  d  function  on  Si  (a  =  0.15).  This  basis  is  known,  from  our  previous 
studies,  to  give  an  accurate  representation  of  the  polarizability  of  SiH4. 

At  the  equilibrium  geometry  of  SiH4,  the  low-lying  excited  state  surfaces  arc  triply 
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Table  10.  Vertical  Excitation  Spectrum  for  SiH4" 


Eigenstate 

Symmetry  [C3v] 

AE(eV) 

Vo 

2Ai 

0.0 

Vl  (3) 

2Ai 

0,59 

V2(3) 

2A, 

2.15 

Vs  (3) 

2A2 

9.10 

(3)  triply  degenerate  in  Tj  symmetry 


degenerate  since  the  full  point  group  with  this  symmetry  is  Tj.  In  the  C3V  subgroup,  the  low-lying 
states  of  SiH4“  are  all  of  symmetry,  which  is  representative  of  a  loosely  bound  electron, 
symmetrically  distributed  about  SiH4.  These  are  pseudoeigenstates  of  the  system,  strongly 
dependent  on  the  chosen  basis  set,  and  are  not  physically  meaningful.  However,  the  third 
excitation  (at  9.20  eV)  is  characterized  by  a  compact  charge  distribution  resulting  from  s->p 
promotion  on  the  Si  atom.  We  identify  this  state  as  2A2  in  C3V  symmetry  (^A"  in  Q)  and  believe 
that  it  represents  the  first  identification  of  the  resonant  excited  states  in  this  system.  Our  location  of 
this  state  at  9.1  eV  is  in  good  agreement  with  previous  experimental  data  placing  the  important 
dissociative  attachment  states  in  the  7-10  eV  range. 

To  give  some  feel  for  the  magnitude  of  these  calculations,  our  configuration  interaction  (Cl) 
expansion  is  composed  of  all  single,  double,  triple  and  quartic  excitations  in  an  active  orbital  space 
of  the  lowest  16  natural  orbitals.  This  results  in  a  Cl  comprising  some  28000  configurations.  The 
identification  of  the  9.1  eV  state  as  2A2  suggests  that  dissociative  attachment  should  not  be  strong 
owing  to  symmetry  restrictions  which  weakly  couple  this  state  to  the  ^Aj  channel  corresponding  to 
e  +  SiH4.  Based  on  this  initial  study,  a  more  thorough  examination  of  excited  state  surfaces  of 
SiH4"  is  indicated. 

Alternate  kinetic  routes  for  the  production  of  negative  ions  of  silane  can  be  derived  from  the 
assumption  of  LTE  (or  at  least  quasi-steady  state)  in  the  cathode  sheath  region.  In  this  situation, 
electron  attachment  can  occur  not  only  to  SiH4,  but  also  possibly  directly  to  the  neutral  silane 
fragments  which  are  in  the  plasma.  In  Table  11,  we  illustrate  electron  attachment  to  SiH3,  SiH2 
and  SiH.  It  can  be  seen  that  a  wide  spectrum  of  threshold  energies  is  spanned  by  these  molecules 
(0.39  -  2.53  eV)  and  that  all  possible  negative  ion  species,  including  H’,  should  be  possible.  The 
relative  importance  of  these  several  processes  is,  of  course,  dependent  on  the  reactant  species 
concentration,  so  that  knowledge  of  the  neutral  composition  of  the  plasma  in  the  sheath  region 
becomes  a  requirement. 

3.4  Negative  lon-MoIecuIe  Silane  Reactions 
3.4.1  Thermochemistry  of  SiHn/SiHn"  Species 

Analysis  of  the  negative  ion-molecule  sik  ^e  reactions  required  development  of  a 
thermodynamic  data  base  for  these  anions.  Data  on  the  SiH^"  anions  have  been  reported  by 
Potzinger  and  Lampe  (Reference  88),  Drzaic,  Marks  and  Brauman  (Reference  94),  Mead,  Stevens 
and  Lineberger  (Reference  95)  and  Nimlos  and  Ellison  (Reference  96).  Using  these  data  as  a 
source  for  comparisons,  a  study  of  the  electronic  structure,  predicted  equilibrium  geometry  and 
resultant  calculated  thermochemistry  for  neutral  and  anion  SiH^  species  was  undertaken.  To 
establish  a  consistent  thermodynamic  base,  an  expanded  oasis  was  employed  at  the  MP2  level  of 
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Table  11.  Thermodynamics  of  Negative  Ion  Formation  in  SiH/SiH2/SiH3 


Reaction  AH  (eV) 


e  +  SiH3[2Ai] 

SiH3-[Ui] 

-1.406 

SiH-[3E-] 

+ 

+0.544 

SiH2-[2Bi] 

+ 

H[2S] 

+1.896 

-> 

SiH2[Ui] 

+ 

H-[lS] 

+2.266 

— ^ 

SiH[2n] 

+ 

H[2S]  +  H-[iS] 

+5.545 

e  +  SiH2[Ui] 

SiH2-[2Bi] 

-1.124 

Sr[4S] 

+ 

+0.393 

SiH-[3E-] 

+ 

H[2S] 

+2.002 

-4 

SiH[2n] 

+ 

H-[lS] 

+2.525 

— ^ 

Si[3p] 

+ 

H[2S]  +  H-[lS] 

+5.502 

e  +  SiHpH] 

— > 

SiH-[3E-] 

-1.277 

Si-[4S] 

+ 

H[2S] 

+1.592 

— > 

Si[3p] 

+ 

H-[lS] 

+2.223 
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theory.  The  silicon  basis  was  the  standard  6-3 IG*  approximation,  augmented  by  a  diffuse  (a  = 
0.025795)  sp-shell  and  a  second  diffuse  (a  =  0.15)  d-shell.  The  hydrogen  atoms  were 
represented  using  the  standard  6-3 IG*  basis,  augmented  by  a  diffuse  (a  =  0.036)  s- shell.  This 
basis  yields  a  good  approximation  for  the  polarizability  of  silane  (a  =  4.1 1  A^)  and  should  provide 
a  reasonable  compromise  for  calculating  the  thermodynamics  for  the  formation  of  SiH3",  SiH2" 
and  SiH“. 

In  Table  12,  we  give  the  calculated  geometry,  electronic  energy  and  predicted  electron  affinity 
for  species  of  the  SiHn/SiHn"  system  at  the  MP2/6-31  (2d,sp,Sn)  level  of  theory.  With  the 
exception  of  SiH3",  our  calculated  electron  affinities  are  in  reasonable  agreement  with  the  known 
experimental  data.  Additional  diffuse  functions  on  the  H  atoms  are  required  to  improve  the 
agreement  for  SiH3'. 

3.4.2  Thermochemistry  of  SiFn/SiFn“  Species 

A  study  of  the  electronic  structure  of  SiFn  and  SiFn",  similar  to  that  carried  out  for  the 
SiHn/SiHn"  system,  was  also  undertaken.  Owing  to  the  greater  complexity  of  this  system,  a  more 
restricted  basis  set  was  examined.  The  silicon  basis  was  the  standard  6-3 IG  split-valence  shell 
basis  augmented  by  d  polarization  functions  (a  =  0.45)  and  a  diffuse  (a  =  0.0331)  sp  shell.  The 
fluorine  atoms  were  also  represented  with  the  standard  6-3 IG*  basis  augmented  by  d  polarization 
functions  (a  =  0.8)  and  a  diffuse  (a  =  0.1076)  shell.  The  calculated  geometries  and  electron 
affinities  are  given  in  Table  13.  Although  F*  is  well  represented  by  this  basis,  it  can  be  seen,  that 
the  electron  affinity  of  both  Si  and  SiF3  are  underestimated.  Experimental  data  on  the  electron 
affinity  of  SiF3  have  been  reported  by  Richardson,  Stephenson  and  Brauman  (Reference  97).  The 
use  of  the  Handy/Amos  basis  for  Si,  which  was  chosen  for  our  studies  of  the  SiHn/SiHn"  system, 
would  yield  improved  agreement  but  with  a  severe  computational  penalty  owing  to  the  extended 
size  of  the  basis.  The  calculated  electron  affinities  in  Table  13  for  SiFp  species,  should  therefore 
be  mterpretcd  as  uniformly  low,  perhaps  by  20  -  25%.  It  is  interesting,  however,  that  all  SiFn" 
species,  with  the  exception  of  the  SiF4"  resonant  state,  are  predicted  to  be  thermodynamically 
stable. 

3.4.3  Silane  Decomposition  Mechanisms 

The  generation  of  the  potential  energy  surface  for  the  H“  +  SiH4  ion-molecule  reaction  was 
approached  in  two  steps.  Since  the  overall  symmetry  of  this  reaction  correlates  with  that  of  H  + 
SiH4,  our  initial  studies  dealt  with  analysis  of  the  thermal  decomposition  of  SiH4. 

Hydrogen  atom  abstraction  from  SiH4  occurs  with  no  decomposition  barrier  as  is  illustrated 
in  Figure  20.  Our  calculated  bond  dissociation  energy  at  the  MP2/6-31++G**  level  of  theory  is 
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Table  12.  Calculated  Energies  and  Geometries 
of  Neutral  and  Negative-Ion  SiKf^  Species 


MP2  Level  of  Theory 

Basis:  6-31G**:  Si:  sp  =  0.025795;  d  =  0.15 
H:  s  =  0.036 


Ea(eV) 

E  (hartrees)  _ 


Species 

State 

r(A) 

0{deg) 

M.P2 

Calc. 

Exp. 

H 

2S 

-  0.4988011 

H- 

1S 

-  0.5036267 

+  0.131 

+  0.754  a 

Si 

3p 

-288.8964559  (P) 

Si- 

43 

-288.9425435  (P) 

+1.254 

+1.385  a 

Ha 

'V 

0.7343 

-  1.1577722 

Ha- 

0.7634 

-  1.0849932 

SiH 

2n 

1.5166 

-289.4967536  (P) 

SiH- 

1A 

1.5385 

-  289.4980076 

+0.034 

+0.109  b 

SiH- 

3£- 

1.5443 

-289.5371537  (P) 

+1.099 

+1.277  b 

SiHa 

1A1 

1.5107 

92.108 

-  290.1083537 

SiHa" 

2Bi 

1 .5357 

92.905 

-290.1426034  (P) 

+0.932 

+1.124  b 

SiHs 

2Ai 

1.4739 

111.109 

-290.7213217  (P) 

SiHg- 

^Ai 

1.5296 

95.401 

-  290.7571167 

+0.974 

+1.406C 

SiH4 

^Ai 

1.4730 

109.471 

-  291.3619018 

SiH4- 

2Ai 

1 .4702 

109.471 

-291.3286672  (P) 

-  0.903 

a  ref.  95 
h  ref.  94 
c  ref.  96 
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Table  13.  Calculated  Energies  and  Geometries 
of  Neutral  and  Negative  Ion  SiF^,  Species 


MP2  Level  of  Theory;  Basis:  6-31 +G‘ 


E  (hartrees) 


Species 

State 

r(A) 

0(deg) 

MP2 

Calc. 

Exp. 

F 

2p 

-  99.5009122 

F- 

1S 

-  99.6260705 

+  3.406 

+  3.399  a 

Si 

3p 

-288.8847310  (P) 

Si- 

4S 

-  288.9233943  (P) 

+  1 .052 

+  1 .385  a 

SiF 

2n 

1.6496 

-388.6061626  (P) 

SiF- 

■I  A 

1.7177 

-388.5897919 

-  0.445 

SiF- 

3r 

1.7392 

-388.6293228  (P) 

+0.630 

SiFg 

'Ai 

1.6351 

100.391 

-  488.3506203 

SiFg- 

2Bi 

1.7184 

96.717 

-488.3573325  (P) 

+  0.183 

SiFg 

2Ai 

1.6101 

108.041 

-588.0224293  (P) 

SiFg- 

^Ai 

1 .7037 

96.739 

-  588.1057925 

+  2.268 

+  2.9b 

SiF4 

'Ai 

1.5892 

109.471 

-  687.7795548 

SiF4- 

2Ai 

1.6140 

109.471 

-687.7424847  (P) 

-  1.009 

(P)  Projected  UMP2  energy 
a  ref.  95 
b  ref.  97 
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3.85  eV,  in  good  agreement  with  the  experimental  value  of  3.90  eV.  Figure  20  illustrates  bond 
dissociation  at  the  HF/6-31++G**  level  where  a  somewhat  smaller  (3.27  eV)  bond  strength  is 
predicted.  Correlated  wavefunction  calculations  were  next  carried  out  for  the  concerted  C2v 
decomposition  of  SiH4  SiH2  + 112.  Since  the  SCF  representation  of  this  reaction  surface  is 
inadequate  for  description  of  the  bond  breaking  in  SiH4  and  formation  of  the  H2  bond,  a  correlated 
analysis  is  required.  In  Figure  21,  we  illustrate  the  decomposition  of  SiH4  in  C2v  symmetry.  The 
parameter  x  represents  the  Si  distance  from  the  center  of  mass  of  the  H2  molecule.  Clearly,  there 
are  two  separate  SCF  descriptions  at  extreme  values  for  x.  At  short  distances,  the  system 
approaches  the  Tj  symmetry  of  SiH4.  As  x  increases,  we  find  a  repulsive  barrier  of  ~  5  eV  before 
the  surface  takes  on  SiH2  +  H2  character.  Starting  from  SiH2  +  H2,  we  find  long-range  repulsion 
in  C2v  symmetry.  The  smooth  curve,  shown  in  Figure  21,  connecting  SiH4  SiH2  +  H2 
represents  a  two  configuration  treatment  of  Si-H  bond  breaking.  For  this  concerted  decomposition 
of  SiH4,  we  find  that  the  initial  ground  state  and  final  ground  state  correspond  to  different 
electronic  configurations  that  cross  with  an  energy  barrier  of  ~  5  eV.  Thus,  decomposition  in  C2v 
symmetry  is  not  a  viable  possibility. 

The  decomposition  of  SiH4  SiH2  +  H2  was  next  studied  in  symmetry  lower  than  C2v 
Gordon,  Geno,  Binkley  and  Frisch  (Reference  98)  has  reported  the  transition  state  for  the 
decomposition  of  SiH4  in  Cj  symmetry  at  several  levels  of  theory.  The  highest  level  of  theory 
used  was  MP2/6-31 IG  (2d,2p)  which  predicts  the  transition  state  shown  in  Figure  22.  Starting 
from  this  transition  state,  the  path  of  steepest  descent  leading  to  the  reactant,  SiH4,  or  the  products, 
SiH2  +  H2,  can  be  followed  analytically  along  the  intrinsic  reaction  coordinate  of  the  potential 
surface.  This  is  illustrated  in  Figure  23  which  shows  the  potential  energy  along  the  reaction 
coordinate.  Points  A,  B,  C  and  D  correspond  to :  [A]:  reactants  (SiH2  +  H2),  [B]:  a  local  van  der 
Waals  well,  [C].  the  transition  state  and  [D];  the  equilibrium  structure  of  SiH4.  The  transition  state 
is  predicted  to  lie  2.8  kcal.mol  above  the  dissociation  products,  SiH2  +  H2.  Thus,  the  minimum 
energy  required  for  this  decomposition  process  is  2.8  +  65.6  =  68.4  kcal/mol.  These  results  are  in 
good  agreement  with  the  study  reported  by  Gordon,  et  al  (Reference  98).  The  initial  SiH2'  +  H2 
pathway  is  also  shown  on  Figure  23.  This  curve  is  drawn  based  on  the  calculated  electron  affinity 
of  SiH2".  This  negarive  curve  must  cross  over  the  SiH2  +  H2  surface  at  some  intermediate 
geometry  since  SiH4"  is  not  a  stable  negative  ion. 

3.4.4  Negative  Ion-Molecule  Reactions:  H'  +  SiH4 

Several  studies  have  been  reported  which  examine  the  mechanisms  of  the  H"  +  SiH4 
reaction.  This  is  believed  to  be  a  classical  second  order  nucleophilic  substitution  reaction  (Sn2), 
similar  to  the  H"  +  CH4  and  F"  +  CH3F  systems  which  involve  substitution  at  the  carbon  center 
with  inversion  of  the  .MO  configuration.  Sheldon,  Hayes  and  Bowie  (Reference  99)  explored  the 
fully  rela.xed  motion  of  H'  +  SiH4  leading  to  three  possible  kinds  of  products,  a)  deprotonation  to 
form  H2  +  Sill3',  b)  equatorial  addition  to  fonn  SiH5“  and  c)  apical  addition  to  form  SiH5".  They 
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Figure  22.  Transition  State  for  SiH4  ->  SiH2  +  H2  Decomposition. 
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find  no  barriers  for  direct  attachment  of  H"  to  silane.  Once  formed,  however,  SiH5"  was  found  to 
exhibit  a  barrier  of  ~  70  kcal.mol  relative  to  SiH3"  +  H2  decomposition.  These  authors  fail  to 
report  their  calculated  stability  of  SiH5‘  relative  to  SiH4  +  H‘  reactants.  Earlier  studies  of  this 
reaction  have  been  reported  by  Wilhite  and  Spialter  (Reference  100),  Keil  and  Ahlrichs  (Reference 
101)  and  Baybutt  (Reference  102).  It  has  been  assumed  that,  the  Sjvi2-reaction  of  H"  +  SiH4 
proceeds  in  €3^  symmetry,  finally  forming  SiH5~  as  a  03^  (trigonal  bypyramid)  structure.  There 
has  been  uncertainty  about  the  stability  of  SiH5"  relative  to  H"  +  SiH4  and  some  evidence  of  a 
barrier  along  this  reaction  path.  Keil  and  Ahlrichs  (Reference  101)  and  Baybutt  (Reference  102) 
predict  no  biurier  for  this  reaction  whereas  Wilhite  and  Spialter  (Reference  100)  predict  a  large  (8.6 
kcal/mol)  barrier  along  the  5^2  reaction  path,  similar  to  that  found  for  the  H"  +  CH4  reaction. 

A  careful  study  of  the  H*  +  SiH4  reaction  shows  that  the  calculated  interaction  energy  is  very 
sensitive  to  the  chosen  basis  set.  In  the  conformational  change  from  H"  +  Silij  to  SiH5"  to  SiH3" 
+  H2,  the  diffuse  character  moves  from  the  hydrogen  center  onto  the  silicon  center.  Thus  a 
balanced  description  of  tiiis  reaction  requires  diffuse  polarization  functions  located  on  all  centers. 

In  addition,  a  diffuse  polarization  function  is  needed  on  the  hydrogen  centers  to  smoothly  permit 
charge  transfer  to  silicon.  In  Table  14  we  give  a  summary  of  the  energetics  of  this  reaction  with 
several  basis  sets  and  approximations.  We  find  that  SiH5"  is  stable  by  ~  20  kcal/mol  relative  to 
SiH4  +  H*,  in  good  agreement  with  the  most  accurate  studies  reported  in  Table  14.  However,  it  is 
clear  that  the  presence  of  a  large  (+39.4  kcal/mol)  reaction  barrier  for  formation  of  SiH3"  +  H2 
rules  out  the  Sfsj2  pathway.  The  energetics  of  the  H’  +  SiH4  reaction  are  illustrated  in  Figure  24. 

Brandemark  and  Siegbahn  (Reference  103)  and  Sheldon,  et  al  (Reference  99)  have  made  the 
interesting  suggestion  that  the  direct  proton  abstraction  reaction  may  be  the  preferred  pathway. 

H-  +  SiH4  [H--H+-SiH3-3  -+  H2  +  SiH3-  (106) 

They  find  a  much  lower  barrier  (tl7.0  kcal/mol)  for  this  direct  abstraction  reaction,  apparently 
ruling  out  the  formation  of  SiH5'  as  an  intermediate.  An  alternate  possible  route  is  to  abstract  H2 
from  the  C4V  pseudorotation  transition  state  geometry  which  lies  only  ~  2  kcal/mol  above  the  SiH5‘ 
minimum  energy  structure  (Reference  104).  At  the  suggestion  of  Mark  Gordon  (Reference  105), 
we  have  examined  dissociation  pathways  starting  from  the  C4V  transition  state  of  SiH5‘  which 
represents  the  intermediate  state  in  the  D3J,  D3h  pseudorotation  process.  As  mentioned, 

the  energy  b.urier  to  fomi  the  C4V  transition  state  structure  is  only  ~  2  kcal/mol,  indicating  that 
dissociation  from  this  geometry  may  be  energetically  possible.  Calculations  to  date,  however,  do 
not  indicate  a  low  barrier  reaction  pathway  from  the  C4V  transition  state.  Both  a  concerted  (C2v) 
and  relaxed  (C^)  dissociation  lead  to  high  energy  barriers  for  the  removal  of  H2.  We  conclude  that 
the  experimentally  observed  formation  of  SiH3'  must  arise  from  a  proton  transfer  reaction  along 
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Table  14.  Energetics  of  the  H  + SiH4  — >  SiH5  SiH3  +H2  Reaction 


AE  (kcal/mol) 


H"  -r  SiH4 

[SiHs  ]  j-gactants 
SiHs" 

(SiHslVoducts 

SiH3  +  H2 


a  h 

0.  0. 

0.  8.6 

- 13.8 

+  71.7 

-  3.1 


£  ^ 

0.  0. 

0.  0. 

-18.6  -20.3 


£  f 

0.  0. 

0.  0. 

-20.7  -10.2 

+  39.4 

-9.5  -4.1 


a)  Ref.  99 

b)  Ref.  100 
c;  Ref.  102 
(5  Ref.  101 
e)  Ref.  103 


HF/6-21G 

CNDO/2;  Basis;  Si;[12s9p2<l/634p2dl.  H:l6s/3sl 

HF/STO-3G  and  sc\cnu  larger  basis  seis  developed  from  Si(I0s6p)  and  H(5s)  primitive  sets. 
CEPA/largc  sp  basis  plus  2  sets  of  d-functions 
Cl/large  sp  basis  plus  2  sets  of  d-funettons 


I)  Ref.  10  MP4/6-31++G**//MP2/6-31-5-+G** 
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the  Si-H  bond  to  form  H2  +  SiH3",  rather  than  abstraction  of  H2  from  SiH5‘  along  the  8^2 
reaction  pathway. 

3.4.5  Negative  Ion-Molecule  Reactions  :  H"  +  SiF4 

The  SN2-reaction  of  H"  +  SiF4  has  been  studied  in  €3^  symmetry.  The  corresponding 
reaction  of  H“  +  SiH4  is  described  above.  For  that  system,  we  did  not  find  any  barrier  along  the 
reaction  path  in  the  entrance  channel.  This  result  contrasts  with  the  H"  +  CH4  system  which 
experimentally  exhibits  a  large  (~  8.6  kcal/mol)  barrier  toward  reaction.  In  Table  15  we  give  the 
enthalpies  for  the  silane/fluorine  system,  including  the  negative  ion  results  calculated  at  the  MP2/6- 
31+G*  level  of  theory.  It  can  be  seen  that  all  SiFj,  (n=l-3)  species  support  a  stable  negative  ion 
and  that  SiF3"  is  particularly  stable  with  an  electron  affinity  of  nearly  3  eV.  The  reaction: 

H- +  SiF4  ->  HF  +  SiF3-  (107) 

is  exothermic  by  19.3  kcal/mol.  At  the  MP2/6-314G*  level  of  theory,  we  find  an  exothermicity  of 
15.0  kcal/mol,  in  fair  agreement  with  experiment,  considering  the  limitations  of  the  basis  set. 

Calculations  of  the  H"  +  SiF4  reaction  surface  were  carried  out  in  €3^  symmetry.  Figure  25 
clearly  shows  that  the  reactants  and  products  are  represented  by  two  different  molecular  orbital 
descriptions,  characteristic  of  a  SN2-reaction  pathway.  The  repulsive  character  along  the  entrance 
channel  is  similar  to  that  found  for  the  H"  +  CH4  reaction.  The  barrier  shown  in  Figure  25  is 
clearly  too  high  since  a  two-configuration  MCSCF  or  GVB  description  of  this  system  is  required 
to  carry  the  description  smoothly  from  the  H"  +  SiF4  reactant  channel  to  HF  +  SiF3"  products. 

We  estimate  that  a  multi-configuration  treatment  will  lower  this  barrier  by  approximately  a  factor  of 
2,  yielding  a  somewhat  larger  barrier,  than  found  for  the  H'  +  CH4  system.  A  more  detailed 
investigation  of  this  system  is  indicated,  particularly  in  light  of  the  great  stability  of  the  SiF3'  ion. 
Based  on  the  results  of  this  initial  study,  we  predict  that  the  H"  +  SiF4  reaction  would  not  proceed 
for  collision  energies  below  ~  0.8  eV.  At  higher  collision  energies,  the  rate  should  approach 
Langevin  (k  ~  5xl0'9  cmVsec)  with  HF  +  SiF3”  as  the  dominant  products. 
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Table  15.  Enthalpies  for  Silane/Fluorine  Gas  Species 


Species 

AHj[0  (kcal/mol) 

AHfO  (eV) 

Ip(eV 

SiF4  [lAi] 

-  384.9 

- 16.68 

[15.7] 

SiF3  [2Ai] 

-239. 

-10.4 

[9.3] 

SiF3  [lAj] 

-  139. 

_6.0 

7.54 

SiF  [3n] 

106.6 

4.623 

8.15 

SiF4+  [^A'] 

[-  24.] 

[-1.0] 

SiF3+  [lA'] 

[-  26.3] 

[-  1.1] 

SiF2+  [2Ai] 

109.0 

4.73 

SiF+  [11+] 

170.4 

7.39 

Si+  [2P] 

294.5 

12.771 

F2 

0.0 

0.0 

15.686 

F2+  [2V] 

361.72 

15.686 

F  [2P] 

18.47 

0.801 

17.423 

F+  [3P] 

420.25 

18.224 

F-  [IS] 

-59.91 

-  2.598 

SiF-  [3Z-] 

[-25.] 

[-1.1] 

SiF2-  [2Bi] 

[-  143.] 

[-6.2] 

SiF3-  [lAi] 

-  306.7 

-13.3 

H-  [IS] 

32.34 

100 


Potential  Energy  (hartrees) 


Figure  25.  Reaction  Path  for  H"  +  SiF4 


101 


SECTION  4 


CONCLUSIONS  AND  RECOMMENDATIONS 

As  discussed  above,  several  qualitative  pictures  of  both  ion  and  neutral  fragment  production 
in  silane  discharges  have  been  developed.  Secondary  ion-silane  reactions  to  fonn  disilane  ions 
were  proposed  to  explain  the  observed  ion  currents.  A  separate  model  was  proposed  by  Turban,  et 
al  (Reference  7)  in  which  primary  silane  ions  are  formed  by  low-energy  electron  impact 
dissociation,  followed  by  secondary  ion-molecule  reactions.  Chapman  and  Gallagher  (Reference 
5),  however,  have  criticized  the  applicability  of  this  model  to  cathode  sheath  conditions  where  it  is 
believed  that  high  energy  electron  processes  (10-100  eV)  govern  the  overall  reaction  kinetics. 

Their  interpretation  of  the  discharge  ion  chemistry  is  that  electron  impact  dissociation  of  silane 
occurs  as  a  primary  process  forming  SiHn'*'  (n=0,3),  with  n=2  or  3  heavily  favored  over  more 
highly  stripped  silane.  Neutral  SiH^  radical  production  as  a  primary  process  must  also  be 
considered.  Based  on  the  known  rate  of  H  atom  abstraction  from  silane  and  the  direct  electron 
impact  dissociation  reaction,  SiH3  radicals  are  believed  to  predominate  in  the  discharge  region. 
SiH2'^,  and  more  highly  stripped  silane  cations,  are  probably  driven  to  SiH3  by  fast  reactions  of 
the  type  SiHn"*"  +  SiH4  SiHn+j  +  SiH3"^.  The  intent  of  this  technical  program  was  to  address 
these  uncertainties  in  the  chemistry  of  silane  ion-molecule  reactions  through  a  theoretical  analysis 
of  the  governing  reaction  rates.  Since  dissociative  reactions,  with  or  without  accompanying  charge 
transfer,  may  occur,  our  approach  was  based  on  detailed  quantum  mechanical  studies  of  the 
reaction  surfaces  to  yield  all  possible  reaction  product  channels.  Subsequent  to  this  analysis  of  the 
reaction  energies,  cross  sections  for  the  various  ion-molecule  processes  were  calculated  within  a 
quantum  framework. 

The  research  effort  carried  out  under  this  program  has  clearly  identified  several  areas  w'here 
further  theoretical  studies  should  be  carried  out.  These  areas  include  a  more  detailed  analysis  of 
ion-molecule  clustering  reactions  and  an  in-depth  study  of  the  role  that  negative  ions  play  in  silane 
plasmas.  The  premise  that  positive  ion-molecule  reactions  in  silane  plasmas  follow  the 
characteristic  growth  pattern,  forming  large  cluster  ions  which  eventually  form  observable  dust, 
has  been  questioned  by  both  experimental  (Reference  1 1)  and  theoretical  work  (Reference  12).  A 
further  argument  against  such  a  mechanism  is  that  positive  ions  have  short  residence  times  in  the 
discharge  since  they  are  swept  to  boundaries  of  the  reactor  by  the  ambipolar  sheath  fields 
(Reference  106).  Garscadden  (Reference  9),  Haaland  (Reference  107),  and  others  have  observed 
that  such  fields  act  to  constrain  negative  ions  in  the  homogeneous  reactor  region,  thereby 
amplifying  the  importance  of  their  gas  phase  chemistry. 

As  an  extension  of  the  research  problems  addressed  under  this  technical  program,  we  suggest 
that  additional  research  be  carried  out  to  gain  a  better  understanding  of  the  role  that  negative  silane 
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ions  play  in  semiconductor  deposition  discharges.  This  research  should  include;  1)  an  examination 
of  the  thermodynamics  of  SiH^",  Si2Hjj‘,  SiFjj"  and  Si2FmHn"  anions,  2)  studies  of  sequential 
clustering  reactions  of  silane  anions,  3)  the  role  of  dissociative  attachment  of  electrons  to  silane 
species  forming  radical  neutral  and  anion  products,  4)  studies  of  the  formation  and  decomposition 
pathways  for  negative  ion-molecule  reactions  such  as  H"  +  SiHjj  and  H“  +  SiFjj,  and  analysis  of 
the  energy  dependence  of  the  kinetic  rate  of  such  processes.  An  outline  of  our  recommendations  in 
each  of  these  areas  is  given  below: 

1 .  The  thermochemistry  of  neutral  SiH„  and  SiF^  species,  and  their  corresponding  cations  has 
been  the  subject  of  several  experimental  (References  108-113)  and  theoretical  studies  (References 
114-119).  In  contrast,  little  is  known  about  the  corresponding  anions,  particularly  the 
thermochemistry  of  SiFj,"  species.  Under  our  current  research  program,  an  initial  estimate  of  the 
anion  thermochemistry  was  made  at  the  MP2  level  of  theory.  Although  this  gives  a  consistent 
picture  of  the  relative  stability  of  these  species,  higher  order  theory,  such  as  the  MP4  or  QCISD,  is 
required  for  predictions  which  are  accurate  to  ±  1  kcal/mol. 

2 .  A  second  recommended  research  task  deals  with  the  analysis  of  sequential  clustering 
reactions  among  negative  silane  ions.  As  discussed  above,  Mandich  (Reference  11)  and 
Raghavachari  (Reference  12)  have  examined  the  Si+  +  SiH4  clustering  sequence,  leading  to 
successive  formation  of  Si2H2''’,  Si3H4+,  Si4H5+  and  Si5Hio+.  They  find  that  Si5Hio+  does  not 
undergo  H2  elimination  and  thereby  becomes  a  bottleneck  against  growth  of  larger  size  clusters. 
The  basic  chemistry  of  these  clustering  reactions  involved  addition  of  .he  SiH4  group  across  the 
Si-Si  bond  which  leads  to  branched  chain  products  rather  than  linear  species.  It  is  unclear  whether 
a  similar  reaction  pattern  would  be  found  for  negative  ion  clustering  where  the  relative  bond 
strengths  are  weaker.  It  is  also  unclear  whether  Si”  and  SijiiH,^”  clusters  would  attack  SiH4  along 
a  Si-H  bond  or  bridge  across  two  hydrogens  in  a  perpendicular  arrangement.  The  pattern  of 
growth  of  negative  ion  clusters  will  be  critically  dependent  on  this  mechanism  of  SiH4  attachment. 

3 .  A  third  recommended  task  is  to  examine  dissociative  attachment  of  e  -h  SiH4  in  greater  detail. 
Our  preliminary  studies  of  dissociative  attachment  (DA)  of  electrons  to  silane  were  concerned  with 
an  analysis  of  the  possible  product  channels  and  quantum  coupling  with  the  reactants.  The  lowest 
adiabatic  attachment  channel  correlates  to  SiH2”pBj]  -f  H2[XlZg+].  The  SiH4”  resonance  state 
corresponding  to  this  channel  appears  to  be  of  the  Feshbach  type  with  a  small  capture  width  from 
reactants.  In  addition,  there  is  a  large  barrier  on  the  reaction  surface  leading  to  SiH2”  +  H2 
products.  The  second  DA  channel  corresponds  to  the  fomiation  of  SiH3”[lAj]  +  H[2S].  This 
channel  corresponds  to  a  shape  resonance  but  also  would  appear  to  have  a  small  DA  cross  section 
for  vibrationally  cold  SiH4.  Based  on  these  known  low-energy  quantum  processes,  our 
conclusion  under  the  present  program  was  that  low-energy  direct  electron  dissociative  attachment 
processes  should  have  low  probability. 
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The  problem  of  DA  of  e  +  SiH4  at  higher  (>  4  eV)  collisional  energies  should  be  examined 
by  looking  at  excited  state  surfaces  of  SiH4".  In  this  case,  where  the  resonance  state  has  a  long 
lifetime,  a  two-step  mechanism  may  best  describe  the  attachment  process.  Step  one  involves  the 
capture  of  an  electron  into  a  resonance  state  of  the  negative  ion  with  lifetime  such  that  this 
resonance  state  has  forgotten  its  mode  of  formation  before  it  decays.  Step  two  governs  the  rate  of 
dissociation  of  the  resonance  state  into  the  possible  product  channels.  This  dissociation  can  occur 
along  d  number  of  exit  channels,  subject  to  the  usual  conservation  laws  of  total  energy  and  angular 
momentum.  The  competition  among  the  various  possible  modes  of  dissociation,  however,  does 
not  depend  on  the  manner  in  which  the  resonance  state  was  initially  formed.  We  suggest  an 
examination  of  this  model  in  detail  for  DA  of  e  SiH4,  based  on  an  analysis  of  the  location  and 
symmetries  of  the  excited  states  of  SiH4"  lying  in  the  7-10  eV  region  above  the  ground  state  of  e  + 
SiH4. 


4.  Finally,  we  suggest  a  study  of  the  formation  of  negative  silane  ions  through  ion- molecule 
reactions,  such  as  H"  +  SiH4  and  H"  +  SiF4.  Oui  initial  studies  of  the  H"  +  SiH4  reaction 
indicated  that  the  most  probable  reaction  pathway  ■  /as  abstraction  of  H'*'  from  SiH4  to  form  H2  + 
SiH3‘.  The  reaction  sequence  :  H"  +  SiH4  -4  SiH5"  ->  SiH3'  +  H2  in  contrast  to  the  H"  + 
CH4  system,  does  not  exhibit  a  barrier  to  the  form  ;tion  of  SiH5’.  However,  the  abstraction  of  H2 
occurs  with  a  57  kcal/mol  barrier,  making  the  overull  reaction  improbable.  In  the  H‘  +  SiF4 
reaction,  we  find  a  large  barrier  in  the  entrance  channel  leading  to  SiHF4",  similar  to  the 
corresponding  hydrocarbon  system.  The  topology  of  this  reaction  surface  is  more  complicated 
than  the  SiH5"  system  and  further  studies  of  the  potential  energy  surfaces  are  required. 

In  summary,  a  quantitative  description  of  all  of  the  major  reactions  that  occur  in  the  discharge 
region,  including  electron  impact  dis.sociation  and  ionization,  and  the  chemical  role  of  both  cations 
and  anions  (formed  by  dissvx:iative  attachment)  is  required  to  permit  reliable  modeling  of  .f 
discharge  sheaths.  A  goal  of  this  research  program  was  to  identify  the  positive  ion  chemistry 
occuning  in  silane  discharges,  and  to  predict  kinetic  reaction  rates  and  product  uistributions. 
Emphasis  was  placed  on  the  development  of  sound  the  ' .  t.eal  me  Jiods  to  predict  the  collisional 
energy  dependence  of  the  reaction  rates  and  on  examination  of  the  role  of  dissociative  charge 
transfer  in  silane  ion-molecule  reactions.  We  suggest  additional  studies  to  examine  the  role  of 
anion  reactions  in  both  charge  treuisfer  and  reactive  rearrangement  collisions.  An  overall  goal 
would  be  to  build  a  data  base,  con^ioting  of  both  theimochemical  and  reactive  cross  section 
calculations,  which  will  help  establish  the  role  that  negative  silane  and  fluonnated  silane  ions  play 
in  discharge  reactors. 
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APPENDIX  B-1 


POTENTIAL  ENERGY  SURFACES  FOR  SILANE 
ION-MOLECULE  REACTIONS 

H.  H.  Michels  and  R.  H.  Hobbs 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 

ABSTRACT 

Commercial  silane  discharge  reactors  usually  operate  with  a  mixture  of  silane  and  a  noble 
gas.  A  reaction  set  involving  collisions  between  hydrogen  ions,  noble  gas  ions  and  silane  (or 
disilane)  must  therefore  be  considered  for  a  complete  description  of  the  ion-molecule  chemistry  in 
such  devices.  Existing  experimental  data  suggest  that  no  simple  reaction  model  for  collisional 
dissociative  charge  transfer  can  be  inferred  for  this  system.  A  theoretical  quantum  mechanical 
study  has  therefore  been  initiated  to  define  the  pertinent  potential  energy  reaction  surfaces.  The 
definition  of  these  multidimensional  reaction  surfaces,  through  electronic  structure  calculations, 
should  yield  information  on  the  energetics  and  branching  probabilities  for  the  various  product 
channels. 

We  have  completed  a  study  of  the  reaction  surfaces  for  the  dissociative  charge  transfer 
reaction: 


H+  +  SiH4  SiH/  +  m  H2  +  ^  H 
(n=0,3;  n+^+2m=5) 

The  dominant  product  ion  is  found  to  be  SiH3'’',  formed  from  a  hydride  abstraction  mechanism. 
Dissociation  channels  connecting  to  SiH2'*'  and  Si"^  exhibit,  surprisingly,  significant  reaction 
barriers  relative  to  the  SiH3'*'  product  channel.  A  complete  set  of  the  potential  surfaces  for  this 
reaction  and  a  preliminary  analysis  of  the  kinetic  rates  and  product  distribution  are  described. 


Supported  in  part  by  WL  under  Contract  F33615-87-C-2718. 

Presented  at  the  Gordon  Research  Conference  on  Plasma  Chemistry,  Tilton,  NH,  1988. 
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APPENDIX  B-2 


DISSOCIATIVE  CHARGE  TRANSFER  IN  H+(H3+)  +  S!H4  REACTIONS 


H.  H.  Michels  and  R.  H.  Hobbs 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 


ABSTRACT 

An  ab  initio  study  of  the  potential  energy  reaction  surfaces  for  the  dissociative  charge 
transfer  reaction:  H+  +  SiH4  SiH^"^  +  mH2  +  £^2*  (n=0,3;  n+£+2m=5)  has  been  undertaken. 
The  basis  set  employed  was  the  standard  6-31G*  representation,  augmented  with  diffuse  functions 
(Si:3d,  H:2p)  for  an  accurate  description  of  the  polarizability  of  SiH4.  Dissociative  channels 
connecting  to  SiH2'^  and  SiH"*"  exhibit,  surprisingly,  significant  reaction  barriers  relative  to  the 
SiH3'''  product  channel.  The  dissociative  channel  leading  to  Si"*"  from  either  SiH2‘*'  or  SiH'*', 
however,  exhibits  no  reaction  barrier.  The  dominant  product  ion  is  found  to  be  SiH3+,  formed 
from  a  hydride  abstraction  mechanism.  A  preliminary  calculation  of  the  rate  for  SiH3'’’  formation, 
using  a  modified  R-matrix  propagator  technique,  yields  results  close  to  the  Langevin  prediction, 
Kl  =  5.2  X  10'^  cm^ec.  The  chemistry  of  the  H3'*'  +  SiH4  reaction  is  qualitatively  similar  with 
SiH3+  +  2H2  formed  as  primaiy  products. 


Supported  in  part  by  WL  under  Contract  F33615-87-C-2718. 

Presented  at  the  41st  Gaseous  Electronics  Conference,  Minneapolis,  Minnesota,  1988. 
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APPENDIX  3-3 


ELECTRONIC  STRUCTURE  AND  THERMOCHEMISTRY  OF  SILICON 
HYDRIDE  AND  SILICON  FLUORIDE  ANIONS 


i 


H.  H.  Michels  and  R.  H.  Hobbs 
United  Technologies  Research  Center 
East  Hartford,  CT  06108 


ABSTRACT 

An  ab  initio  study  of  the  electronic  structure  of  the  SiH^"  and  SiPp"  anions  has  been  carried 
out  using  M0ller-Plesset  perturbation  theory.  Optimized  geometries  were  calculated  at  the  SCF  and 
MP2  levels  of  theory  using  several  basis  sets.  Correlation  energy  treatments  included  both  MP4 
and  QCISD  methods.  The  anion  thermochemistry  was  calculated  through  several  isogyric 
processes  for  internal  consistency.  Excellent  agreement  with  experimental  electron  affinities  is 
found  for  the  SiHn"  anions.  We  find  that  all  SiFn"  anions  are  stable  with  the  exception  of  SiF4" 
which  exhibits  two  low-lying  resonance  states.  Corresponding  resonance  states  are  also  found  for 
the  SiH4"  anion.  These  states  adiabatically  correlate  to  SiH3"  +  H  and  SiH3  +  H"  products  but 
exhibit  a  potential  energy  curve  crossing  at  ~  2.1  A.  The  role  of  these  anions  in  plasma  chemical 
vapor  deposition  of  silane  and  fluorosilane  is  examined. 


Supported  in  part  by  WL  under  Contract  F33615'87-C-2718. 

Presented  at  the  Annual  American  Chemical  Society  Meeting,  Atlanta,  Georgia,  1991. 
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DISSOCIATIVE  ELECTRON  ATTACHMENT  OF  e  +  SiH4 


H.H.  Michels  and  R.  H.  Hobbs 

i  United  Technologies  Research  Center 

East  Hartford,  CT  06108 


J.  M.  Wadehra 

Department  of  Physics  and  Astronomy 
Wayne  State  University 
Detroit,  Michigan  48202 


ABSTR.1CT 


Dissociative  electron  attachment  of  e  +  SiHj  has  been  studied  as  a  mechanism  for  forming  the 
negative  fragment  ions:  SiH3~,  SiH2“,  SiH“  and  Si“  in  silane  discharge  plasmas.  Detailed 
calculations  of  the  low-lying  dissociative  reaction  surfaces  for  SiH4~  indicate  that  low-energy 
direct  attachment  processes  should  I.  ve  low  probability.  At  higher  electron  collision  energies  (7- 
10  eV),  a  long-lived  resonance  state  of  SiH4~,  which  statistically  decomposes  into  the  accessible 
negative  fragment  ions,  appears  to  be  the  dominant  attachment  process. 


« 


Supported  in  part  by  WL  under  Contract  F33615-87-C-2718  and  by  AFOSR  under  Contract 
F49620-89-C-0019  and  Grant  AFOSR-87-0342. 

To  be  submitted  to  Chemical  Physics  Letters,  1991. 
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1.  Introduction 

Dissociative  electron  attachment  (DA)  of  e  +  SiH4  is  being  studied  as  a  mechanism  for  forming 
the  negative  fragment  ions:  SiH^' ,  SiH2~,  SiH~  and  Si“  in  silane  discharge  plasmas.  Owing  to 
the  constraining  effects  of  the  sheath  fields  of  typical  plasma  reactors,  it  is  believed  that  clustering 
of  negative  ions  may  occur  in  such  devices.  In  contrast,  positive  ions  typically  have  a  short 
residence  time  in  the  discharge  region,  leaving  some  question  to  their  importance  as  nucleation 
centers.  ‘ 

The  first  measurements  cf  DA  in  silane  are  those  of  Ebinghausen,  et  al  [1]  who  report  the 
formation  of  SiH3“,  SiH2“  SiH“  and  Si~  as  products,  all  with  appearance  potentials  between  7  / 

and  8  eV.  Similar  results  were  found  by  Potzinger  and  Lampe  [2]  who  first  suggested  that  the 
observed  negative  ion  products,  SiH3~  and  SiH2“,  may  correspond  to  decay  channels  from  a 
single  resonant  stnicture  of  3iH4~  lying  at  7-8  eV.  More  recently,  Srivastava  [3]  has  reponed  DA 
cross  sections  for  the  fomiation  of  Sin3~,  SiH2“  SiH"  and  Si~,  all  peaking  in  the  8-9  eV  region. 

Potzinger  and  Lampe  [2]  and  Srivastava  [3]  also  report  the  appearance  of  a  much  weaker 
attaci'.ment  process  at  2-3  eV  for  the  formation  of  SiH3“  and  SiPl2~.  In  the  present  study,  the  low- 
lying  dissociative  reaction  surfaces  for  Sill4~  have  been  examined  in  an  effort  to  identify  the 
dissociative  attachment  mechanism. 

2.  Theoretical  Considerations 

The  traditional  DA  process  occurs  by  the  following  mechanism: 

e  +  AB  [AB**]-  A-<**t)  +  b(*.-)  (1) 

where  AB  is  any  diatomic  or  polyatomic  molecule,  [AB**]  is  a  resonant  state  of  the  corresponding 
negative  ion  and  and  represent  the  possible  set  of  dissociation  products,  possibly 

with  internal  eleclronic^'^^  or  lOvibrationaK”)  excitation.  The  cross  section  for  such  an  attachment 
process  is  governed  by  th,  chtiracteristics  of  the  resonant  state  [AB**].  In  resonant  scattering 
theory  [4],  assuming  a  local  width  approximation,  the  scatteri.ng  wovefunction,  ^(R),  satisfies: 

{ -  i  V2  +  jg+n  +  v-(R)  -  i  r(R)  -  e  ]  ^(R)  =  -v(R)  Cvj(R)  (2) 

2|i  2iiR2  2 

where  Cvj('f^)  *•'»  nucle:ir  motion  wavefunciion  of  the  neutral  target  molecule,  E  is  the  total 
(conserved)  energy  of  the  system,  and  V“(R)  -  i/2  r(R)  is  the  complex  interaction  potential.  The 
envelope  of  l^(R)|2  decreases  with  R(t),  owing  to  the  possibility  of  autodetachment,  and  the 
probability  for  dissociative  attachment  is  thus  governed  by  the  asymptotic  limit  of  l^(R)|2 : 

=  (3) 

p.  k  R-»« 
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where  k2/2  is  the  energy  of  the  incident  electron  and  is  the  relative  kinetic  energy  of  the  ion- 
atom  or  ion-molecule  product  state.  Eq.  (2)  is  developed  on  the  basis  that  the  outgoing  wavepacket 
is  directly  coupled  to  the  incident  channel  comprised  of  a  free  electron  and  the  neutral  target 
molecule.  In  such  a  case,  DA  occurs  as  a  sequence  of  dissociation  reactions,  each  process 
occurring  with  maximum  cross  section  at  the  energy  threshold  corresponding  to  the  opening  of  a 
i  particular  product  channel. 

In  the  case  where  the  resonant  state  has  a  long  lifetime,  a  two-step  mechanism  may  best 
^  describe  the  attachment  process: 

e  +  AB  <=>  [AB**]-  :  <\Kab°  Vab**>  (4) 

[AB**]~  A“^*>')  +  B(*d)  :  (5) 

Eq  (4)  governs  the  rate  of  capture  into  the  resonant  state  of  the  negative  ion  and  the  corresponding 
rate  of  autoionization.  Eq.  (5)  governs  the  rate  of  dissociation  of  the  resonant  state  into  the 
possible  product  channels.  This  dissociation  can  occur  along  a  number  of  exit  channels,  subject  to 
the  usual  conservation  laws  of  total  energy  and  angular  momentum.  The  competition  among  the 
various  possible  modes  of  dissociation,  however,  does  not  depend  on  the  manner  in  which  the 
resonant  state  was  initially  formed.  It  follows  that  the  overall  attachment  probability  can  be  written 
as  a  product  of  two  terms: 


W[t}/’j(£j^),  \|/{jj]  =  Wpgp[\}/j(ej^)]  WjJgpgy  [X/ccl 

(6) 

(7) 

WdccaytVa'i  “  !  ^lot  > 

(8) 

where  is  the  energy  corresponding  to  a  particular  resonant  state  in  the  incident  scattering 
channel,  E^gp  is  the  energy  width  for  coupling  the  incident  channel  into  the  resonant  state,  and 
is  the  partial  width  corresponding  to  the  probability  that  the  resonant  state  will  dissociate  into 
products  a.  Obviously,  the  total  decay  width,  Hot  =  S  Fjj,  is  the  sum  of  the  partial  widths  for  all 
open  dissociation  channels.  This  mechanism  is  often  described  in  the  literature  of  nuclear  physics, 
where  the  resonant  state  is  referred  to  as  a  "compound  nucleus"  [5,6].  The  independent 
probabilities  for  the  fonnation  and  subsequent  decay  of  the  resonant  state  are  thus  interpreted  as  a 
model  whereby  the  resonant  state  has  "forgotten"  its  mode  of  formation  before  it  decays. 

This  model  for  DA  can  also  be  understood  in  terms  of  a  semiclassical  analysis  of  scattering 
theory  [7,8].  In  such  a  description,  the  cross  section  for  DA  can  be  written  as  a  product  of  two 
factors: 


C-3 


(9) 


where  o^^gpCe)  is  the  cross  section  for  the  capture  of  an  electron  to  form  the  resonant  negative  ion 
state.  S  is  the  survival  probability,  which  is  classically  related  to  the  decay  width  of  the  resonant 
state: 

i 

S  =  exp(-l^s  r(R)dR)  (10) 

r 

where  the  region  of  interaction  extends  from  the  classical  capture  radius,  R(.,  to  the  effective  radius 
of  stabilization  of  the  ion-atom  or  ion-molecule  products,  for  which  V'CRj)  =  VO(Rs).  Thus  the 
resonant  state  decays  exponetially  with  mean  life  l/T. 


3.  Results  and  Discussion 

Calculations  of  the  dissociative  reaction  surfaces  for  negative  ion  formation  in  silane  were 
undertaken.  In  Table  1,  we  show  the  possible  product  branching  fore  +  SiH4.  Two  possible 
kinetic  routes  must  be  considered  for  conditions  that  exist  in  the  cathode  sheath  region  of  a  silane 
glow  discharge.  The  first  assumes  that  a  Jahn-Teller  unstable  PT2]  resonance  state  of  SiH4"  is 
formed  by  primary  electron  capture.  Subsequently,  this  unstable  negative  ion  can  decay  into  the 
thermodynamically  accessible  product  channels  which  correspond  to  dissociation  products.  Table 
1  illustrates  that  reactions  (2)  -  (6)  fit  this  criteria  but  that  (7)  and  (8)  do  not.  Dissociation  of 
SiH4"[-T2l  should  exhibit  thresholds  of  1.27, 2.44, 2.78  and  3.09  eV  for  formation  of  SiH2~, 
SiH3“,  Si“  and  H“  respectively.  Although  taken  with  limited  energy  resolution,  the  data  of 
Poizingcr  and  Lampe  [2]  appears  to  be  in  agreement  with  these  predictions  for  formation  of  SiH2“ 
and  Sin3~  Srivastava  [3]  has  also  observed  resonances  at  ~2.5  eV  which  he  attributes  to  SiH2“ 
and  SiH3~  formation.  This  process  of  negative  ion  formation  following  decay  of  the  2T2 
resonance  state  of  SiH4~  was  first  described  by  Haaland  [9]. 

Tlie  reaction  surfaces  governing  the  low-energy  dissociative  attachment  (DA)  of  elecU'ons  to 
SiH4  have  been  analyzed  using  the  intrinsic  reaction  coordinate  procedure.  In  Fig.  1,  we  show  the 
energetics  of  the  reaction  coordinate  leading  to  dissociation  of  SiH4  to  SiH3  +  H.  At  the  HF/6- 
31+-fG**  level  of  theory,  we  find  a  bond  dissociation  energy  of  3.27  cV.  At  the  MP2/6-31-t-+G** 
level,  a  bond  energy  of  3.85  eV,  very  close  to  the  experimental  value,  is  predicted.  Also  shown  on 
Fig.l  are  preliminaiy  data  on  the  51113"  +  H  and  SiH3  +  H"  reaction  surfaces.  We  find  an  avoided 
crossing  of  these  two  surfaces,  very  close  to  the  intcmuclear  separation  where  the  ion  curves  cross 
the  neuual  Sill3  + 11  +  e  surface,  but  probably  in  the  region  where  the  ion  potentials  are  complex. 
Thus,  the  excited  state  ion  channel,  leading  to  Silly  +  H~  may  depopulate  through  this  avoided 
crossing  leaving  Silly"  +  II  as  the  only  significant  products  of  the  DA  reaction.  Experimental  data 
by  Haaland  [9]  and  by  Potzinger  and  Lampe  [2]  indicate  that  H"  is  probably  not  fomied  by  direct 
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DA  of  e  +  SiFi4.  These  dissociative  attachment  studies,  including  the  work  reported  by 
Ebinghaus,  et  al  [  1 1,  all  indicate  that  low  energy  (<  3  eV)  DA  to  SiH}  must  occur  with  a  small 
cross  section  (a  <  lO'^O  cm^).  Our  calculations  of  the  low-energy  structure  of  the  SiH4~[2T2l 
resonance  suggest  that  autoionization  is  strongly  formed  over  DA  since  the  resonance  state 
geometrj'  is  very  close  to  that  of  neutral  SiH4.  Quantum  calculations  of  the  cross  sections  for  low 
energy  DA  of  e  +  SiH4  are  currently  in  progress  [10]. 

In  contrast,  if  the  resonance  state  is  long-lived,  a  statistical  breakup  may  occur  with  nearly 
equal  probabilities  (a  priori)  for  all  accessible  product  channels.  This  alternate  kinetic  route  can  be 
derived  from  the  assumption  of  LTE  (or  at  least  quasi-steady  state).  In  this  case,  electron 
attachment  can  occur  not  only  to  SiH4,  but  also  directly  to  the  neutral  silane  fragments  which  result 
from  the  decompositions  of  the  negative  ion  resonance  state.  The  product  reaction  channels  are 
most  likel>  those  correlating  to  repulsive  Rydberg  or  high  lying  valence  states  of  SiH4.  In  Table  2, 
we  illustrate  electron  attachment  to  SiH3,  SiH2  and  SiH.  It  can  be  seen  that  a  wide  spectrum  of 
threshold  energies  is  spanned  by  these  molecules  (0.39  -  2.53  eV)  and  that  all  possible  negative  ion 
species,  including  H~,  should  be  possible.  The  relative  importance  of  these  several  processes  is, 
of  course,  dependent  on  the  branching  probabilities  for  decomposition  of  the  resonance  state  and 
direct  attachment  process  to  the  various  neutral  silane  fragments.  Knowledge  of  the  neutral 
composition  of  the  plasma  in  tfie  sheath  region  now  becomes  a  requirement  for  analyzing  the 
relative  importance  of  direct  versus  statistical  breakup  electron  attachment. 

References 

[1]  V.  H.  Ebinghaus,  K.  Kraus,  W.  Miiller-Duysing  and  H.  Neuert,  Z.  Naturforschg.  19a 
(1964)  732. 

[2]  P.  Potzinger  and  F.  W.  Lampe,  J.  Phys.  Chem.  73  (1969)  3912. 

[3]  S.  K.  Srivastava,  private  communication. 

[4]  J.  N.  Bardsley,  A.  Herzenberg  and  F.  Mandel,  Proc.  Phys.  Soc.  (London)  89  (1966)  305. 

[5]  N.  Bohr,  Natur  137  (1936)  344. 

[6]  G.  Breit  and  E.  Wigner,  Phys.  Rev.  49  (1936)  519. 

[7]  J.  N.  Bardsley  and  F.  Nandi,  Rep.  Prog.  Phys.  31  (1968)  471. 

[8]  I.  S.  Elets  and  A.  K.  Kazanskii,  Sov.  J.  Chem.  Phys.  2  (1985)  1439. 

[9]  P.  Haaland,  Ion  Kinetics  in  Silane  Plasmas,  AFWAL-TR-88-2043,  AFWAL,  Kinland  AFB, 
NM  April,  1988. 

[10]  J.  M.  Wadehra  and  H.  H.  Michels,  in  preparation. 


r 


C-5 


Table  1 .  Thermodynamics  of  Negative  Ion  Formation  in  Silane. 


e  -r  SiH4[iAi] 


Reaction  AH  (eV) 


— ^ 

SiH4-[-T2l 

<4.0 

(1) 

SiH2-[2Bi] 

+  H2[^5:g+] 

41.266 

O) 

SiH3-[lA,l 

+  H[2S] 

42.442 

(3) 

-4 

Si-[4S] 

+  2H2[‘Sg+] 

42.783 

(4) 

—4 

SiH3[2A,] 

4  H-[*S1 

43.094 

(5) 

-4 

Si[3p] 

-  H2[Jrg+] 

4  H[2S] 

4  H-[IS] 

43.269 

(6) 

SiH-[3Z-] 

4  H2[lSg+J  4  H[2S] 

44.392 

(7) 

SiH[2ni 

4  H2[*5:g+]  4  H-[»S] 

44.915 

(8) 

afw3l22.1.5/S9 
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Table  2.  Thermodynamics  of  Negative  Ion  Formation  in  SiH/SiH2/SiH3 


Reaction  AH  (eV) 


+  SiH3[2Ai] 

SiH3-[lAi] 

-1.406 

SiH-[3Z-] 

+ 

H2['2g+] 

+0.544 

-4 

SiH2-[2Bi] 

+ 

HpS] 

+1.896 

— ) 

SiH2[lAi] 

+ 

H-(1S] 

+2.266 

“4 

SiH[2n] 

+ 

HpS]  + 

+5.545 

+  SiH2[lAi] 

— > 

SiH2-pBi] 

-1.124 

-4 

Si-[4S] 

+ 

+0.393 

-> 

SiH-[3L-] 

+ 

HPS] 

+2.002 

-4 

SiH[2n] 

+ 

H-IIS] 

+2.525 

-4 

Si[3p] 

+ 

HpS)  + 

H-[1S] 

+5.502 

+  SiH[2n] 

-> 

SiH-[3E-] 

-1.277 

-4 

Si-[4S] 

+ 

HpS] 

+1.592 

-4 

Si[3p] 

+ 

H-[1S] 

+2.223 

afwal22-2.5/89 
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